gr 


ur 
r 
ö 
2 3 
ä 
% 
; 


ZEITSCHRIFT FÜR ANGEWANDTE 
MATHEMATIK UND MECHANIK 


INGENIEURWISSENSCHAFTLICHE FORSCHUNGSARBEITEN 


Band 39 Mai/Juni Heft 5/6 


Einige analytische Ergebnisse über die Wasserbewegung 
in einem untiefen Meere®”) 
Von D. van Dantzig 


Die freundliche Einladung des örtlichen Tagungsleiters dieser Versammlung, einen Haupt- 

vortrag zu halten über einige in den letzten Jahren in den Niederlanden durchgeführte For- 
schungen auf dem Gebiete der angewandten Mathematik, bietet mir — obwohl ich eigentlich kein 
angewandter Mathematiker bin, sondern ein Mathematiker der nur gelegentlich angewandt wird 
— die Gelegenheit, Ihnen einen, wenn auch sehr unvollständigen, Überblick zu geben über einige 
Arbeiten, die unter meiner Leitung seit 1953 in der Abteilung für angewandte Mathematik des 
Mathematischen Zentrums in Amsterdam vollbracht wurden. Sie wurden veranlaßt durch die 
Sturmflut, die unser Land am 1. Februar 1953 getroffen hat, beruhen auf einem Auftrag des 
Centrale Studiedienst van de Waterstaat, und bezwecken, einiges Verständnis zu gewinnen für 
die Weise, in der einer über ein Meer, insbesondere die Nordsee, sich bewegender Sturm die Höhe des 
Meeresniveaus beeinflußt. Obwohl die endgültige Lösung des Problems noch in unabsehbarer 
Ferne liegt, möchten unsere Bemühungen, die sich gezwungenermaßen nur auf wesentlich ein- 
fachere Probleme beziehen, vielleicht doch einiges Interesse erwecken können. 
Bis auf einige wenige durch die Umstände erzwungene Ausnahmen sind die Ergebnisse 
nicht mir selbst, sondern meinen Mitarbeitern zu verdanken. Dies waren anfangs Herr T. C. 
Braakman und Herr G. W. Veltkamp, später Herr Dr. H. A. Lauwerier, und jetzt auch 
Herr Dr. D. J. Hofsommer. Wenn auch gelegentlich Anregungen, Ansätze und weitere Sugge- 
stionen meinerseits in den Berichten verarbeitet wurden, so haben doch meine Mitarbeiter zum 
größten Teil ihre Ergebnisse selbständig erhalten. 

Zeitbeschränkung nötigt mich, mich auf einen Teil — wenn auch einen wichtigen Teil — der 
Arbeiten zu beschränken, alle Ansätze zu den Eigenwertproblemen außer Betracht zu lassen, und 
auch über die zu behandelnden Arbeiten nur einen manchmal oberflächlichen und sehr unvoll- 
ständigen Überblick zu geben. 


1. Formulierung und erste Reduktion des Problems 


Die meisten Arbeiten, die seit etwa 5 Jahren in der Abteilung für angewandte Mathematik 
des Mathematischen Zentrums in Amsterdam unter meiner Leitung ausgeführt wurden, beziehen 
sich auf die Bewegung des Wassers in einem rotierenden Becken unter Einfluß von gegebenen, von 
einem Windfelde verursachten Kräften an der Oberfläche, sowie von Bodenreibung. 

Ein Nachteil der analytischen Behandlungsweise besteht darin, daß man erhebliche, und 
im ‚Grunde nicht gerechtfertigte Vereinfachungen anzubringen genötigt ist. Als Modell für die 
Wasserbewegung in einem Meer unter Einfluß von Gravitation, Wind und Reibung betrachten 
wir diejenige in einem mit konstanter Winkelgeschwindigkeit rotierenden Bassin von konstanter 
Tiefe. Die nichtlinearen Bewegungsgleichungen werden durch die entsprechenden linearisierten 
ersetzt, d. h. wir betrachten in erster Näherung kleine Abweichungen vom Gleichgewichtszustand. 
In der linearen Annäherung (nicht in Wirklichkeit) kann deshalb die Gezeitenbewegung außer 
Betracht gelassen werden, weil die berechnete Bewegung auf sie superponiert wird. Der Einfluß 
des Windes wird als eine gegebene Schubspannung!) an der Oberfläche betrachtet; die Reibung 
(insbesondere die Bodenreibung) wird durch eine der lokalen Geschwindigkeit proportionale 
Massenkraft ersetzt. Entlang einer jeden Vertikalen wird eine quasistatische Druckverteilung 
angenommen, und die Geschwindigkeitskomponenten werden durch ihre Integrale über eine 
solche Vertikale ersetzt. Schließlich werden als Randbedingungen nur undurchdringbare Küsten 
oder Begrenzungen durch Ozeane in Betracht gezogen, welche letztere als unendlich tief ange- 
sehen werden. 2 

Wir stellen das Bassin durch das ebene Gebiet D, seinen Rand durch Cdar, D+ C durch D, 
die über die Vertikale integrierten Geschwindigkeitskomponenten im Punkte (, y) durch 


*) Hauptvortrag, gehalten auf der Wissenschaftlichen Jahrestagung der Gesellschaft für Angewandte 


Mathematik und Mechanik vom 8. bis 12. April 1958 in Saarbrücken. - 
1) Das barometrische Feld kann als ein Gradientfeld in der äußeren Kraft aufgenommen werden, und 


wird deswegen nicht ausdrücklich erwähnt. 
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2= (u,, u,), die Erhebung des Wasserspiegels über das Gleichgewichtsniveau durch w, die Kompo- 


nenten der durch den Wind verursachten Kraft durch X = (X, X,), den Reibungskoeffizienten 
durch A, den Coriolis-Parameter (= 2 wsin ß) durch 2, die Tiefe durch A, die Gravitationsbe- 
schleunigung durch g, die stationäre Wellengeschwindigkeit Yghdurchc und die Dichte durch 1. 
Die linearisierten Bewegungsgleichungen und die Kontinuitätsgleichung lauten dann be- 
kanntlich (H. Lamb, 1932): un 
— +41 —- 2ı%+ ex, 2uw+ = Haut = = + 2-0 (1.1) 
oder in Vektorform 
ow 
ie 
wo Q ein senkrecht nach oben gerichteter Vektor der Länge 2 ist. Entlang einer Küste C, ver- 
schwindet die normale Geschwindigkeitskomponente u,„, entlang eines Ozeanrandes C, die 


Erhebung w. 
Wir nehmen eine Laplace-Transformation vor, oder, was auf dasselbe hinauskommt, wenn 
wir die alten Bezeichnungen für die Laplace-Transformierten beibehalten, wir ersetzen den 


du =, P7 > > R > 
Zy rt 2 xuregradu=X, div ır+ 0, 


Differentialoperator & durch einen Faktor p. Selbstverständlich wird dadurch eine der wesent- 


lichen Schwierigkeiten bis zum Schluß hinausgeschoben, nämlich die Umkehrung der Laplace- 
Transformation (Vgl. 3.). 


Die Determinante der Koeffizient-Operatoren der linken Glieder von (1.1) ist, mit p = = 


ot 
ae 
p 0X 
0 
2 p+4 Re = 2(p+9)A— vp{P +90... .. (12) 
ro & 
0x oy p 
def 9? 0579) SER. £ s 
od Fr + a Elimination von u, und u, ergibt also für w eine Gleichung 
a EEE (1.3), 
wo 
def 
am HEFTE I FEINE Re da 
‚def IX oX 2” 10X oX er. > 
F= a h rs D le = div Kpitgyirot X SEES) 
mit 
def Q 
| Der (1.6). 
Die Randbedingungen lauten dann: 
w=( entlang eines Ozeans 
ow 1) Br | 
an +1tgy , — le entlang einer Küste 2) 
mit 
Tex, ey ET BER BT, 


wo X,, X, die normale und die tangentiale Komponente des Vektors X darstellen, erstere beim 
Durchlaufen des äußeren Randes von D im positiven Sinne nach außen gerichtet. 


Die Geschwindigkeitskomponenten U), U,, welche Gleichungen derselben Form (1.3) (mit 
demselben x*, aber mit anderen rechten Gliedern) genügen, ergeben sich sofort aus ı: 


+++ (ne) + ol, 2). 
. (1.9) 


u, = {(p +2)? + O9. I-2(& Pe = +(p+A eo 4) 


def 2 kr s 2 
?) Das Symbol = bezeichnet eine Gleichheit, die das linke Glied durch das rechte definiert. 
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oder in Vektorform 


a = {(p +92 + 99 (p+4+2x)X — @gradu). 


Wir bemerken noch, daß, bei konstantem (p + 4)? + Q2, ı, = = ‚u=-— OU, ist, 
Y 

Die größten Schwierigkeiten werden von der „schiefen‘‘ Randbedingung (1.7), also von der 

Corioliskraft verursacht, die weiteren von der komplizierten Form (1.4) von x2, das für 2=/=0 

. einfach p?/c? wäre. Man möchte deswegen versuchen, eine Störungsrechnung hinsichtlich 2 vor- 

zunehmen. Dies geht aber nur für viel kleinere Gebiete als die Nordsee. Im allgemeinen kann 


man sagen, daß anfänglich, wenn ein Sturm einsetzt, nur die äußeren Kräfte X von Bedeutung 
sind; nach einiger Zeit beginnt sich der Einfluß von 2 zu zeigen, und noch später macht sich die 
Reibungsdämpfung bemerkbar. Auf dem Höhepunkt eines Orkans läßt sich keiner dieser Fak- 
toren vernachlässigen. Bei einem einigermaßen ‘lange anhaltenden Sturm wäre eher p im 
Vergleich zu 2 (nicht aber zu A) als 2 gegenüber p vernachlässigbar. Auch macht es die Coriolis- 
kraft unmöglich, die Wasserbewegung eindimensional zu behandeln, sogar wenn man die Nordsee 
alseinen unendlich langen Kanal betrachtet; dieQuerbewegung ist dafür von zu großem Einfluß. 


2. Zurückführung auf Fundamentallösungen 


Zur Diskussion der Lösungsmethoden der Differentialgleichung (1.3) mit den Randbe- 
dingungen (1.7) benutzen wir den Greenschen Satz. Einfachheitshalber setzen wir das Gebiet D 
als einfach zusammenhängend voraus, und sein Rand C sei stückweise glatt und bestehe aus zwei 
Teilen C, (Ozeanrand), C, (Küste) die je aus endlich vielen Bögen bestehen. Weder die Konstanz 
von x? aus (1.4), noch die stückweise Konstanz von 1(s) = tgy aus (1.7) auf C,, ist für das nächst- 
folgende nötig. Punkte (x, y), (%,, Yo) usw. seien kurz mit P, P,, usw. angedeutet. Die Entfernung 


Ya, 2° + y,)? zweier Punkte P, Pzs sei mit o(P,, P,) bezeichnet. Es sei G,(D) die 
Menge aller Funktionen in D, die dort zweimal stetig differenzierbar und am Rande stückweise 
glatt sind, und S(D, M) die Menge aller Funktionen, wo dies für das Komplement bzgl. D einer 
endlichen Menge M eD gilt, und die in jedem ce M wie — (2 n)-" In o(P, c) unendlich werden. 

Für je zwei Funktionen ge S(D, M,), ye S(D, M;), wo M, NM, = 0 sei, gilt dann für die 
„Dirichleische Entfernung“ (S. Bergmann und M. Schiffer, 1953) 


ar ( (po pop, , 
2 De HE He pp)dady 


0) 
-— [pw #yaray— | pH as+ 5 9 GE (2.1). 
"D G on aeM, 


Weil E{p-y} offenbar symmetrisch ist, erhalten wir für den Fall, wo M, nur aus einem Punkte 
P, besteht: 


Be h A 
y(P,) - | [var naar | | p (Ay — x: Day + | vn | P2,ds 
D [6] 


o &% 


MT. 0 
a ‚vie vl) | El + as + |* Womit (2) 
Cr 
Es sei jetzt WB) —= w(P) eine Lösung e C,(D) von (1.3), (1.7) na Y(P) = G(P, P) eine beliebige 
Funktion € S(D, z,). Dann wird (2.2), unter der Voraussetzung, daß 9 y auf ganz C stetig und 
stück weise glatt ist?): 


N Ypenns! | GFaxay+ | Gfas+ | | w (AG — 26) de dy — [o®as 
D- E ei } 
ee) ; 
— a ul 2 7 0,3; 
+| w(> 7) G)ds (2.3) 


C% 
A. Falls eine verallgemeinerte Greensche Funktion G gefunden worden ist, d.h. daB G 
a) der homogenen Differentialgleichung 
EEE 0 WED Wi en = Eee), 
b) der Ozeanbedingung 
> Sl) Hanno) [OR Fo . (2.4b), 


s) Also, wenn @ und y stetig sind auf C, muß in jeder etwaigen Unstetigkeitsstelle von r(s) entweder 
y oder y verschwinden. 
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c) der homogenen adjungierten Küstenbedingung 


86.79 
ee = ET, ne EEE EEE 
pe )G=0 auf G 


genügt, ergibt (2.3) eine explizite Darstellung von w(P,) mittels Quadraturen. 

B. Falls aber eine Funktion G bekannt ist, die wohl der Differentialgleichung (2.4a) und 
der Ozeanbedingung (2.4b), aber nicht der Küstenbedingung (2.4c) genügt, so stellt (2.3) eine 
für Pe C; singuläre Integralgleichung für w(P) entlang der Küste C, vor. Falls @ überdies der 
Küstenbedingung auf einem Teile von C, genügt, so bleibt eine ebensolche Gleichung für den 
Rest von C, übrig. Diese singulären Integralgleichungen können mittels bekannter Methoden, 
die zum Beispiel im Buche Muskhelishvili’s über singuläre Integralgleichungen dargestellt 
worden sind, durch Ausübung des adjungierten Integraloperators in reguläre verwandelt und 
sodann nach klassischen Methoden aufgelöst werden. y 

C. Man kann aber auch sofort eine reguläre Integralgleichung erhalten, wenn auch eine zwei- 
dimensionale, falls man eine Funktion G=G, bestimmen kann, die den sämtlichen Rand- 
bedingungen (1.13b, c), aber nicht der Differentialgleichung (1.13a) genügt, jedoch stattdessen 
harmonisch ist (AG, = 0). Für den Fall wo C, (also auch C;) ein zusammenhängender Bogen ist, hat 
H. A. Lauwerier (TW31, 1955) eine solche Funktion G, konstruiert mittels einer konformen Ab- 
bildung von D auf eine Halbebene, bei der C, auf die negative reelle Achse I an wird. 

i 

Man kann G, durch eine Summe von vier Gliedern von der Form + ge H ) 
darstellen, wo H(w) eine unvollständige Beta-Funktion ist. Für den Spezialfall eines Rechtecks 
von der Breite 2 a und der Länge 2 b, bei dem C, eine Rechtecksseite ist, und das ein Modell der 
Nordsee darstellt, ist M(z) = (2 — 2a —2ib) — p(2ib). 

Es ergibt sich dann für w unter Vertauschung von P und P, aus (1.12) 


w(P) = — x? Sf wP) GP, P,) dx, dyo a Go(P, Po) Fo) N G,(P, P,) F(z,) dx, dy, (2.5), 
k 


aus welcher sich w mittels der Neumannschen Reihe als eine Potenzreihe nach Potenzen von 
»2 ergibt. 

Weil dabei der Koeffizient von x?” u.a. noch ein 2n-faches Integral enthält, worin der 
Integrand a) vom Aufpunkte P, b) von den Randbedingungen, c) von der Form des Gebietes, 
d) von dem gegeben äußeren Felde abhängt, ist es von Bedeutung, die Fälle zu betrachten, in 
denen sich die Greensche Funktion explizit bestimmen, und zwar mittels einer beschränkten 
Anzahl von Integrationen herstellen läßt. Dies ist gelungen für diejenigen Gebiete, die nur von 
höchstens zwei Geraden begrenzt werden. 

Die zu bestimmende Greensche Funktion sei, unter Vertauschung der Punkte P und P,, 
mit G(P, P,) bezeichnet. Sie genügt dann hinsichtlich P und P, der Differentialgleichung 
(2.4a) als der Ozeanbedingung (2.4b). Der Küstenbedingung (2.4c) genügt sie nur hinsichtlich 
P,, während mit Bezug auf P die homogenisierte Bedingung (1.7), d. h. allgemeiner 


0G 0%G 
a + 7) Fa) aut“ Ca 5 20 ee De 


gilt. Die Funktion G(P, P,) ist also nicht symmetrisch. In unsrem Fall, wo z(s) = tgy auf der 
Küste konstant ist, geht sie in sich über, falls die beiden Argumentpunkte vertauscht werden, und 
gleichzeitig das Vorzeichen von tgy, d.h. von 2 umgekehrt wird. 


3. Bestimmung von Fundamentallösungen für von höchstens zwei Geraden begrenzten Meeren 


Ein wesentlicher Beitrag zur Bestimmung der Fundamentallösungen der Differential- 
gleichung 


AC-HESO,.. RUE TERBERTEN 


mit homogenen schiefen Randbedingungen (vgl. 1.7) wurde von G. W. Veltkamp (1955) in 
einer kurzen internen Note für das Mathematische Zentrum gegeben, über deren Grundgedanken 
ich jetzt kurz berichten werde. 


Die Fundamentallösung G(P, P,), die in der ganzen Ebene der Gleichung (3.1) genügt, 
im gegebenen Punkte P, = (x, y,) eine Singularität G = — = In o(P, P,) + O(1) aufweist und 
7 
für P— P, ins Unendliche verschwindet, ist gegeben durch die Besselsche Funktion 


+ © 


1 
5, Ku) = z- | erwechäiggniic, 1 te ale log I (3,2). 


er, 


j 
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Wir werden sie, ähnlich wie in der Theorie der harmonischen Funktionen, kurz einen 
Pol (in P,) nennen. Durch Differentiation nach dem „Aufpunkte“ P, in einer vorgegebenen 
Richtung ergibt sie einen Dipol, durch zweifache Differentiation in zwei Richtungen einen 
Quadrupol, durch k-fache Differentiation einen 2*-Pol. Wir betrachten den Richtungsdifferen- 
tialoperator 


. (8.3). 


In Polarkoordinaten ist 


1 
fallsytpo=a+ 57 nr (mod r) ist. 
Die Operatoren D, kommutieren für alle Werte von &, und sie erfüllen, wie man leicht 
einsieht, bei Beschränkung ihres Anwendungsbereiches auf Lösungen der Gleichung 
Au—Ru=0, 
die Identität 


DD 5, D, = uno —=d)sina an Fi), 
oder auch 
De D, Bra Dia sin? DE 
>. 
Bürz = Ost Days. Dy — Des mitin Ska) &, so daß es für jedes k bis auf Bewegungen nur 
1 


eine Art von harmonischen 2*-Polen gibt; für x = 0 aber ist das nicht der Fall. Zufolge (3.5) 
. ist die schiefe Randbedingung D,;f = 0 entlang der x-Achse erfüllt, falls {= D_;g und g als 
Funktion von y ungerade ist. Es ist dann nämlich g = 0 für y = 0, folglich auch D5g = N 0 
entlang der x-Achse, also nach (3.5) ebendort D,D,9 = 0. 

Um eine Fundamentallösung von (3.1) in der Halbebene y > Omit Pol in 2, = (%, Yo) 
zu erhalten, die der Randbedingung 


DiE=V et Daran een rer 
genügt, setzen wir 0* = {(e — 2)? + (y + yo}? und 279 = Kyu(l# 0) — K,(x o*). Dann ist g 
ungerade in y; folglich genügt 
1 


Les; DE Kl DR AO a ee) 
sowohl (3.1) als (3.6). Diese Lösung besteht aus Fe, 
zwei entgegengesetzten parallelen Dipolen in den ER 
Spiegelpunkten z, und 2, = (&, — Yo), deren Fe 


Richtung spiegelbildlich zu & ist (vgl. Bild 1). 

Weil diese Lösung für jeden Aufpunkt z, der 

Bedingung (3.6) genügt, bleibt dasselbe der Fall, Br ws 
wenn man z, durch (@, — cos, yy + £sin a) für 

allei>0 ersetzt und über £ integriert. Dabei 

geht der Dipol in z, in einen einfachen Pol über, ki 

während der Dipol in z, in eine mit gleichen und \ 

parallelen Dipolen gleichmäßig besetzte Halb- Rn Ay 

gerade in der Richtung & übergeht. Zerlegt man an 

diese Dipole in ebensolche senkrecht und parallel 

zu dieser Halbgeraden, so erhalten erstgenannte die Größe sin 2a, während letztere sich zu 
einem einzigen Pol in z, der Stärke — cos 2% zusammensetzen (Bild 2, Bild 3). 


Analytisch wird die Lösung durch 


+© 


Iapipns We ee IM 
En, ar We — Le r(y+Y) cht—ixr (&—%o) sh dt, 
ee sh (t— io)“ 


3,5 


+@ 
1 1 e 
Fir 42 —x |y—yo| eht—ir (e—x,) shit 
an Kb re J : . | 


—o 


gegeben. 
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Dieses Ergebnis, das Veltkamp auch für beliebige Richtung der begrenzenden Geraden 
herleitete, wurde auf etwas anderem Wege auch von H. A. Lauwerier bewiesen und zur Lö- 


IS „+ 


—————, AR yeo AR 20 


+ PA; x 
RR DASWEIT TER 
Ni N 
> NN- 
an & > Be 
Bild 2. (Die beim oberen Halbstrahl eingezeichneten + - Bild 3 


und — -Zeichen sind zu streichen. Dieser Halbstrahl ist aus 
„Kopf zu Schwanz“ liegenden Dipolen zusammengesetzt.) 


sung des Randwertproblems für einen (beiderseits unendlichen) Parallelstreifen ausgenutzt 
(TW 31). Für den Streifen |r| < a (unendlich langer Kanal) mit 


0 eos y +5, siny = 01. fs za 
wird 
ıl il u 5 ch (n + iy) 
ur pin 87 | | —+(2a+2,—e)chn »(2a—2,+2)chn ——l 20%, AS ACH 
pe ot 2 e te ea 
ch (n et! 2) —x(z+7 nn] ei»Ww—yo)shn 
ch (m + ip) sh(@wachn) 1 et) 
wofür auch 
+%© 
1 © eben ER 
types erh d | ern 2najxchh, Alm + isgn n —-1)"y)| 
x . ei ch(n+isgnny) | (3.10) 


geschrieben werden kann. (Vgl. Bild 4.) 


Te 


ı +1 } + EN \ 


Bild 4 


: ch(n Fiy) , : a ih \ 
Die Faktoren Po (7 Fan ändern einen Pol in einen „Dipolschwanz‘‘. Für den Streifen 
0<y<bmit 


DG=-I ru, G=0 für y=b, 


zug 1 
so daß y=- 0 eine Küste (a Ser y) und y= b eine Özeanbegrenzung darstellt, erhält 
Lauwerier die Lösung ” 


+ © 
1 
G= 5 Kuo) — % Kolo*)+2 | dn”® nn ı uZychn 
An am 1 shm—iy) e?xachn 


. e-ix(e—x,)shn 


{a 


| 3 | sh(n+iy) (3.11), 
deren Schema in Bild 5 skizziert worden ist, 


NE 
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Wesentlich schwieriger gestaltet sich das Problem schon im einfachsten Falle, wo der 
Rand eine Ecke aufweist, nämlich im Falle, wo D ein Winkel ist, z.B. in Polarkoordinaten: 
r>0, 91 <9P<g, Wir werden da die Lösungen, insbesondere die Fundamentallösung der 
Gleichung (3.1) mit den Randbedingungen | 


er 0G £ 0G 0G n 
en, it o=9, en für =  (.12) 


betrachten, obwohl wir eigentlich nur die beiden Fälle.y, = y, (zwei Küsten) und coty, = 0 
oder coty; = 0 (eine Küste und ein Ozeanrand) (im letzteren Falle mit der Zusatzbedingung, 
daß dort w = 0 sei) bedürfen. Schon 1956 hat 
G. W. Veltkamp (TW 40) das Verhalten einer für 
r—0 endlich bleibenden Funktion w in der Umge- 
bung der Ecke studiert, die in D definiert ist und 
entweder der Potentialgleichung Aw=0 genügt (oder 
allgemeiner Aw = O(r”2) für r—0) und den Rand- N 
bedingungen (3.12) (mit w anstatt G; im Falle y, =; 
ist r2 durch’ re, 5 >0'zu ersetzen), oder für 
x” =0 den speziellen Randbedingungen (3.12) mit 
yı = ya Oder cos y, COS yy = 0 genügt. Die Notwendig- 
keit, mehrere verschiedene Fälle zu unterscheiden, ließ 
wenig Hoffnung, auf diesem Wege durch Heranziehen 
höherer Näherungen weiterzukommen. Beiläufig sei 
bemerkt, daß ein Spezialfall, nämlich y —yı = 9 
(mod x) auf Grund von (3.8) (Bild 1) trivial ist. 


Jüngstens aber ist es uns gelungen, das Problem 
für den Winkel vollständig zu lösen. Herr Lauwerier’s 
Krankheit zwang mich, eine schon 1956 in groben Zü- 
gen skizzierte Methode selbst auszuarbeiten; seit der 
Wiedergewinnung seiner Kräfte hat Lauwerier so- 
dann unabhängig von der meinigen das Problem nach 
einer anderen Methode gelöst, und später mehrere 
wertvolle Zusätze zur meinigen angebracht. Lau- 
werier’s Methode beruht, ähnlich wie es bei den 
Wiener-Hopf-Methoden geschieht, auf der Zurück- 
führung des Problems auf ein Hilbert-Problem. In 
seinem Vortrag (TW 50) behandelt er in dieser Weise 
den speziellen Fall eines Vollwinkels. Übrigens lassen 
sich die wichtigsten Elemente jeder der beiden 
Methoden auch bei der anderen wiederfinden. Ich werde deswegen meine eigene Lösung kurz 
skizzieren. 


! 
8 
| 
is - 
i 
u 
ln 
ö 


Bild 5 


% Ist’G ı0.D rt, o)Ir>0, 9, <gP<Y,} mit Ausnahme höchstens einer endlichen 
Menge von logarithmisch singulären Stellen eine Lösung der Gleichung (3.1) die für r— 
genügend schnell gegen Null geht, so ist 


VER ea Melange.) 
0 


eine harmonische Funktion, also in der Form 


In = Un ti) +, n tip) ..- nee. « (8.14) 
zu schreiben. Dabei sind U, und U, stückweise holomorph für 9, Sp = 9,; auf den Geraden 
p = konst., die den Singularitäten von G entsprechen, können sich U, und U, unstetig ändern. 
Sind andererseits U, und U, stückweise holomorph und ist fürn — + ®, (ch n)® U — 0 für eine 
gewisse Konstante c, so genügt, der Lösung von (3.13), (3.14) entsprechend, 

+ 
13. uf } h ’ 

G(r, 9) Eye | ewehn(U,(n +ip+U,(—-n+ip)}dn..... . 8.15) 
in D der Gleichung (3.13). Die Transformation (3.13) ist im wesentlichen einem schon von 
Peters (A.S. Peters, 1952) für das „sloping beach“ Problem benutzten analog. 

2. Den Randbedingungen (3.12) entsprechen zwei Funktionalgleichungen für U, und U,, 


die eine spezielle Lösung a 
UA (a DE u IP (0 ER 1 (e BEER a 13,0) 
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besitzen, wo D(£) eine in der ganzen Ebene meromorphe Funktion ist. Und zwar ist 


72V N a a un 


wo e(&,y) für |Im£|< 0 + 5 |Re y| durch 


\ +@ 
def 1 dt shyi(l— cosöh) 
e(£, y) “nn | nr ee hot r u ET TE (3.18), 
tsh — rt 
— (60 2 
eye glaub (b WIER Kun TYENBEEE 
definiert, aber in der ganzen Ebene meromorph ist. Mn 


Gemäß (3.15) entspricht ® einer überall in D, außer vielleicht in r = 0 reguläre Lösung 


von (3.1), (3.12). Für — 5 


Für y, Sy, dagegen gibt es keine überall, einschließlich r = 0 reguläre nichttriviale Lö- 
sung.?) 


Sp<qmS 5 ist sie auch für r = 0 regulär.*) 


3. Um die Fundamentallösung zu finden, müssen wir die Funktion K,(z o(P, P,)) in der 
Form (3.15) darstellen: 


es 
| EUGENE A AMEIRT, ERIER L (3.20). 


—-& 


K,= Kolz o(P, P,)) = 


Je nachdem 9 > 9, oder 9 < 9, ist, entspricht also U, oder U, der Null, und U, oder U, der 
Funktion (4)! ei*rsshw+ile—eo), Fordern wir jetzt, daß G— (2r)!K, überall regulär ist. 
Es entstehen U, tınd U, aus ® durch Multiplikation zweier im Streifen g, < Im <s», stück- 
weise holomorphen Funktionen, die periodisch mit der Periode 2i © sind und an den Geraden 
9=9+2mO bzw. 9=29,—9+2mOi (m gerade) einen Sprung (4)! ei*rshn bzw. 
— (4 n)t ei*rsshn aufweisen. Mittels einer Plemelj-Formel lassen sich diese Funktionen dazu 
bestimmen. Das Endergebnis ist 


so 25 

iu gg ar 
lan Er zu ur 
[) 0 12) Ss 


=] | osersunsintenshn 


de a n 
für -Z Sy Ssp=s DE und 
0) e) TT 
sh—=& 
E j ’ ® oO 
Gr,p)=, Pr) sin («rshn)cos (x r,shn,) & d£,dE (3.22) 
We Po) ch &,—ch =» 
on 62 oO 0 ö [@) > 
7 7 n 5 \ 
für — Jen<n 5 wwi=n +19, o=Mm + ip undg, =0 m,= © gewählt ist. 


Gemäß Lauwerier’s Bemerkung ist die Lösung für y, < y, die einzige, die für r—0 
endlich bleibt; für y, > y, kann man noch eine überall in D reguläre Lösung hinzufügen. Da- 
durch kann man erreichen, daß sie für r = 0 Null wird. Die Lösung (3.22) ist dieser Forderung 
schon angepaßt worden. 

Weitere Untersuchungen dieser Lösung sind augenblicklich noch im Gange. Leider ist 
nicht zu erwarten, daß eine entsprechend einfache explizite Lösung für beliebige Form des Ge- 


bietes D anzugeben ist. 
4. Zeitabhängige Lösungen 


In kurzer Form möchte ich noch einiges berichten über die schon anfangs erwähnten 
Schwierigkeiten und über die Möglichkeiten bei der Umkehrung der Laplace-Transformation. 


*) Ich verdanke diese Bemerkung Herrn Lauwerier, 


. 
& 
5 
3 
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Herr Lauwerier hat dieses Problem in Angriff genommen, und wir verdanken ihm mehrere 


. recht verschiedene, wenn auch immer noch einfache Fälle vollständiger Lösungen, manchmal 


bis zu numerischen Resultaten durchgerechnet. 


Zunächst (TW 35) betrachten wir den Fall, wo D die ganze unendliche Ebene ist, wo also 
weder Küsten noch Ozeane vorhanden sind. Ersetzen wir in 


+%© 


n durch n + In Bi und setzen wir für x? den Ausdruck (1.4) ein, so ergibt sich 
+ oo 


K,(& r) = 5 J e 


—o 


al au 1 y22 
en NE Mr 
=D 2 r(p-+A)e Sen, e 


ae Ei Saat (4.2). 


Zur Lösung kann man p mit dem Operator identifizieren und sodann (4.2) auf eine 


& 
gegebene Funktion der Zeit X(f) (die auch von den räumlichen Koordinaten Bon kann) 
wirken lassen. Im einfachsten Falle, wo weder Reibung noch Drehung vorhanden it A= 2 = 0), 


ist x = p, und (4.1) ergibt 


oo 


Kar) X) = y xumrenan= | Xn,e 


r 


. (4.3). 


Yar — 


Falls wohl die Reibung, aber noch nicht die Erddrehung in Betracht gezogen wird, erhalten wir 


RR 
Kulx r) X) -|« r ch (4 year) X DIAS HEN ER (4.4) 
NE 078 
und endlich im allgemeinen Fall 
[) a dr 
KR, T) &o-= fie ch (5 Ayı2 — D X(t—?)- 
er 
rt af 
[e0) 2 44 r2 
RS e*dıX(t—r) | dp J,(Q c0sp-Yr?— 72). exp— as Fate (4.5), 
% % 
2 T+sing- YRa—r2 


r ee 


ein Ausdruck, der noch verschiedentlich umgeformt werden kann. Lauwerier betrachtet 


t 
speziell den Fall, wo, bis auf eine zeitlich konstante Funktion, X(f) = «!) = h vr 


d.h. wo ein Windfeld plötzlich einsetzt und dann konstant bleibt. Insbesondere hat er die 
folgenden Spezialfälle untersucht: 


ist, 


a) Kein Wind; rotationssymmetrisch verteilte Depression. Es handelt sich dann lediglich 


um ein Barometergefälle, d.h. der Vektor X in 1 ist ein Gradient. _ 
b) Rotationssymmetrisches wirbelfreies Feld. (Auch hier ist X ein reiner Gradient.) 
c) Rotationssymmetrisches Wirbelfeld. 
d) Punktquellenfeld konstanter Richtung. 
e) Beliebiges Windfeld konstanter Richtung, das bis auf einen beliebigen zeitabhängigen 
Faktor, eine zeitlich konstante rotationssymmetrische, aber sonst beliebige Verteilung hat. 


Schon vorher (TW 32) hatte Lauwerier für den Fall einer Halbebene ein räumlich kon- 


stantes Windfeld untersucht, das von der Zeit mittels eines Faktors X(f) Z mtem: (>0) 
abhängt, also anschwillt, ein Maximum erreicht und wieder abflaut. Die Erhebung w läßt sich 
mittels elementarer Methoden nur an der Küste berechnen, und wird dort dargestellt durch 
die Konvolution der beiden Funktionen 


eek nd. a a, 
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und 
t 


Ar 1 Fe | 1 { 
v)=e ? n(521)—2e 2 E J(@ 9) (514-9) de... 
e 0 


Die numerischen Rechnungen hierzu wurden von der Rechenabteilung des Mathematischen 
Zentrums ausgeführt und gestatten eine gute Übersicht über die Weise, in der die Erhebung w 
vom Winkel zwischen Windfeld und Küste abhängt. RL 

Es ergibt sich, daß die maximale Erhebung für etwa « = 170° maximal ist, so daß die 
Erhebung an der Küste viel stärker von einem zur Küste parallelen als von einem senkrecht 
einfallenden Wind beeinflußt wird, eine qualitativ von vornherein einleuchtende Tatsache. 
Nachdem w(f) sein Maximum erreicht hat, dämpfte es allmählich in leicht oszillierender Weise 


ab; die Oszillationen haben eine Periode von etwa = und korrespondieren analytisch mit den 


Singularitäten der Laplace-Transformierten bi p=—A+tiQ2. Das Erhebungsmaximum 

kommt viel später als das Windmaximum. Das Erhebungsmaximum hängt der Größe nach 

kaum von Q ab, in seiner Abhängigkeit von & aber (offensichtlich) sehr stark. Auf die wichtigen 

Untersuchungen Lauweriers bezüglich eines streifenförmigen Meeres („unendlich breite Nord- 

see‘‘) einzugehen, sowie auf verschiedene spätere Ergebnisse, verbietet mir leider die Zeit. ü 
Ich möchte noch einige Bemerkungen über die Inversion der Laplace-Transformation 

im Falle des Streifens, also der Lösung (3.10) hinzufügen. Der exponentielle Faktor unter dem 


Integralzeichen läßt sich ganz ähnlich wie oben bei K,(z r) durch Ersetzung von „7 durch „ + In = 

behandeln. Weiter ist P 
ch(n-+ iy) 2i@thn 

u ——=l+ — ——— nn ..2.2:.2....0. 0. (48). 

ch (n—iy) T nd ihn un, 

Hierin aber muß noch die obengenannte Verschiebung von n durchgeführt werden, wodurch 


th) in th (m + In = übergeht. Nun wird 


chm+ip) _P+A—iQ er(p+i+i9®+pP+M 

hm—iy) "p+i+iQ erp+r—iN®+pP+A) 
Zerlegung in Partialbrüchen ergibt eine ziemlich unangenehme Abhängigkeit von n, da e?” unter 
dem Wurzelzeichen vorkommt. 


Man muß aber bedenken, daß bei den wirklichen Anwendungen / nur sehr ungenau be- 
kannt ist. Die Berechnung der Reibung beruht nämlich auf Turbulenztheorien, die nicht alle 
dasselbe Resultat geben (Schönfeld 1954). 


Überdies ist die Reibung proportional zum Strom und unabhängig vom Winde angesetzt, 
was in Wirklichkeit nicht der Fall ist. Schließlich ist A abhängig von, und zwar näherungsweise 
umgekehrt proportional mit der Tiefe. Im Falle der Nordsee nun ist die Tiefe keineswegs kon- ° 
stant. Sie nimmt gen Norden allmählich zu und ist in der Mitte durch die Doggersbank stark 
beeinträchtigt. Es ist also der Wert von A nicht nur nicht genau bekannt, sondern nicht mal 
genau definiert, ja sogar nicht genau definierbar.5) 

Nun hat es in der angewandten Mathematik nicht den mindesten Sinn, im rein mathe- 
matischen Teil vollkommene Exaktheit zu verlangen, wenn man bei der Anwendung auf sehr 
unexakt bestimmte Größen stößt. Man wird also im letzten Glied der Formel (4.9) (die ich 
lediglich als Beispiel gewählt habe) bedenken, daß Vermehrung des Zählers sowie des Nenners 


ee 


mit 7 42 lediglich auf eine kleine Änderung des Wertes von A hinauskommt, also infolge der 


geringen Bestimmtheit von / „dasselbe“ Resultat ergibt, vorausgesetzt, daß lediglich große 
Werte von p einen wesentlichen Beitrag geben, das heißt, daß man nur kurz dauernde Stürme 
in Betracht zieht. Damit geht (4.9) über in 
ch(n+1y) | 210 G C 
Y . 4 1 > = E 2 h 1 4 1 + 2 
ch (n—iy) p+A+iQ 2 ‘ A 
rider) Rail] 


| | 4.10), 
wo C, und (C, nicht von p abhängen. 


| 5) Übrigens gilt Entsprechendes auch, wenn auch in geringerem Maße, für die Fortpflanzungsgeschwindig- 
keit c= Vs h der ungestörten Wellen, die wir sogar als Einheit der Geschwindigkeit gewählt haben. 
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Welche Möglichkeiten unsere Lösung des Problems vom Winkel bieten wird, läßt sich 
jetzt noch nicht sagen. Einige Bemerkungen aber kann man schon machen. 


Die Nordsee hat bekanntlich in grober Annäherung die Form eines Rechtecks. Für die 
holländische Küste aber werden bei den häufigst vorkommenden Stürmen weder die Ostseite 
(Deutsche Bucht, Dänemark, Kattegat, Norwegen), noch die Nordseite (Ozean „rand‘‘) großen 
Einfluß haben. Ein rechter Winkel ist deswegen als Modell nicht gar so schlecht, um so mehr, 
als über den Einfluß des als Senke auftretenden Kanals, der genügend kleine Dimensionen hat, 
um dort die Corioliskraft vernachlässigbar zu machen, schon eine ausführliche Untersuchung 
von Lauwerier vorliegt (TW 36). 


Eine über die Nordsee sich bewegende Depression wird durch eine Verteilung von Polen wie- 
dergegeben, die sich hauptsächlich bei der Bahn des Zentrums der Depression anhäufen. Daneben 
haben wir die durch Reflexion an den Küsten entstehenden Pole samt ihren Multipol,,schwänzen“, 
die hier in Form der Funktionen ® auftreten, in Betracht zu ziehen. Die größten und heute 
noch nicht übersehbaren Schwierigkeiten wird wohl die Laplace-Inversion der Faktoren ® 
bereiten. Wegen der Dämpfung werden wiederholte Reflexionen wohl kaum einen merklichen 
Beitrag ergeben, es sei denn, das Zentrum der Depression ziehe ganz nahe an den singulären 
Ecken vorbei, ein Tatbestand, der sich nur selten (aber doch gelegentlich mal) vortut. Die Re- 
flexionen an der schottisch-englischen Küste werden sicher wohl nicht vernachlässigbar sein, 
insbesondere wenn die Depression sich ziemlich nahe und parallel zu dieser Küste südwärts be- 
wegt, oder wenn sich nahe dieser Küste ein starker barometrischer Gradient befindet. Beide 
Erscheinungen traten beim Februarsturm 1953 stark hervor. Die Reflexionen an unserer Küste 
schließlich sind — für uns! — die wichtigsten, insbesondere wenn sie nicht zustande kommen, 


denn eben das bedeutet den Deichbruch. 


Obwohl also die analytische Behandlung von Problemen wie dem hier besprochenen ganz 
erheblichen Einschränkungen unterliegt und keinesfalls eine maschinell numerische vollständig 
ersetzen kann, so hoffe ich doch in diesem Vortrag gezeigt zu haben, daß die klassischen Me- 
thoden der angewandten Mathematik noch keineswegs als überholt zu betrachten sind, in enger 
Verbindung mit den numerischen einen nicht unerheblichen Beitrag zu unserem Verständnis 
bieten können und methodisch interessante und für die Weiterentwicklung der reinen Mathe- 
matik nicht unbedeutende Ergebnisse vorzuzeigen vermögen. 
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Ein Beitrag zur exzentrisch belasteten Kreismembran 
Von Max Hieke 


Wie im Falle der konzentrisch angeordneten, kreisförmig begrenzten, gleichförmigen Belastung der 
schwingenden Kreismembran wird auch bei der exzentrisch gelegenen yleichförmigen Belastung die Webersche 
Beziehung zur Integration der inhomogenen Schwingumgsgleichung herangezogen. Mit Hilfe des a 
positionsprinzips gelingt es, die Lösung der inhomogenen Schwingungsgleichung für die am Ra ein- 
gespannte Kreismembran herzustellen, wenn das belastete Kreisgebiet ganz dem Innern der Membran an- 
gehört. 


Weber’s relation, which has been used to integrate the inhomogeneous equation of motion of a circular 
membrane acted upon by a uniform load of circular boundary concentric with the membrane, is now being 
used to solve the corresponding problem for a non-concentric load that does not extent beyond the periphery 
of the clamped circular membrane. T’he solution of the problem is constructed with the help of the principle 
of superposition. 


Comme dans le cas de la charge uniforme de la membrane circulaire oscillante qui est &tablie concen- 
triquement et limitee circulairement, aussi dans le cas de la charge uniforme ewercde excentriquement la 
relation de Weber pour lintegration de lequation non-homogene des oscillations est appliquee. O’est ü 
V’aide du principe de superposition qu’on reussit ü &tablir la solution de l’&quation non-homogene des 
oscillations pour la membrane circulaire fixee au bord si la region du cercle charge appartient toute entiere 
a l’interieur de la membrane. 


Tarıke Kak U B CAIy4ae KOHUEHTPHYecKU PAcupejleJleHHON, OTPAHHYEHHOH OKPYFKHOCTEIW, 
PaBHOMePHOÄ Harpy3KM KosleÖlllomelich MeMÖPaHbI, MI HHTETPHPOBAHHA HEOAHOPOAHOTO 
mmuddepeHnmasnbHorO YPABHeHMA Kosle0aHnü U B CiIy4ae 3KCHEHTPHYECKH PACcHOJIOFKEeHHON 
PaBHOMepHoä HaTpy3KH IIPuUMeHAeTcH COOTHOMIeHHe Bedepa. Ilpıu moMmoINm NHpHHINMTMa HaJlo- 
JKeHUA YNaeTca HOAIy4YUTb PellieHme He0AHOPONHOTO YPABHEeHuN KoleÖaHnii IH MeMÖpaHhı, 
3akpelIlIeHHOÜ BAONb Kpafl, B IIPe]NIIOJIOSKEHUH, YTO BeChb KpyT, IIOABepraeMklä Harpyake, 
IPHUHAJLIeKHT BHYTpeHeii yacTu MeMÖPAHnI. 


Für die beabsichtigte Analyse der Schwingungen einer an ihrem Rande eingespannten 
Kreismembran vom Radius r=a sei der Mittelpunkt O der Mittelebene der Membran zum 
Ursprung eines ebenen Polarkoordinatensystems {0,r,9} gewählt. Dann genügt die vertikale 
Verschiebung w(r, o, it) eines beliebigen Punktes der Mittelebene unter transversaler Belastung 
der Membran der inhomogenen partiellen Differentialgleichung 


1 ow(r, 9, 


Aw(r, 9, 1) — 2 EB = Bpfr, 9), 
mit der Randbedingung 
lw(r, 9, !lr-a = 0. 
Die Konstante ist B—= — n . Dabei bedeutet c = £ die Schallausbreitungsgeschwindigkeit 
& 


innerhalb der Membran und 4 = «&- h deren Masse je Flächeneinheit. T ist darin die Spannung, 
& die Massendichte und h die Dicke der Membran. Die Funktion p(r, 9, t) ist die auf die Flächen- 
einheit der Membran bezogene Last. 

Es sei nun eine zur Zeit {= 0 plötzlich einsetzende Belastung der Membrangrundfläche 
angenommen, die der Verteilungsfunktion 


Yac Ipı(r,9) für >0, 
AR ER A Pa a 
entspricht. 

Als eine einfache Belastungsfunktion, mit der wegen des hier gültigen Superpositions- 
prinzips eine Reihe andrer Belastungsfälle erzielbar ist, werde die gleichförmige Belastung 
eines bestimmten Teilgebiets @ der Membranfläche gewählt. Das Teilgebiet G werde als ganz 
dem Inneren des Randkreises r = a der Membranfläche angehöriger, zu dessen Zentrum O 
exzentrisch gelegener Vollkreis vom Radius R=b angenommen (Bild 1). Es kann stets so 
eingerichtet werden, daß der Mittelpunkt M des kleineren Kreises auf der von O auslaufenden 
Halbgeraden 9 = 0 gelegen ist. Sein Abstand OM von O werde mit r, bezeichnet. Ein auf M 
als Ursprung bezogenes neues Koordinatensystem {M, R, y} geht durch Translation längs der 


Halbgeraden @ = 0, um die Strecke OM = r, aus dem vorigen hervor. Zwischen den Koordinaten 
beider Systeme bestehen die Relationen 


TC Ssp=Rcosy-+ T,, 


rsinpo=Rsiny, 
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eu nn 


woraus folgt 


in nn nn sine en m nn nm 
- ana jun 


iR = (r@ 89 — 1, + rin: p, 


Bezeichnet G die belastete Kreisfläche im Membraninneren und G das restliche Teilgebiet der 


Membranmittelfläche, so ist die Lastverteilung i 
an ung im Koordinatensystem {M, R,y} durch die 


in G, 


gegeben. In dieser Funktion bedeutet p, die kon- 
stante Flächenbelastung. Diese ist nur innerhalb 
des m Me7as auf O exzentrischen Vollkreises R= b 
von ’Null verschieden, und wegen der Radialsym- 
metrie vom Winkely nicht abhängig. 


pHR,Yy) = pi(R) 1% in G, 


Bild 1 Bid z 


Zur Integration der inhomogenen Schwingungsgleichung wird zunächst ein Integralsatz 
verwendet, der unter bestimmten Voraussetzungen über die Funktion F(£) des komplexen 
Arguments {=£-+ in für einen positiven reellen Parameter t folgendermaßen lautet 


v iv 
1 RAN © < iu (i y ij = 
Sri j e® FOL = 3% Ak (4—1)! u Ze 


c 0 für t<0. 


In dieser Formel bedeutet F(£) eine rationale oder transzendente Funktion des komplexen Argu- 
ments &=&+ in mit den m Polen Z, von der Ordnung I, Die AY sind die Koeffizienten 
der Laurentschen Entwicklung der Funktion F(£). Die Funktion selbst genüge der Forderung 
F(Ö) — 0 für {> . Das Integral ist als Hakenweg C in der komplexen Z-Ebene nach dem 
Bild 2 zu erstrecken, der alle PoleZ, zur Linken läßt und über einen aus dem Ursprung Z = 0 
geschlagenen Halbkreis mit unendlich großem Radius geschlossen wird. 

Die obige spezielle Belastungsfunktion 


a vr = |Ps für R< b, 
PKRY)=PKR=-PÄRM=-|g für R>b, 


ist nur innerhalb des zu 0 exzentrischen Vollkreises R = b von Null verschieden. Mit diesem 
von y völlig unabhängigen Störungsglied geht nach der Ausführung der Translation die ur- 
sprüngliche Differentialgleichung in die inhomogene, partielle, radialsymmetrische Differential- 


gleichung 
Fw*(R, t) 1 5w*(R, h) 1 Zw*R, ty 


SR? R 3R @ 9 
mit der Randbedingung 


B p*(R, I) 


[vH R, th. = 0 
und der Störungsfunktion 
[pi(R) für 1>0, 


PIEA=| 9 für 1<0, 


über. Zur Vorbereitung der Lösung dieses Problems geht man sachgemäß von der folgenden 
inhomogenen, radialsymmetrischen, partiellen Differentialgleichung aus 
en 


Su —i Bi er. 
atom ae BbnJa Me 
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Hierin ist die reelle Zahl s ein willkürlicher Parameter und J.(s R) die Besselsche Funktion 
nullter Ordnung des Arguments s R. Man findet als partikuläres Integral der inhomogenen, 
partiellen Differentialgleichung die Funktion 


u*(s, Rd) = — Sr her ei _ Je di 


Dabei ist zur Abkürzung k =£/c gesetzt worden. Der zeitunabhängige Faktor des Integrals 
u*(s, R, t), nämlich 


2— 
genügt der gewöhnlichen, inhomogenen Differentialgleichung 
DIPS AUL/ACER " e 


Bedeutet ferner J,(s b) die Besselsche Funktion erster Ordnung des festen Arguments sd, so 
ist die Funktion y* = J,(s b) J* ein partikuläres Integral einer gewöhnlichen, inhomogenen 
Differentialgleichung, die zwar den gleichen Differentialoperator wie die vorige, aber das Stö- 
rungsglied B b p, Jı(s b) J,(s R) besitzt. Nach dem Superpositionsprinzip ist dann das uneigent- 


liche Integral 
f Jı(s b) Jo(s R) 
* sLY en 
I =—2p0]| en ds, 


wenn es existiert, ein partikuläres Integral der inhomogenen, gewöhnlichen Differentialgleichung 
+RI*=Bbp, | 69) Is R)ds. 


0 


daI$ e .dI® 
dR®  RdR 
Wegen der Weberschen Beziehung 


f Pie meh, 
| 166) Js R) ds = |" a 
ö e 


ist aber die obige genau die Störungsfunktion, die der im Zeitpunkt ? = 0 aufkommenden sta- 
tischen Belastung mit der Verteilungsfunktion 


ix Ti, REF, 
er 


entspricht. Sofern das uneigentliche Integral existiert, liefert es demnach das partikuläre In- 
tegral der inhomogenen Differentialgleichung mit der vorgegebenen stückweise stetigen Störungs- 
funktion. Zur Berechnung des uneigentlichen Integrals 

| Jı(s d) Jos R) p“ 


ug für R>b 


N) 


geht man sachgemäß von dem in der komplexen z-Ebene (<=s- in) erstreckten Linien- 


integral 

P 1 

| ur ai dz für R<b 
aus. Darin bedeutet Hi(z b) die erste Hankelsche Funktion erster Ordnung. Diese Funktion 
ist in der komplexen z-Ebene unendlich vieldeutig und nähert sich in der positiv-imaginären 
z-Halbebene mit wachsendem z— oo asymptotisch dem Werte Null. Unter der Annahme 
R< b verschwindet die Funktion Hi(z b) J,(z R) in der oberen z-Halbebene mit nach dem 
Unendlichen gehenden z gleichfalls asymptotisch. Um die Hankelsche Funktion eindeutig zu 
erklären, werde die z-Ebene längs der negativen reellen s-Achse aufgeschnitten. Der Integrations- 
weg C ist dann so zu erstrecken, daß er beständig der positiv-imaginären z-Halbebene angehört. 
Dabei müssen die bei z= + k liegenden Pole des Integranden sowie der Pol erster Ordnung 
der Funktion Ai(z b) selbst durch kleine Halbkreise, die ganz der oberen -Halbebene angehören, 
umlaufen werden. Der Integrationsweg ist in Bild 3 dargestellt. 
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Da der Integrand in der positiv-imaginären z-Halbebene asymptotisch verschwindet, 
erhält man nach der Cauchyschen Formel 


Be (s b) J(s R) 
R Er weg ds =ni{ZRes). 
Insbesondere erhält man das zur Stelle z= 0 ’ge- 
hörende Residuum anhand der Definitionsgleichung 
der Hankelschen Funktion 


Hi(zb) = J,(zb) +iN,(zb) 


durch eine Untersuchung der darin enthaltenen 
Neumannschen Funktion N,(z b) in der Umgebung 
ihres Poles. 

Damit erhält man nach der Berechnung der 
übrigen Residuen die folgende Beziehung 


Ib) Is R) 
2| * 2 rer 


Hi(k b) — 


+ ED 7 Ai für R<b. 


er 


0 
Berücksichtigt man hierin nochmals die Definitionsgleichung der Hankelschen Funktion 
Hi(kb) = J(kb) +iN,(kb), 
so erhält man 


(A6b)I6@R), 1 Ny(kb) Ik R) „. 
j ae = — pa Ze Mür R<b, 
Für die Berechnung des uneigentlichen Integrals 
[ J(sb) Js R) 5, 
_R 


im ba b< R ist der ET RIDE Ausgang das geschlossene Kurvenintegral 


er b) — ai dz für b<R. 


Auf einem Halbkreis mit ih großem Radius in der positiv-imaginären z-Halbebene ver- 
schwindet die Funktion Hi(z R) J,(zb) unter der Annahme b < R asymptotisch. 

Die erste Hankelsche Funktion nullter Ordnung Hi(z R) ist bei z= 0 logarithmisch 
verzweigt. In der längs der negativen reellen Halbachse aufgeschnittenen z-Ebene ist die Funk- 
tion eindeutig, und der Integrationsweg kann in den Nullpunkt hineingeführt werden. 

Mit dieser Abänderung an der Stelle z = 0 verläuft der Integrationsweg nach Bild 3, wie 
im vorhergehenden Falle beim Auftreten der Hankelschen Funktion erster Art und Ordnung. 
Die Cauchysche Formel ergibt hierfür 


= 


da —-nifZRes} für bD<R, 
oder B 
fe R) Jı(s b) FR DD REN S en Jkb) für DR. 
Aus den Relationen 
Jı sb) = — Jı(s b),, 


Hi(-s b) = His R)— 2 Ju(s R) 
und der Definitionsgleichung der Hankelschen Funktion H}(k R) folgt dann 


Je EACDILFCE Erg 
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oder 
Er o(s R) Jı(sb) ._ er Jı(k b) No(k R) = 
j ar ee Ben Tree tr DH, 
) 


Auf dem Kreis vom Radius R = b gehen die in den zwei verschiedenen Gebieten erklärten 
Funktionen ; 
NEE ud 


| A — Ik N für R<b D 
NR) AR), pe Bo, 
2k z 
stetig ineinander über, da für R = b die Gleichung besteht 
2 


— p  Nilk b) Jo(k b) = — No(k b) Jı(k b) . 
Setzt man hierin vorübergehend kb = o, so erkennt man in dieser Gleichung die Funktional- 
gleichung der Zylinderfunktionen des Arguments o für den Spezialfall v = 1, nämlich die Relation 
2 
Jo Na — Jo N = 70 = 
Um der Randbedingung 
lw(r, 9, Ylr-a = 0 

zu genügen, braucht man ein partikuläres Integral der homogenen Differentialgleichung für 
w(r, 9,1) = J(r,p) e''®. Man findet ein solches als Integral der zuvor homogen gemachten 
partiellen Differentialgleichung für diese Funktion. Sie lautet, wenn die Zeitvariable zuvor 
separiert wird 

2J(,9) 100,9) 1 8%J(r,Y) 

or? Ta or r 0 


+RJ,g9)=0, 


2 
wobei k? — 5 gesetzt wird. 


Dieser homogenen Gleichung genügt, von einem konstanten Faktor abgesehen, formal das 
in Gestalt einer unendlichen Reihe mit beliebigem Koeffizienten A, angesetzte Integral 


IKn)= E A,Jlknew. 


t=—o 


Unter der Annahme, daß sich die Konstanten A, so bestimmen lassen, daß die Reihensumme 
existiert, erhalten wir das gewünschte Integral der homogenen Differentialgleichung. Mit dieser 
Reihe und dem partikulären Integral der inhomogenen Gleichung machen wir den Ansatz 


Bbp,Jı(sb u ur 
y(s 7,9) = — en (I R)+ 3 AnJılkn) er). 
Aus der Randbedingung 
[y6s, T, P)]r=a mas 0 » 
folgt 


JG R)=— & AJılka)err, 
mit | 
R,=yR +@—2ar,cosp. 
Die gesuchten Koeffizienten A, erhält man nun aus der Koeffizientendarstellung 
A,n=— BR. | 36 R0 etrrap. 
2nJ,(ka) ve 
0 


Sofern sie Konvergiert, liefert dann die Reihe 


_ Bbpo SS Julkr) | Js DIRRE R) er dp 


27 „AT lka) u — u 


S 


0 
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ein partikuläres Integral der homogenen Differentialgleichung für die früher definierte Funktion 
I(s, 1,9) = I*(s, R), welche die Ausbiegung der eingespannten Membran auf dem Rande zum 
Verschwinden bringt. Da stets b< R, ist, folgt wegen der Vertauschbarkeit der Integrations- 


reihenfolge für jedes einzelne Reihenglied eine Gleichung 


© 2n 2n 
[ Jsb)S Jos Ro) ei"? dp Jı(k b) [ No(k Rı) ei"? dd 
ieeund nee ud oil en nbenine gene, 0 — seen 
} 2 — KR? Ei an" Ik 


wobei noch wegen der Freiheit in der Bezeichnung der Integrationsvariabeln unter dem zweiten 
Integral für die Winkelveränderliche 9 der Buchstabe # gewählt ist. 
Damit sind die zeitunabhängigen Integrale des vorgelegten Problems gefunden. Für diese 


Integralfunktionen gelten dann die folgenden Definitionen 
27 


1 1 Nı(kb) J(kR) mb =, Ju(k r) eiv® | 
EI (IN near 0 rl, ONE HD —ivd 
Bp, mt" 2k red 2 Tka)2n has 
o für ADS, 
und 

2 
ei _,Jı(kb)NokR) mb = Ju(kr) eiv» En 
Br aaa N Nee N,(k R,) ei? dd 

. fur DENE 


Die Eindeutigkeit dieser Funktionen in der {-Ebene erkennt man mit Hilfe der Additionstheoreme 
für Zylinderfunktionen. Diese lauten 


a) = 2 Julkır) Jullers) e®0®r für aller, 
v.—o© 
NER) = AI UHREN Nylkr,) e@r, fürar rn: 


oo 


Ni(kR)= 3 Jukr)Nkr)eer, für n<r, 
vw—=— om 


und für die Werte r auf dem Rande r = a, wenn man zugleich die Veränderliche p durch die 
Integrationsveränderliche # ersetzt, insbesondere für die Neumannsche Funktion unter dem 
Integralzeichen 
NK.) = 2, der) N.(ka) €? fire re 
cv _—o© 
Mit dem letzteren Theorem erhält man nach Ausführung der Quadraturen für die Integrale die 
Reihendarstellungen 


(kb) JkR) md J,(k b) >23 Ju(kr) Ik ro) Nokk a) , 


1 1 N 3 
7 =urte ” ur Di, 


2k 2k „. J,(k a) 
und 
1 Jı(kb)NykR) mb  J(kr) Julkr) Nylka) . F 
ie . en itich ia ER A re 
Hr 1,(k) = mb 3% 5 1928 T.ko eivo für D<R. 


Ersetzt man nun in diesen Reihenausdrücken die entsprechenden Funktionen durch ihre Ad- 
ditionstheoreme und nimmt für die Neumannschen Funktionen deren Potenzreihenentwick- 
lungen hinzu, indem man für die logarithmischen Glieder derselben wegen der Gleichung k = £/e 
die folgende Schreibweise benutzt, 


Ping yr 
In (%) — Inö+InS,: 


2C 


wobei die Ausdrücke In er ‚In r, In b, usw. reelle Zahlen sind, so bleibt als unendlich vieldeutige 
c 


Funktion jeweils die Funktion In{ übrig. Durch gesonderte Betrachtung der einzelnen Koeffi- 
zienten der Potenzen von e'?, bemerkt man, daß deren mehrdeutige Anteile sich in der unten- 


stehenden Anordnung zusammenfassen lassen 


In) ine + ind, USW. 


Ik b) Inkk 1) Inlk ro) Ink @)) 1; 
In(k a) 


’ 


oe 1 I1(k b) Inkk FI In(k To) — 


14 
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und ersieht daraus, daß deren jeder für sich identisch verschwindet, und dadurch die logarith- 
mische Verzweigung aufgehoben wird. Die Funktion /,(£) ist danach in der $-Ebene eindeutig. 
Auf gleiche Weise zeigt man die Eindeutigkeit der Funktion I;(2). Beide Funktionen sind dem- 
nach im Integrationsgebiet eindeutig. . 
Jetzt soll gezeigt werden, daß beide Funktionen im Unendlichen der {-Ebene verschwinden. 
Dazu geht man sachgemäß wieder von der Reihenentwicklung der Integralfunktion aus, indem 
man die Koeffizienten gleicher Potenzen von e'? zusammenfaßt. Das ergibt für den Koeffi- 


zienten der Potenz e'*r der Reihe für EA 1,(£) den folgenden allgemeinen Ausdruck 


B po 
NIDEATDPAUTRERLLDEIT 2 De ro) Nu(k a) | 


ab 


Dabei ist der allen Reihengliedern gemeinsame Faktor By weggelassen worden. Man kann 


dafür schreiben 


Ju(kr) Jn(k ro) Ni(k b) I _N,(k b) J„(k a) 


Führt man im Faktor innerhalb der Klammer für die Zylinderfunktionen ihre asymptotischen 
Werte 


Jı(k b) Nn(k a) 


10) = |, ein + =) 2]. 


Ba l\r 
x = zn en +3) 5]. 
ein, so läßt sich zeigen, daß, wenn darin das Argumento den Erfordernissen der Rechnung 


entsprechend gleich ka bzw. kb gesetzt wird, der eingeklammerte Faktor asymptotisch durch 
den Ausdruck 


hir alsnsve—(n+ 2) 3] 


angenähert werden kann. Für die Tangensfunktionen gelten die folgenden Periodizitätseigen- 
schaften 


3 1 
alas La] = cotg er i 
und 
N k 
1 —tglka — — | für ungeradesn, 
m 4 
cotg |k a — Ink Sl = 
2/2 ı 
cotg ka EZ | für geradesn, 


und für ihre Fortsetzungen in die komplexe &-Ebene, mtl =&+ in, Additionstheoreme in 
der Form 


r b nr 
sin 9E —16©in?2 


nb 
Sb c c 


cotg — = 5 iR; x 
ut ) 
Co] 2 NT _ 082 & 
C C 
a Be 
ein 2 
sa C c 
tg — a Ser san - age $) 
[6 € 


cos? 


ae für große positive Werte von n in der positiv-imaginären £-Ebene die asymptotischen 
erte 


5b i Ca 
Yoga le er 1 
a a 0 tg Bann 
folgen. Entsprechend erhält man für große negative Werte von n die asymptotischen Werte 


b 
org Kr, NE 


ir 
S 


nn BEL =: Tz Fra 
a u 3 Ee 
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Mit Verwendung dieser Werte zeigt man das asymptotische Verschwinden der Faktoren in der 
Klammer, und damit jedes beliebigen Koeffizienten der nach ganzzahligen Potenzen von ei» 
fortschreitenden Reihe für das Integral I/,(£) im Unendlichen der komplexen £-Ebene. 

Bei der Untersuchung des asymptotischen Verhaltens des zweiten Integralanteils I,(£) 
ergibt sich insofern eine Komplikation, als für die in ihm enthaltene Funktion N,(k r), wie aus 
der Aufstellung hervorgeht, zwei verschiedene Additionstheoreme zu verwenden sind, je nachdem 
nämlich r<r, ist. 

Die Zusammenfassung der Koeffizienten einer Potenz von e'? in diesem Integralanteil 
ergibt dann, wenn der allen Gliedern gemeinsame Faktor weggelassen wird, gemäß den zwei 
Reihenentwicklungen die beiden Darstellungen 


[7.0« b) Ile r) Nalkry — ED J N To) Nn(k a 
— (kb) Inkkr) Na(k ro) 11 are n für r<r, 
und 
[Walk 8) Indk ro) Nacken — EN I u ro) Nu(k a)| 
= Ik Ile N | fü n<r. 


Mit Hilfe der vorhin angeschriebenen Additionstheoreme für die Tangensfunktion, zu denen 
noch ergänzend das folgende 


kr = für ungeraden, 


— cotg 7 


uln-+4)3]- 
4 


hinzuzutreten hat, zeigt man auch das asymptotische Verschwinden der Koeffizienten des 
Integralanteils /,(£) in der komplexen {-Ebene und damit das asymptotische Verschwinden 
beider Integralanteile überhaupt. 

Damit ist die Möglichkeit gegeben, durch Anwendung des eingangs angegebenen Integral- 
satzes, in einer dem gestellten Problem angemessenen Form, die zeitabhängigen Faktoren der 
Integrale zu bestimmen. 


Wir wählen den Integralsatz in der folgenden Form 
Er 
w(r, 9,0) = i er v=1)2. 


2ri e 
1) 


tg Ik LE A für geradesn, 


Der Integrationsweg C ist der in Bild 1 dargestellte. Die Voraussetzungen für die Anwendung 
des obigen Integralsatzes zur Berechnung der Durchbiegung w(r, 9, f) sind erfüllt, da die Funk- 
tionen /,(£), die mit den Funktionen I,(k) über den Zusammenhang k = £/c sofort gegeben sind, 
diesen genügen. Wir beginnen mit der Bestimmung des Residuums an der Stelle © = 0 als des 
umständlicheren Teils der Aufgabe. Sie wird am einfachsten erledigt, indem beide Funktionen 
I,(£) e'* in der Umgebung £ = 0 in eine Laurentsche Reihe entwickelt werden. Wir schreiben 
zur bequemen Übersicht die zeitfreien Integrale noch einmal in der für unser Vorhaben zweck- 
mäßigsten Form an 


2n 


} —, J,(kr) eiv» 
5% I,(k) = z LIE NE BEN PN A | Ny(k R,) e'?? dd 
0 vw=-——o 


2k 212% Jılka) , 
0 fun RD), 
2n 2 
1 Jı(kb)Ny(kR) mb | Jı(k r) e'’® J ee 
bl 22 u 10 ET (Iah —— —— [| Ns(k R dd 
na zu I Te | NkRde 
0 AUrN DER: 
1 * ee 
Eingangs beschäftigen wir uns mit dem zur Funktion Er I,(£) gehörenden Durchbiegungs- 


0 
anteil. Diese eindeutige Funktion hat nur einfache Pole, sowohl an den Nullstellen der Bessel- 
schen Funktionen aller Ordnungen, als auch an der singulären Stelle & = 0 der Neu mannschen 
Funktion erster Ordnung. Wir erledigen die Bestimmung des Residuums an der Stelle ln 
14* 
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zwei Schritten. Zuerst betrachten wird den der inhomogenen Differentialgleichung entsprechenden 
! 1 2 F 

Integralanteil, also die ersten zwei Summanden der Funktion Bn I,(£). Von diesen Gliedern 
0 ” 

setzen wir nur den Anteil der Laurentschen Entwicklung an, der tatsächlich zum Residuum 


- 1 - 
etwas beiträgt. Unter Berücksichtigung der Tatsache, daß die Funktion Bp I,(£) insgesamt 
0 


keine logarithmische Verzweigung mehr hat, lautet ihr dem Integral der inhomogenen Differen- 
tialgleichung entsprechender Anteil folgendermaßen 

1 [b® Ir na . H 

— 1 - —_ ... f RR: 

£ zind+z(R b2)t + E ür < 
Dabei sollen die Punkte solche Glieder der Laurentreihe andeuten, die zum Residuum nichts 
beitragen. Die Bestimmung des Residuums des dritten Summanden, der dem Integral der homo- 
genen Differentialgleichung entspricht, lassen wir noch beiseite, da er auch im Funktionszweig 


Er 1;(£) und zwar als zweiter Addend auftritt. 
0 


Wir lösen zuvor die einfachere Aufgabe, das Residuum für den Summanden dieses im Ge- 
biet b< R definierten Integrals zu finden, der dem Integral der inhomogenen Differentialglei- 
chung entspricht. Man erhält dafür, da auch hier keine logarithmische Singularität für die ge- 
samte Funktion existiert, die Laurentreihe 

1 5 ” 

Fi Herrn für u: <ichr 
Dabei ist nur das Glied derselben ausgeschrieben, welches einen tatsächlichen Beitrag zum Re- 
siduum liefert, die übrigen, die das nicht tun, sind durch Punkte angedeutet. Man bemerkt, daß 
die bisher erzeugte Funktion 


2 
ind + R— für Rep: 
und die Funktion 
2 
in für WR 


auf dem Kreise R = b stetig sind. 


Nun erübrigt noch die Berechnung des Residuums für das beiden Teilfunktionen zugleich 
angehörende partikuläre Integral der homogenen Differentialgleichung an der Stelle = 0. 
Das partikuläre Integral der homogenen Differentialgleichung ist 
2n 


| No(k R,) e-'r® dB. 
0 


Um die in der vorausgehenden Argumentation öfter verwendete Tatsache zu begründen, daß die 
logarithmische Verzweigung in der Integralfunktion selbst nicht mehr auftritt, soll zunächst 
gezeigt werden, daß das logarithmische Glied der Reihendarstellung separiert werden kann. Mit 
Hilfe der Potenzreihe für die Neumannsche Funktion nullter Ordnung 


ykR, 
2 


_abJ(kb) DS Jılkr) 
2n 2k „= _Jı(ka) 


2 | 
Nolk Ro) = = Ik Röin E84 2 (HR) — ZIRR) +4 IKER)—::), 


und mit Hilfe des Additionstheorems für die Besselsche Funktion läßt sich dieses partikuläre 
Integral nämlich folgendermaßen darstellen 


kb) &, Iulkr) Julkr) Ilka) u...yk 
af 4 v v {1} v i 
Et J,(k a) 7207 


va— 


2n 


er | {Jo(k R,)In R, + 2 (Jg(k Ro) —  - -)} erir? dd. 
0 
Bezeichnen wir in dieser Entwicklung die zweite Reihe der Kürze wegen mit $*, so lautet sie 


mit nochmaliger Verwendung des Additionstheorems für die Besselsche Funktion nullter 
Ordnung 


b J(kb) Ss Jılkr) 
2n k Jy(k a) 


v=—o 


EITALDELLRONEL IE" 


k 
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Davon hebt sich aber das erste Glied, also die logarithmische Singularität gegen die logarithmi- 
schen Glieder der in den beiden Teilgliedern der Integralfunktion, nämlich die, welche 
Ni(kb)J(kR) . 1 P | 
ab yk sn Bp, Türe =D, 
enthalten und diejenigen, welche 


Jıkb)N.(kR) . 1 N 
IT b 5) k ın Bp I;(k) für , = R 
enthalten, heraus. Man bemerkt das, wenn mian in den letzgenannten Funktionen, die darin auf- 
tretenden Zylinderfunktionen J,(k R) und N,(k R) gleichfalls durch ihre Additionstheoreme aus- 
drückt. Die logarithmische Singularität ist also gemäß der Behauptung im Integral nicht mehr 
vertreten. 


Aus diesem Grunde wurden die vieldeutigen logarithmischen Teilfunktionen in den vor- 
angehenden Rechnungen, die für die zwei verschiedenen Definitionsgebiete gesondert durchge- 
führt wurden, garnicht erst angeschrieben. : 

Für die Berechnung des Residuums an der Stelle & = 0 ist es deshalb nur notwendig, vom 
Integral, das der inhomogenen Differentialgleichung zukommt, den Anteil 


2n 


Zr -2r0) 5 ns ’ a | [ol Ron R, + 2(Ik R)— Z Ik R) + -- \ewao 


vw —o 


0 


zu betrachten. Für die Ausrechnung der Integrale ist festzuhalten, daß bei ganzzahligen Werten 
des Index v, die hier lediglich auftreten, für die Besselschen Funktionen die lineare Abhängig- 
keit besteht 


Jo(kr) = — 1)" Ju(kr), für ganzzahlige Werte n. 
Unter Beachtung der Gleichung k = £/c folgt dann für das Residuum 


2n 
b2 & r\e . 
Bst > = ee [mr re dl: 
52 ) An, —, a / 


denn die unter dem Integralzeichen in der runden Klammer stehende unendliche Reihe aus 
Besselfunktionen gerader Ordnung n >, liefert dazu keinen Beitrag. 


Dieser aus der Funktion 3/* gewonnene Ausdruck stellt bis auf den Faktor b?/2 eine Reihen- 
entwicklung dar, die sich in reeller Schreibweise in die Gestalt bringen läßt 


[6,0] 


oo [ol en v 
1 DS: (2) | I) ei? dd — er sn Ex (2) (&, c0Svp + f,sinvp), 
0 


an v= —o v=1 


mit den Koeffizienten 


2n an 
= al 9) cosn dd, Br -- | no) sinnddd, 
7 
0 0 


und dem Randwert 

f#) = In R,=Inya Hr —2ar,cos®. 
Diese Reihe konvergiert für alle Punkte r< a innerhalb des Kreises r = a gleichmäßig und 
stellt dort eine harmonische Funktion des Arguments z = re‘? dar. 


Andererseits aber stellt die linke Seite der Gleichung, die im Innern des Kreises r< a 
gültige Fourierentwicklung des von O nach 2 r längs des Kreises r = a erstreckten über den 
Randwert 


f#) =1InR, 
genommenen Poissonschen Integrals 
2n 
ade I) dd erh = (2) dr 
DER ea 7) @eosvp +5 0) 


0 
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dar. Insbesondere hat das Poissonsche Integral über den Randwert f(9) = In R, genommen 
und längs des Kreises r = a erstreckt, den speziellen Wert 


w.(r, 9) = In Ya + AR — 2a@rr,cosp—Ina. 


Demnach lautet der mit der eingangs definierten Belastungsfunktion p(r, p) zur Zeit t = 0 auf- 
kommende Beitrag des Integrals, der von der Stelle £ = 0 herrührt, 


> b? 2 

ind +4 (RB) — zul, für R<b, 
w(r,9)=Bp B2 ja 

„In R—- 5 wir, 9) für b<R. 


Dabei ist w,(r, 9) der vorhin ermittelte Wert des Poissonschen Integrals. 


Nun hat aber das Integral der homogenen Differentialgleichung, das aus Gründen be- 
quemerer Rechnung in entwickelter Form 


mb 705) Sr Fkn) Jdkr) Nikko) a, 
ar) 3,7 ke | 
benutzt werden soll, außer der Singularität bei £ = 0 noch weitere unendlich viele einfache Pole 


an den einfachen Nullstellen der Besselschen Funktionen n-ter Ordnung Jn(k a). Diese unend- 
‚ich vielen einfachen Pole des Integranden liefern gleichfalls Beiträge zum Linienintegral 


ee. ee 
Imi De ca) | 


eve dk. 
Er A 


[0] 
L Car: 
Die von Null verschiedenen einfachen Nullstellen der Besselschen Funktionen J, ) sind 


reell, und wir schreiben für die einzelne Bessel-Funktion ganzzahliger Ordnung n vom Argu- 
ment o die einfachen Nullstellen in der Weise 

0m 9; a EB. 
Diesen Zahlen entsprechen die von Null verschiedenen Z-Werte 


ee, RER hen 
Die Nullstellen liegen auf der reellen s-Achse der komplexen {-Ebene und sind bei der Ausrech- 
nung des Linienintegrals längs des Hakenweges nach Bild 2 auf kleinen Kreisen zu umlaufen. 
Man erhält auf diese Weise den Beitrag der einfachen Pole zum Linienintegral 


(R bon r\ Cora) Cu. a 
Bi; —J. N} N 
b 6) 2a ı( 2 ) .[ ii Ai zi ; 
e 7 er 9 “= —&0 2 2: # Eu s i ern? eoxt ”. 
Sax ge v Zr 


Da der Beitrag der bei & = 0 gelegenen Nullstellen n-ter Ordnung der Besselschen Funktion 
n-ter Ordnung bereits ermittelt wurde, ist bei der Summation über x der Wert x = 0 ausgelassen. 
Diese Vorschrift ist durch den Strich am betreffenden Summenzeichen ausgedrückt. Beachtet 
man die Beziehungen 


ferner 
La e) 2) 
N, | BRIAESN, (0, ya la A 
\ c u (On) T Onk In(Onx) 


welch letztere aus der Wronskischen Determinante folgt, so läßt sich dieses Resultat zusammen- 


fassen zu 
Hi Be, >) J, Ei | 2‘ 


s-ab > N VER . E 


=—O x=—o 


eivp Eon! 
8 B 


% ” ie, A| 
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Wir schreiben die Reihensumme reell, nachdem wir sie zuvor mit der Konstanten B P, multipli- 
ziert haben. Das Produkt B’p, S sei mit W;(r, 9, t) bezeichnet, womit folgt 


Cu b\ in (C, 6 
wi, De g. c EA )2.( ) cos & 


ct 


nn 0x [Jo@,)P? a 
ee 
+2 r cos dv cos * ct 2 
Dr 2 019,1 ale 

Dabei sind, wie herkömmlich, die Nullstellen der Besselschen Funktion nullter Ordnung mit 

00x == 0x 
und entsprechend die Frequenzen mit 

as = &, 


bezeichnet werden. 

Mit Verwendung der obigen Bezeichnung Wj(T, 9, i) für die freien Schwingungen läßt sich 
die Lösung des Schwingungsproblems der eingespannten Membran mit plötzlich einsetzender 
gleichförmiger Belastung, die auf einen zum Randkreis exzentrischen Kreis verteilt ist, in der 
Form schreiben 


Ei 9-1) für t<O 
Ilm) für 1>0, mit v=12, 
wobei sind 
b2 1 b2 E 
w(,90) = Bps In b+ a er walr, 0) +Wn 9b), für R<b, 
b2 b2 
w(,9,D)=Bp > In R— 7 wo(T; o)) + Wil, 9 D) , für BER. 


Daraus ergibt sich wegen w(r, o, f) = 0 für die Zeiten t< 0, daß die Funktion [W,(r, 9, Hl-o = 
W;(r, 9, 0) gleich der Reihenentwicklung der stetigen Funktion 


b 3 
Lanfennad mung) nr Reh 
w(r,g) = p2 p2 
— Bp,( m R— 5 (6, 9)) für b zul; 


zur Zeit t = 0 ist. 

Mit dieser Lösung hat man durch die Anwendung des Überlagerungsprinzips auf diese, als 
der Grundlösung, das Mittel in der Hand, um für eine weitere Anzahl komplizierterer Belastungs- 
funktionen die Lösungen der betreffenden Differentialgleichungen zu berechnen, oder um dort, 
wo eine exakte Bestimmung auf diesem Wege nicht durchführbar ist, mindestens eine Näherungs- 
lösung des betreffenden Problems herzustellen. Die Möglichkeit der Approximation ist für solche 
Probleme geraten, deren Lastverteilungen nicht axialsymmetrisch um den Mittelpunkt M eines 
zum Ursprung des Koordinatensystems {0, r, 9} exzentrisch gelegenen Kreises angeordnet ist. 

Damit die gegebene Lösung als Ausgangslösung für diese Fälle brauchbar wird, ist sie zuvor 
zu modifizieren. Die Änderung der Grundlösung besteht nun darin, daß man in ihr an Stelle der 
Winkelveränderlichen 9 das Argument  — 9,; einzuführen hat. Die darin eingeführten Größen 
Po, bezeichnen die festen Winkelkoordinaten der Mittelpunkte M, der einzelnen Lastverteilungen, 
aus denen die Gesamtbelastung aufgebaut zu denken ist. 

Ergänzend sei noch erwähnt, daß man für r,— 0 aus der obigen die Lösung der Differen- 
tialgleichung erhält, welche zu einer zum Ursprung O0 konzentrischen axialsymmetrischen Stö- 
rungsfunktion gehört. 

Für eine vorgegebene rhythmische Belastung aber, mit dem Störungsglied 


9! für t<0, 
DE Pır,g)sinot für 1>0, 


ist der Integralsatz sachgemäß etwa in der Form 


np) =— 5, DER den D= 1,2; 
zu wählen. Mit dieser Wahl ist aber schon von vornherein über die Anfangsbedingungen verfügt. 
Der von den Lagekoordinaten abhängige Faktor der Störungsfunktion p,(r, p) ist der gleiche wie 
vorhin und ® die neu hinzukommende Störungsfrequenz. 
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Die Randbedingungen sind gegenüber den vorigen unverändert beibehalten. Weiter ver- 
wenden wir in den folgenden Formeln die Bezeichnungen 


@ 1 4 
93, =, Ks BF Eu 


Man erhält dann mit diesem Ansatz nach längerer Rechnung, die bis auf ‚einige Veränderungen 
der für die statische Belastung durchgeführten analog verläuft, die Ausbiegung 


Bpo ir BA! Zen + War, n) sin ot + Wı(r, 9, ) für R<b, 
w(r, ®, )) = ; 


Bp Ihr b a 3 + Wulr, o) sinot-+ W(r, 9, p) für b<R. 


Die darin auftretenden Funktionen 


kb) Si Julk Ze 
we DT kr Nike, 


oder reell geschrieben 


Wir, g) Be As b) [doc 1) Ny(k a) Jolk r) + 2 5; J(k r) Jo(k 70) Ny(k a) cosD e) - 


k J,(ka) J,(k a) 
bedeuten einen Anteil der im Rhythmus der Frequenz w erzwungenen Schwingung, und 
x Jı(k b) = Jolkı r) we 
Weapdnb 170 ro) No(k a) 2 BE 


oo oo T (ko r) P) #. | 
> Tre UFER ut; 
+2 ps Jı(k ro) Ny(k a) cos v ap (2, — RK) Js(o, sin Zu, 


2 


die freie Schwingung der Membran. Die Konstanten k,, = = sind die den einfachen Nullstellen 
der Besselschen Funktionen J,(e,,) zugeordneten Wellenzahlen und k, = 2 die den Null- 


stellen der Besselschen Funktion nullter Ordnung entsprechenden Wellenzahlen. 
Manuskripteingang: 6.6. 1958 


Anschrift: Dozent Dr. Max Hieke, Institut für Theoretische Physik der Universität, Halle (Saale), 
Friedemann-Bach-Platz 6 


Der Spannungszustand einer Reihe von Kräften, 
die in der unendlichen Halbebene angreifen 
Herrn Prof. Dr. Ludwig Föppl zum 70. Geburtstag in Dankbarkeit und Verehrung gewidmet 
Von Walter Engl 


Der ebene Spannungszustand für die Belastung einer unendlichen Halbebene mit einer im Innern an- 
greifenden unendlichen Reihe von Einzelkräften wird behandelt. Die Lösung wird als hyperkomplexe Funk- 
tion angegeben, deren Teile die Spannungskomponenten darstellen. Mittels Residuumrechnung wird die 
Eindeutigkeit der Verschiebungen bewiesen. Mit einem überlagerten einachsigen Spannungszustand kann 
die Lösung auch als erste Näherung für einen Streifen (Lasche unter Zug) dienen. 


The two-dimensional state of stress in an infinite semi-plane with an infinite series of concentrated forces 
acting in its interior is being discussed. T’'he final solution is said to be a hypercomplex function, the parts 
of which represent stress components. By a residuum calculation it is shown that the displacements are unique. 


With a superimposed, one-dimensional state of stress, the solution can also be taken as the first approximation 
for a strip under tensile force. 


Il est question de l’ötat de tension d’un plan a deux dimensions pour la charge d’une suite infini de forces 
independantes agissant a linterieur d’un demi-plan infini. La solution est donnde par une fonction hy- 
percomplexe dont les parties forment les composantes de la tension. A l’aide d’un calcul de residu on de- 
montre la netteld du deplacement. 


Avec un etat de tension uniaxe superpose la solution peut de möme servir comme approximation pour une 
bande. 


PaccmatpuBaeTcH InIocKoe COCTOAHNE HANPFAHKEHHA IA HATPy3ku ÖeCKOHEeYHOH MOJIYIRIOC- 
KOCTU, HONBEepTaeMmoi N3 HYTpu ÖecKoHeYHOMY PAAy ONHHOUHEIX cu. PemmeHme naetca B Bune 
TUTIEPKOMIMIEKCHON PYHKUUN, YACTH KOTOPOH PaccMaTpuBawTcH Kak KOMNOHEHTEI HAIPsTKe- 
una. Ilpm moMoLM BbIUMCIEHUA BbIYUeTOB MOKA3LIBAeTCH EHUHCTBEHHOCTB CMEeIEHHA. C HAIIo- 


}XCHHBIM OJ[HOOCHBIM COCTOAIHHEM HAUPSBKCHHUA, PEINEHNE MOHMET CHYIKHTb KAK IIEPROE IIPHÖNN- 
FKCHHE JWIA IOJIOCHI (HAKJIAKU, HONBePpraeMmoH Tare) 
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Im Innern der unendlichen Halbebene greift im Abstand a von deren Rand eine unend- 
liche Reihe von Einzelkräften gleicher Größe P in Richtung der negativen y-Achse an. Die Kräfte 
haben untereinander alle den Abstand b. 


Bild 1 


Der von dieser Belastung hervorgerufene elastische Spannungszustand läßt sich durch 
folgenden hyperkomplexen Spannungstensor darstellen: 


OEL a es Dez 
= 1 Om+Di+ mn Pete | @ — an] + I6Em Dim Yrlctg|7 e+ an) 
n a/b 


+[@ m—2) + (6m +2) 72] = 
Teta) 


+ 4A(m + Diocg|” (z+aj) N 


sine] 7 (z-+aj) 


sin? 


Wir wollen zuerst zeigen, daß diese Lösung mehrere schon bekannte Spannungszustände 
als Grenzfälle enthält. 


a—o— 
ZUR; 
(o) = er a Er [5 (1). 


Unendliche Halbscheibe mit einer unendlichen Reihe von Einzelkräften an deren Rand an- 
greifend (vgl. Münchner Dissertation: W. Engl, S. 38) 


a—o— 
bo 
4.08] a | 
E20 we . 42): 
a akt. eier, 2) 
Durch Einzelkraft senkrecht belastete Halbebene (vgl. z.B. L. Föppl, D. &Z. III, S. 26) 
b> © 
= lem + my + Em Ni-m—dA, a; 
Arm z—aj z+a]j 
2m 2) 4 (6m 1 age 4 Er LE (3) 


Unendliche Halbscheibe mit Einzellast im Innern (vgl. E. Melan: Zamm 12 (1932); Sobrero: 
Theorie der ebenen Elastizität, Hamburgische Math. Einzelschriften (1934), siehe auch L. Föppl: 
D. &Z. 111, S.;101). 


Die Randspannung (o,),-o erhält man leicht durch Zerlegung der Funktion in ihre reellen 
Teile: 
2» a 
n-m T“ 
r=@+n. 


(0z),—0 = ! (m } 1) | 


Nachdem wir die einzelnen Grenzfälle der Lösung aufgezeigt haben, wenden wir uns ihr 
selbst zu und betrachten zuerst die einzelnen hyperkomplexen Funktionen, aus denen sie aufge- 
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baut ist. Durch elementare Rechnung erhält man deren reelle Teile: 
1 (Sin2 1— cos 2£ C0j 2 $ 
-| ie ni _— a4... 


jtgt= N 
1 (Sin2n 1— cos2£Cj2n 
a a 


1 sin2F6©in2n 
er ee Mi . 


1 [&in2n ee ah ek : 
En EM 


% 


5 Fe 1— cos25C0]2n 
PR Le le a 5 yarcaaz 
1sin2& nsin2E-Ein2n]l, 
+Pl2 ie IE N 
1 11—cos2£b0j2n n Sin2n(l —cos2560j2n +sin?25) _ 
im 2 mM u N? 
EaTZ N ua 22 Bee 
4 N® 
“ 1sin22-Gin?2n 7 sin?2&(1—cos2£60j 27 — Sin?2n) j 
ee 
KR n Sin2n(l—c0s25C0j2n +sim22) 
sin?l 4 N® 
pl? Sin2n( — 0082£00j27 + sin 22 Alam a 
4 N3 2 N? 
+ | 1 sin2E-©in2n N 
a N: . EEE ir 


Daraus berechnen sich die Spannungen zu: 


2 er 
® | Sin Wo) Sin" W+® 


Tim N 
S 1— 005 77.00 ya) 
= m>41) ay; y—a) = Nr 
1 eos. Sof - — W u om 
— [m —)y+(m+1ba] 2 N 
zuä 1— cos . ie b ey +4a+ an“ — En- u +4) 
ir y(m +1) NG. | 
mit der Abkürzung N#) — sin? 7 ag Sin: 7 ‚eo +4); 
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1 — cos _ -Eoj Fr (y — a) 
+(m+ D5, u od IN wre 
1 — cos TE Go I (y+ a 
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T sin Aug [m — 1, 1 m sr Ins] 
a Abm b | 3 NS” 771 
. an ; 2 
7 uw „ Ein) Y+0 
Fam +) a + Em Dy—n + 1a], — — 


zu 1— cos - y SO y+n Em y+a 
Gy m+)) I 


Aus diesen Formeln ist leicht erkennbar, daß die Spannungen o, und r am Rande der Halb- 
ebene y = O identisch verschwinden, womit die Lastfreiheit nachgewiesen ist. 

Unsere Funktion muß aber noch eine Reihe anderer Bedingungen erfüllen. Schneiden wir 
aus der Halbebene ein beliebig geformtes Teil heraus, welches den Angriffspunkt einer Kraft P 
enthalten soll, so muß dieses Teil für sich im Gleichgewicht sein, d.h. das T 
Integral über Normal- und Tangentialspannungen längs des Randes muß 
gleich der im Innern angreifenden Kraft P sein. 

Die Berechnung dieses Integrals erfolgt am einfachsten mit Hilfe 
des Residuensatzes. 


A I b 
= N — > 1 u—— 1 = 
Ara 2rri Res (o),-.; 2ri,,,){@m+Di+m Er y 
Als einzige Funktion wird ctg (z2— aj) für z = aj singulär. Bild 2 
Nach alten Regeln ist das Residuum einer Funktion na an der Stelle gegeben durch a 
7 - 
cos 7 (.—aj) 1 h 
Res a — = — 
2 3 7 Ä 
sin „ (2@—-aj) 2 ee Mai 
z=(] 
jan ae = 2 sp 
Dane 5- (Bm diirm Yif=—-BAr- = 2 RE) 
(©) m m 
zufolge der Umformung 
ee N Bam el er 3.0 Be A zen5 
ij=ili 2\- a ar ’ a a 
Also: ZP, =—=P; > Pr=0, 


Damit ist gezeigt, daß die Singularität im Punkte y = a und wegen der Periodizität der trigono- 
metrischen Funktionen auch in den entsprechenden Punkten von der Art der Einzelkraft in 
Richtung der negativen y-Achse ist. 

Schließlich gilt es noch zu beweisen, daß die Verschiebungen &;n eindeutige Funktionen 
sind. Wollte man diesen Nachweis nach der üblichen Methode im Reellen führen, so wäre dies 
mit reichlich umständlichen Rechnungen verbunden, man betrachte nur die komplizierten 
Formeln für die Spannungen. 

Die hyperkomplexe Schreibweise läßt eine elegante Lösung zu. Die Verschiebungen sind 


gegeben durch m-+2 | ; 
el rel joR=-E+--, 
m+2 1 | 
nz el | o@=- Hint. 


Erstreckt man das Integral längs einer vorgegebenen Kurve, so erhält man die Differenzen aus 
den Verschiebungen der Endpunkte der Kurve. Betrachten wir insbesondere einen geschlossenen 
Weg, welcher eine Einzelkraft umzingelt, so müssen die Differenzen verschwinden, wenn & und 
eindeutige Funktionen sind. 


nr=E:2” |, 1 Pr #1 | F|=- nr u 
a P+J° 1 1 1 “ 
=» ee A DOR. SEI . SEIEN Tai? AM . .iepe ee —( RE 2: 
a m(m+1l) m-+I ni). J “ I 
‘ j2 
ln Ho) den 1-70 er. = nur gerade Potenzeninj. 


Nachdem die angegebene Funktion allen geforderten Bedingungen entspricht, stellt sie die strenge 
Lösung der Aufgabe dar. 
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Die Verteilung der Randspannungen ist gegeben durch: 


| een ge 
mg RR et: ü b 
yet n — ee Er — 
2bm im”? ı ml b (sin: ** + Ein m 
b b b b 
0 


Für da Flldel,d)=m,m= n haben wir die Randspannungen numerisch berechnet. 


Setzen wir das Verhältnis von Abstand der Kräfte vom Rande a zu Abstand der Kräfte 
untereinander b gleich A, so können wir das Abklingverhalten der Randspannungen bei wachsendem 


4 gut verfolgen, wenn wir (o,),-o und (0,),—o Bhc 
z=0 a—=b| 7) 
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Bild 5 
Um uns einen weiteren Einblick in die Lösung zu verschaffen, geben wir die Spannungen 
in verschiedenen ausgezeichneten Schnitten an: 
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I 
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Ein weiterer Symmetrieschnitt ist x = + 
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Überlagert man dem angegebenen Spannungszustand einen zweiten, der durch o,, = T3,, = 0; 
= —y gegeben ist, so wird Ga. für y> aNull. 


b 
Es gilt dann längs der Geraden x = + 5 


0 ’ 
Wir haben damit eine 1. Näherung für den Streifen 2 „ mit Einzellast 
— PA <T< 0/2 


im Innern, die gegenüber einer durch Reihenentwicklung gewonnenen Näherung den Vor- 
zug hat, daß sie in der Umgebung des singulären Kernes a F : genaue Werte liefert, die zudem 


von der strengen Lösung nur unerheblich abweichen dürften. 
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Die strenge Lösung selbst erfordert die Angabe eines Spannungstensors, ‚der (Oz)z=oj2 
selbst zu Null macht, ohne andere Randbedingungen zu verletzen. Wie man diesen findet, 
wurde in meiner Dissertation für den Fall a=0 gezeigt. Bei der Durchführung kommt man 
auf eine Fredholmsche Integralgleichung 2. Art, deren Lösung sich formal mit Hilfe der 
Neumannschen Reihe anschreiben läßt. Die numerische Durchführung dürfte, da man wegen 
der Darstellung der Singularität in geschlossener Form eine vorzügliche Konvergenz erwarten 
darf, keine grundsätzlichen Schwierigkeiten bereiten. 


Manuskripteingang : 10. 6. 1958 
Anschrift: Dr. Walter Engl, Karlsruhe, Germersheimer Str. 20 


Du 


Spannungsoptische Untersuchung des Querstoßes 
auf eine Kreisplatte*) 
Von H. Schwieger und V. Reimann **) 


Mittels des spannungsoptischen Reflexionsverfahrens wurde der zentrale Biegestoß auf eine am Rande 
freigestützte Kreisplatte untersucht. Dabei bestand die Kreisplatie aus zwei Schichten gleichen Glases und 
hatte in der Plattenmittelfläche eine Reflexionsschicht. Bei einseitiger Einstrahlung mit! inear polarisiertem 
Licht und durch eine geeignete photoelektrische Analyse des reflektierten elliptisch polarisierten Lichtes 
wurden die maximalen Torsionsmomente punktweise in Abhängigkeit von der Zeit bestimmt. Mit diesen 
Größen war es möglich, durch Integration die einzelnen Hauptbiegungsmomente und auch die Platten- 
durchbiegungen zuberechnen. Die experimentellen Ergebnisse wurden den aus der Theorie erhaltenen zum 
Teil gegenübergestellt. 


A photoelastic reflection method was used to investigate the effect of a central transverse impact on a circular 
plate simply supported at the edge. T’he plate consisted of two plates of the same glass cemented together with a 
thin layer of reflecting material beiween them. T'he used linearly polarised light beam entered the plate at one 
surface. After the light had transversed twice one half of the full plate thickness, the principal twisting mo- 
ments could be determined point by point in their dependence on time by analysing the reflected elliptically 
polarised light by means of a suitable photoelectric method. Hence the principal bending moments and the 
deflections of the plate could be computed by integration. T’he experimental results are compared partially 
with the theoretical ones. 


C’est au moyen d’un procede de reflexion photo-elastique que l’on a examine le choc de traverse central 
sur une plaque circulaire qui etait simplement appuyee le long du contour. La plaque circulaire consistait 
de deux couches adherentes du m&me verre et avait dans son plan moyen une couche mince r£flechissante. Un 
rayon lumineux polarise rectilignement entrait perpendiculairement a la surface de la plaque. Apres la 
lumiere a traversee deux fois la meme moitie de la plaque on a determin& les moments de torsion mazximaux 
point par point dependant du temps par une analyse photo-electrique appropriee de la lumi£re refletee et 
polarisee elliptiquement. A laide de ces resultats et par l’intögration on a pu caleuler chaque moment de 
flexion principal ainsi que les flexions de la plaque. Les resultats experimentaux furent opposes en partie ü 
ceux obtenus de la theorie. 


Ilpr IoMolıMm MeTona IA U3MepeHuA HANPAHKeHHÜ, OCHOBAHHOM HA HCCJIeNOBAHHH OITH- 
YECKOTO OTPArKeHNMA, ObLI UCCHIENOBAH HEHTPAABHLIH U3rNHÖamımmi ynap HA KPyToBylo INIACTHH- 
Ky, CBOÖ0]H0 Hasleraiolol HA ONOPY BIIONIb CBOeroO Kpaf. Ilıacrunka IpH 3TOM CocTonma 
u3 ABYX CJIOeB ONMHAKOBOTO COoPpTa CTeKa MH Conep;Kala B CpenHeli CBOeli IOBEPXHOCTH OT- 
pakarımä cıoA. Ilpm ONHOCTOPOHHEM OCBelleHMM JIMHEeÄHO IOJIAPH3OBAHHLIM CBETOM H IIPM 
HOMOIIM HAanerKameroO POTOIAIEKTPHYECKOTO AHAaım3a OTPArKaeMoOrTO YILIHNTHYECKH TOAIAPH- 
30BAHHOTO CBeTa ÖbLIM OIpelejIieHbI MAKCHMAABHLIE KPYTHJIbHbIE MOMEHTEI MIA OTNEIBHBIX 
TOYeK B 3ABUCHMOCTH OT BpeMeHu. IIpu IOMOLIM 9TUX BEJIHYMH IIYTeM HHTErPHPOBAHHNA YAAIOCh 
BbIYMCHUTb OTIEJIDHLIE TJIABHLIE MOMEHTEI U3THÖAHHA, A TAR>Ke U IIPOTHÖ INTACTUHkU. Pesy- 
JIBTATbI OIIBITA OTYACTH ÖbLIM CONOCTABJIEHLI C Pe3yAIbBTATaMM, HONYYCHHBIMM IIYTeM Teoperuye- 


CKUX Paccy;kıeHnd. 
1. Einleitung 

In der vorliegenden Arbeit wird über spannungsoptische Untersuchungen an einer zentral 
gestoßenen, elastischen Kreisplatte berichtet. Dabei handelt es sich um ein dynamisch-mechani- 
sches Problem, das sowohl in seiner experimentellen wie theoretischen Behandlung Schwierig- 
keiten mit sich bringt. Da eine homogene Platte bei Durchstrahlung mit linear polarisiertem 
Licht in Richtung ihrer Normalen keinen spannungsoptischen Effekt liefert, wurde von dem 
Reflexionsverfahren Gebrauch gemacht. Danach wurde die eigentliche Kreisplatte aus zwei 
Teilplatten spannungsfreien Glases zusammengesetzt. Vor dem Zusammenkleben der beiden 
Teilplatten erhielt die eine Platte einseitig eine Silberschicht und wurde so mit der anderen zu- 
sammengefügt, daß sich die reflektierende Schicht in der Plattenmittelfläche befand. 

Durchstrahlt man nun diese auf Biegung beanspruchte Kreisplatte einseitig mit linear 
polarisiertem Licht, so erhält man einen spannungsoptischen Effekt, da die Biegungsspannungen 
in der einen Plattenhälfte und die durch diese Spannungen verursachten optischen Doppel- 


x) Auszugsweise vorgetragen auf der Spannungsoptiker-Tagung in Leningrad, 13.—21. Febr. 1958 
und auf der GAMM-Tagung in Saarbrücken, 8.—12. April 1958. 2 s = ; 
**) Deutsche Akademie der Wissenschaften zu Berlin. 
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brechungen ihre Vorzeichen nicht ändern. Bei der Durchstrahlung einer homogenen Platte 
dagegen wird der in der einen Plattenhälfte erzielte optische Effekt wegen des Vorzeichenwechsels 
der Biegungsspannungen an der Plattenmittelfläche durch den in der anderen Hälfte auftreten- 
den Effekt kompensiert. 

Zur Messung der einzelnen maximalen Torsionsmomente an den verschiedenen Platten- 
stellen wurde eine photoelektrisch-polarisationsoptische Methode [1] benutzt, die schon bei vor- 
angegangenen Längsstoß- und Biegestoßuntersuchungen [2], [3] und [4] an Balken erfolgreich 
zur Anwendung kam. Unter Benutzung der maximalen Torsionsmomente ließen sich mittels der 
elastischen Gleichungen für die einzelnen Hauptbiegungsmomente sowohl die Plattendurchbiegung 
wie auch die Hauptbiegungsmomente selbst durch Integration ermitteln. Bei den Stoßunter- 
suchungen handelt es sich um den Stoß einer Stahlmasse auf die Mitte einer Kreisplatte, die an 
dem Rand freidrehbar gestützt war, wobei die maximalen Torsionsmomente insbesondere in den 
ersten Phasen der Biegestörungsausbreitung also kurz (bis ca. 6 - 10-*(s)) nach dem Aufprall der 
Fallmasse, in ihrem zeitlichen Verlauf punktweise gemessen wurden. 

Bei der theoretischen Behandlung dieses Stoßproblems geht man meistens von der ver- 
einfachten Differentialgleichung für die Plattenschwingungen aus, in welcher sowohl die Quer- 
kraftverformung wie auch die Rotationsträgheit der einzelnen Plattenelemente unberücksichtigt 
bleiben. Als Erregungskraft geht die Stoßkraft in die Schwingungsgleichung ein, wobei sie dem 
Betrage nach gleich der Trägheitskraft der abgebremsten Fallmasse ist und entsprechend der 
Wechselwirkung zwischen beiden Stoßpartnern, im Sinne der Hertzschen Theorie, von den 
elastischen Eigenschaften beider Körper und hier noch von den Schwingungseigenschaften der 
Platte abhängt. Die Stoßkraft selbst hat dabei einer Integralgleichung zu genügen, so daß unter 
Einbeziehung der Randbedingungen das Stoßproblem gewisse theoretische Schwierigkeiten be- 
reitet. Eine vereinfachte angenäherte Berechnung des Maximums der Stoßkraft und der Stoß- 
dauer ist leicht möglich, wenn man nach Eringen [5] die ‚‚collocation method“ benutzt, wobei 
man für kleine Massenverhältnisse von Fallkörper und Platte die Stoßkraftkurve durch eine 
einfache trigonometrische Funktion annähert. Die nachfolgende kurze Wiedergabe der Theorie 
verfolgt das Ziel, die Berechnungsgrundlagen für obige Stoßgrößen wiederzugeben. Gleichzeitig 
werden die Eigenwerte für freie Plattenschwingungen, die wegen der Randbedingungen von der 
Poissonschen Konstanten » (< 0,5) abhängen, für » = 0,25, 0,30 und 0,35 mitgeteilt. Aus den 
Versuchsergebnissen werden dann die Querkräfte kurz nach dem Stoß in Stoßstellennähe abge- 
schätzt und mit den nach der ‚collocation method‘ berechneten gegenübergestellt. 


2. Theoretische Zusammenhänge 
Beschränkt man seine Betrachtungen auf ein einzelnes Plattenelement, dessen Mittelfläche 
mit der x, y-Ebene bzw. r, g-Ebene eines gewählten Koordinatensystems zusammenfallen soll, 
so herrscht unter Vernachlässigung von Querkraftdeformation und Rotationsträgheit in z-Rich- 
tung in jedem Moment t ein Gleichgewicht zwischen der elastischen Kraft B- AAw, der Träg- 
heitskraft u w und der auf die Flächeneinheit bezogenen äußeren Belastung q: 


B An 1a EU ee ne pe ee kenne (1). 


Hierbei sind w = w(z, y, ft) (cm) die Durchbiegung der Plattenmittelfläche, A = 02/02? + 02/dy? 

bzw. 02/ör?2 + 1/r - 9/Or (für den rotationssymmetrischen Fall), u die Plattenmasse pro Flächen- 
3 

einheit, w = 92w/df? (cm/s), B = 1/12. ( - er (kp - cm) die Biegungssteifigkeit mit h(cm) der 

Plattendicke, E (kp/cm?) dem Elastizitätsmodul und » der Poissonschen Konstanten. Bei einem 

Einzelstoß ist qg nur an der Berührungsstelle vorhanden. Dabei ist die örtliche Verteilung von q 

an der Berührungsstelle unbekannt. Jedoch muß das Integral über die Verteilungsfläche gleich 


der Stoßkraft sein, so daß sich ergibt: 
FOÖ= [aa y,)dF=dqd/dF, 4 == Rlattenflache er se (02), 
A ) 


Hierbei denkt man sich anstelle von g als Ortsfunktion diese gleichmäßig über eine kleine Fläche 
f verteilt und die Stoßkraft selbst durch einen Grenzübergang gewonnen, wonach bei kleiner 
werdender Fläche die Größe q zunehmen soll. Bei dieser Grenzbetrachtung ist die Stoßkraft: 


Eee tnae es a. 
1>0 
9>o 


Zum anderen ist die Stoßkraft dem Betrage nach gleich der Trägheitskraft der abgebremsten 
Fallmasse. Es gilt mit w, als Verschiebung des Schwerpunktes der Fallmasse m 
dw, 


mia EO een nen 
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und nach zweimaliger Integration erhält man 
t(t 
mu —S||Fodeldr +mot ee (5) 
ob 


oder umgeformt 
t 


w, ni |» Fe) ee (6). 


v, ist hier die Auftreffgeschwindigkeit der Fallmasse. 

Beim Zusammentreffen der beiden Stoßpartner erfahren sie jeweils eine Zusammendrückung 
in der Nähe ihrer Berührungsflächen, so daß sich die beiden Schwerpunkte von der Fallmasse 
und dem Plattenelement an der Stoßstelle gegenseitig nähern. Ist w, = w(r,, 9) die Platten- 
durchbiegung an der Stoßstelle, so ist die gegenseitige Annäherung | w,— w,. Nach der Hertz- 
schen Stoßtheorie steht diese Annäherung mit der Stoßkraft in folgendem Zusammenhang: 

w—w,=kF(%B bzw. w=KFOP+W,....- 00. oil: 


Die Hertzsche Konstante k ae sich aus: 
BR — =D, +D,)- (a Fa)" 
2 
wobei 
D, = ME) A—HWM), D= (4) (l—%) 
und r, bzw. r, die Krümmungsradien der beiden Stoßflächen sind. 
Verknüpft man Gleichung (6) mit Gleichung (7), dann erhält man: 


kurlDae— in j dr Fi) dem ua. 1 (8). 
ö 


Diese Gleichung wird vollständig, wenn man die Plattendurchbiegung w, da sie von der Stoß- 
kraft und den Schwingungseigenschaften der Platte abhängt, mit diesen Größen in Verbindung 
bringt. Ohne die theoretische Ableitung, die von Timoshenko [6] unter Benutzung von Glei- 
chung (1) ausgeführt wurde, im einzelnen wiederzugeben, ist schließlich 


t 
2 ; 
= >, | (Lan) FOisn nd Hi.. TR, 
und somit 
t 
1 2m c x 
k FR = v,t— jr |r@ len ur; en SONO E—Didr . . (10). 
1) m,n 
Im einzelnen day. für eine zentral gestoßene Kreisplatte (r = 0): 
Cno = Imo (0,0) . a m=1,2,3,..., n=0 (endliche Lösung für r = 0), 
ar! p Pm @ \ Jı(Pm 4) 
by =20a2 J? m - —2| de *) 2 a ne M 
0 ee) T (e2)5 Ilpm al’ 
9mo(0,0) = J,(0) — J(Pm @) I,(0) Eigenfunktion, 
Io(Pm 4) 


J, und Jı: Besselsche Funktionen 1. Art Nullter und Erster Ordnung mit reellem Argument, 
I,und I,: Besselsche Funktionen 1. Art Nullter und Erster Ordnung mit imaginärem Argument 


: 0: h\ı12 j 
Dn = | B (em 2): Quadrate der Eigenwerte, 


%m(s): Eigenfrequenzen, 

o (g/cm?): Dichte der Platte, 
a (cm): Plattenradius, 

M (g): Plattenmasse. 


Die Eigenwerte p,„ müssen dabei aus folgender Frequenzgleichung bestimmt werden: 
JPm a) , I(Dm a). 2 Pmm A 


Js(Pm ©) Tlpae) "15 Ver 11). 


ö 
> 
- 
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Diese Frequenzgleichung gilt insbesondere für den Fall, wie er zunächst von uns untersucht 
wurde, daß die Kreisplatte am Rande frei drehbar gelagert ist. Durch die Randbedingungen 
kommt die Poissonsche Konstante ins Spiel, so daß die Bestimmung der Eigenwerte lang- 
wieriger als bei Kreisplatten mit eingespannten Rändern ist und sich somit für jeden Unter- 
suchungsfall etwas verschiedene Stoßbedingungen ergeben. Für v = 0,25, 0,30 und 0,35 be- 
kommt man folgende (p,, a)-Werte und c„o-Werte: 


Pm @-Werte 
m 1 2 3 | 4 
(Pros | 232045 5,4462 8,6082 | 1,7586 
(0), 00 Anz 5,4516 86118 | 11,7607° 
(a or 2,2378 508 | SE 11T 
m 5 6 7 8 
(Pu Y),-o | 149051, | 18,0498 21,1936 24,3368 
De 14,9068 18,0512 | 21,198 24,378 
(Pd 10086 | 18,0527 | 21,100  24,3380 
Cmo- Werte 
m l 2 3 4 
(an 1,7460 4,3290 6,7883 9,2598 
en no, Ang“ | 6,70 NL 9,201 
m | 5 | 6 m | 8 
ee: 11,7389 14,2942 16,7340 19,1500 
nos else, 108er 


Mit diesen Größen wird es möglich, nach der ‚„collocation method“ das Maximum der Stoß- 
kraft und die Stoßdauer zu berechnen. 

Führt man zum Zwecke der Berechnung dieser speziellen Stoßgrößen dimensionslose Ver- 
änderliche und verschiedene Abkürzungen ein, so geht Gleichung (10) in folgende Form über: 


4 
FAR DPESTFT- ALT KEN An) it —4SO . x Den). 
0 
Dabei sind 
L=tIsn-= te: T, willkürliche Bezugszeit, £, bezogene Stoßdauer 
a a A rel brndhit 


A®X(&) = (Tofmv,)-FET) mit X) =Xl)=0, AK) =1, 


Ken) = (Cm) +2 mim) I)" sin on Cm). 


m 


Zur Berechnung der maximalen Stoßkraft differenziert man obige Gleichung und bekommt unter 


der Beachtung, daß im Maximum der Stoßkraft E- ke, ® — ( ist, 
[AS (2) _ cr def 
A® = ı ——— = BEE Ce LE (Lil); 
a a A e)) (13) 
woraus sich, eingesetzt in Gleichung (12), die halbe bezogene Stoßdauer 
Be SCAN nenn N 


ergibt. Nun liegen verschiedene Berechnungen darüber vor, durch welche Funktionen man einen 
symmetrischen Stoßkraftverlauf am besten annähern kann. Benutzt man die als geeignet be- 
15 
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fundene Funktion X(£) = sin z & ‚ dann wird 
8 


sm - Elan +2 HR Te EL. eure 
sam = ll t2luRg. RER A 
mit is 
Rı = Cmo (1 An) (aut sin "ie — 1), 
Be, To = . 
ag 1%; Sn T=o, 
IT 


50 


40 


w 
b>} 


Hilfsgrößen R, und Rz 
n 
oO 


Bild1. Die Hilfsgrößen R, und AR, in Abhängigkeit vom Eigenwertparameter 4, 


In Bild 1 sind nun R, und R, in Abhängigkeit von A, für eine freigestützte Kreisplatte aufgetragen. 
Mit diesen Hilfskurven und den Gleichungen (13) und (14) ist es leicht möglich, die Stoßdauer und 
das Maximum der Stoßkraft zu berechnen. Diese Berechnung wird brauchbare Werte liefern, 
solange das Massenverhältnis m/M klein bleibt, da dann die Stoßkraftkurve etwa einen symme- 
trischen Verlauf zeigt. Bei größeren Massenverhältnissen wird der Verlauf zunehmend unsym- 
metrischer und kann mehrere Maxima haben. Für diese Fälle versagt die „collocation method“ 
in dieser Form und auch eine schrittweise Berechnung der Stoßkraftkurve nach Timoshenko [7] 
wird langwieriger. 


3. Versuchsdurehführung und Ergebnisse 


Bevor mit den eigentlichen Messungen der maximalen Torsionsmomente begonnen wurde, 
wurde zunächst festgestellt, inwieweit das elastische Verhalten der zweischichtigen Kreisplatte 
mit dem einer homogenen Platte übereinstimmt. Dazu wurden an der am Rande frei gestützten 
Kreisplatte die Durchbiegungen in der Mitte für verschiedene Stufen der ebenfalls in der Mitte 
aufgebrachten Einzellast gemessen. Das Ergebnis dieser Messung findet man in Bild 2. Hierbei 
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wurden die gemessenen Durchbiegungen mit der 3. Potenz der jeweiligen Plattenhöhe multi- 
pliziert, um eine Steigung der Geraden zu erhalten, die nur von den elastischen Materialgrößen 
abhängt. Bei genau gleichem elastischem Verhalten beider Platten müßten beide Geraden die 
gleiche Steigung haben. Die Abweichungen können nun darauf hinweisen, daß doch einige Unter- 
schiede im Biegungsverhalten vorliegen. Jedoch ist dabei zu beachten, daß die Platten nicht 
durchweg jeweils die gleiche Höhe hatten und hier mit einer längs des Plattenumfanges gemittelten 
Plattenhöhe gerechnet wurde. Diese Höhen können durchaus von den über die Plattenfläche ge- 
mittelten abweichen und, da sie in der 3. Potenz eingehen, größere Steigungsunterschiede ergeben. 
Im folgenden wurde nun angenommen, daß die Platten in ihrem Biegungsverhalten nahezu gleich 
sind. Zur Erhärtung dieser Annahme, wurde aus dem vorne angegebenen Eigenwert die Eigen- 


4 3,0 frequenz der Grundschwingung berechnet, diese der polari- 
Pr sationsoptisch-elektrisch bestimmten und der aus der experi- 
; mentell ermittelten Zeit bis zur maximalen Durchbiegung 
' abgeschätztenEigenfrequenz gegenübergestellt. 
$; 24 | 


o]] 


| 


Hg. 


Bild 3. Schema des polarisationsoptischen Ver- 
suchsaufbaues zur Untersuchung des Platten- 
“ biegestoßes; 


3 BER ( kp) 9 EM: Elektromagnet, FM: Fallmasse, Pl: Platte, 
b Ag: Silberschicht, A/4-Bl.: Viertelwellenlängen- 
Bild 2. Die mittige Durchbiegung bei verschie- blättchen, P: Polarisator, A: Analysator, Sp: Spie- 
denen Laststufen, gel, L: Linse, Zg,.L.: Quecksilberhochdrucklampe, 
©: für eine homogene Platte, F: Filter (A = 546 ma), Ph: Photozelle; Einfalls- 
°; für eine Zweischichtplatte winkel des polarisierten Lichtes ca. 1,5° 


Dem errechneten Wert von 261 (Hz) stehen. experimentell gemessene von ca. 275 (Hz) und 
278 (Hz) gegenüber. 

Zur punktweisen Bestimmung der zeitlich sich ändernden maximalen Torsionsmomente 
beim Plattenbiegestoß wurde folgende, in Bild 3 festgehaltene Anordnung verwendet. Die Kreis- 
platte (M = 2,94 (kg), oe = 2,5 (g/cm?)) wird längs eines Kreisringes mit einem Durchmesser von 
40 (cm) freidrehbar gestützt. Eine Fallmasse aus Stahl mit einer kugelig abgedrehten Stoßfläche 
(Krümmungsradius = 5 (cm)) und einer Masse von m = 200 (g) fällt nach dem Abschalten 
eines Elektromagneten auf die Mitte dieser Platte. Dabei ist die Höhe des Elektromagneten 
so eingestellt, daß die Masse eine Fallhöhe von 4 (cm) durchfällt, was einer Auftreffgeschwindig- 
keit von v, = 88,95 (cm/s) entspricht. 

Bei der Berührung der Fallmasse mit der Kreisplatte wird sie selbst abgebremst, wodurch 
eine Stoßkraft entsteht, deren Maximum um vieles größer (Faktor 500—1000) ist, als das Gewicht 
des Fallkörpers ausmacht. Durch diese Stoßkraft wird die Platte durchgebogen; es entsteht also 
ein Biegespannungszustand, der sich zu den Plattenrändern hin schnell ausbreitet. Auch ent- 
stehen in der Stoßstellennähe gegenüber einer statischen Belastung mit diesem Gewicht viel 
größere Querkräfte. 

* Durchstrahlt man nun bei der gestoßenen Platte die untere Hälfte mit linear polarisiertem 
Licht, so wird beim Eintritt in die nunmehr optisch doppelbrechend gewordene Glasplatte die 
einfallende linear polarisierte Lichtwelle in zwei Lichtkomponenten aufgespalten. Die Richtungen 
dieser Lichtkomponenten fallen dabei mit den Richtungen der Hauptbiegungsmomente zusam- 
men. Unter der Annahme, daß diese Richtungen über die halbe Plattendicke (h/2) konstant sind, 
addieren sich für jede differentielle Plattenschicht die durch die optische Doppelbrechung hervor- 
gerufenen optischen Phasenunterschiede dieser Lichtkomponenten angenähert algebraisch. 
Nimmt man ferner an, daß auch bei einer dynamischen Beanspruchung die Biegungsspannungen 
linear von der Plattenmittelfläche aus ansteigen, dann beobachtet man einen resultierenden opti- 
schen Gangunterschied, welcher der Spannung in der Mitte (h/4) zwischen Oberfläche und Plat- 
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Bild4. Schema der photoelektrischen Registriereinrich- 
tung zur Messung der optischen Phasendifferenzen; 


EM: Elektromagnet, FM: Fallmasse PI: Platte, KO: Zwei- 

kathodenstrahl-Oszillograph, Ph: Photozelle, V: Breit- 

bandverstärker, ZA: Gerät zur einmaligen Ablenkung der 
Elektronenstrahlen, 7’G: Tongenerator 


Sı- Kurven S2p-Kurven 


r=7,5cm 


Einige Beispiele aufgenommener Oszillogramme zur Bestimmung der optischen Phasendifferenzen; 
8/-Kurven: ohne phasenändernde Wirkung des A/4-Blättchens; 
5>-Kurven: mit einer zusätzlichen Phasenverschiebung von d= + n/2 bzw. 4° = + 90° 
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tenmittelfläche proportional ist. Bringt man die Polarisationsfilter, Polarisator und Analysator, 


in gekreuzte Stellung und ihre Schwingungsrichtungen unter 45° zu den jeweiligen Haupt- 
momentenrichtungen, dann erhält man bei einem zeitlich veränderlichen Spannungszustand 
eine Intensitätsänderung, die proportional zu sin? A/2 ist. Hierbei ist A die oben erwähnte 
optische Phasendifferenz der beiden Lichtkomponenten. Diese optische Phasendifferenz steht 
mit der mittleren maximalen Schubspannung 7,,.z und dem Durchstrahlungsweg 2 - h/2 und der 
wellenlängenabhängigen spannungsoptischen Konstanten K, in folgendem Zusammenhang: 


Al = BR du. K ne 22 De te De ma en (17): 
Zum weiteren ist 2 Tyar = m ‚ so daß sich bei Kenntnis von A schließlich das maximale 
(auf 1 (cm) Schnittlänge bezogene) Torsionsmoment wie folgt berechnen läßt: 
h 
T eu A (a a AS Ä 
mul) A (d 360° .6 K, s (18) 


Die Bestimmung von 4 erfolgt nun polarisationsoptisch-elektrisch (siehe [1]). Dazu bringt man 
zusätzlich ein A/4-Blättchen in den Strahlengang zwischen Polarisator und Analysator. Am 
besten wird es vor dem Analysator aufgestellt. Liegen die Schwingungsrichtungen des A/4-Blätt- 
chens parallel zum Polarisationskreuz (1. Stellung), so ist bei veränderlichem A die Intensitäts- 
änderung hinter dem Analysator ebenfalls proportional zu sin? A/2. Dreht man dagegen das 
ı/4-Blättchen um 45° (2. Stellung), so daß also die Schwingungsrichtungen desselben parallel zu 
denen des optisch doppelbrechenden Plattenelementes (also parallel zu den Hauptmomenten- 
richtungen) verlaufen, dann ist bei einer Änderung von A die Intensitätsänderung hinter dem 
Analysator proportional zu (sin? (A/2 + 45°) — 1/2). Diese Intensitätsänderungen werden nun 
photoelektrisch gemessen, wozu eine Photozelle in den an der Silberschicht reflektierten Licht- 
strahl hinter dem Analysator gebracht wird. Diese ist wiederum mit einem Breitbandverstärker 
und einer Braunschen Röhre verbunden. Auf diese Weise werden die Intensitätsschwankungen 
in proportionale Stromschwankungen umgeformt und die am Photozellenwiderstand entstehen- 
den Spannungsänderungen verstärkt den Meßplatten der Braunschen Röhre zugeführt. 

Eine Übersicht über die elektrische Anordnung findet man in Bild 4. Dabei wird die Ab- 
lenkungsschaltung für die Elektronenstrahlen durch den Kurzschluß zwischen Fallmasse und 
einer an der Stoßstelle aufgeklebten dünnen Aluminiumfolie (ca. 15 # dick) ausgelöst. Zur zeit- 
lichen Festlegung der Zeitablenkung der Elektronenstrahlen des Kathodenstrahloszillographen 
(mit zwei Ablenksystemen) dient hier ein Tongenerator. Die Sinus-Spannungen des Tongene- 
rators werden an die 2. Meßplatten gelegt, so daß man neben dem zeitlichen Verlauf der Inten- 
sitätsänderungen die Sinusschwingungen beobachten kann. Die Oszillogramme werden dann zur 
Bestimmung der optischen Phasendifferenzen mit einer Kleinbildkamera photographiert. Einige 
Versuchsbeispiele, bei denen sich das A/4-Blättchen in der 1. bzw. 2. Stellung befand, sind in dem 
Bild 5 festgehalten (Tonfrequenz 20 (kHz)). Aus den für jede Meßstelle in Abhängigkeit von 
der Zeit sich ergebenden Ausschlägen s, und s, errechnet sich die optische Phasendifferenz aus: 


10 32 = S1/s, Keindeutig für AS 180°), .n . sulor an (19). 


Zur Untersuchung der Kreisplatte an den einzelnen Meßstellen wurde längs des Durchmessers die 
Plattenoberfläche mit schwarzem Papier abgedeckt, in welchem sich zur Ausblendung der Meß- 
punkte kleine Öffnungen (2 ca. 0,28 (cm)) befanden. Die einzelnen Blendenöffnungen werden mit 
dem Lichtstrahl nach und nach ausgeleuchtet, indem man die ganze optische Anordnung parallel 
zum Durchmesser verschiebt. 

Vor der eigentlichen Messung werden zunächst die Richtungen der Hauptbiegungsmomente 
an den verschiedenen Meßstellen festgestellt. Dazu dreht man das Polarisationskreuz schrittweise 
so lange, bis man beim Stoß keinen Ausschlag am Kathodenstrahloszillographen beobachtet. In 
diesem Fall fällt dann die Schwingungsrichtung des Polarisators mit einer der Richtungen der 
Hauptbiegungsmomente zusammen und es findet demzufolge keine Komponentenzerlegung statt. 
Im Anschluß daran stellt man das Polarisationskreuz unter 45° zu diesen Hauptmomentenrich- 
tungen und photographiert in der vorne geschilderten Weise die s, und s,-Kurven. 

Damit der Elektronenstrahl in seiner Ausgangsstellung nicht das Bild überstrahlt, wird 
gleichzeitig mit dem Abschalten des Elektromagneten eine Helligkeitsschaltung betätigt. Da- 
durch wird die elektrische Spannung am Wehnelt-Zylinder der Braunschen Röhre so verändert, 
daß erst kurz vor dem Stoß die Stärke des Elektronenstrahls und damit des Lichtpunktes auf dem 
Schirm der Braunschen Röhre auf die gewünschte Größe ansteigt. Ferner ist es noch möglich, 
nach der in Bild 6 dargestellten Weise die mittige Plattendurchbiegung zu messen. Dazu wird 
an der Plattenunterseite eine dünne Rasierklinge angeklebt, die beim Stoß auf eine Meßspindel 
stößt. Durch die Berührung wird eine Gleichspannung an die Meßplatten des Kathodenstrahl- 
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FM 


Bild 6. Elektrische Anordnung zur Messung der 
Durchbiegung in der Plattenmitte; 

FM: Fallmasse, Pl: Platte, KM: Kontakt-Mikro- 
ımeterspindel, KO: Kathodenstrahl-Oszillograph, 
U: Gleichspannung, ZA: Zeitablenkung, 
tp: Durchbiegungszeit 


EikunerT 
ann 


d), 
Die aus den s,- und sg,-Ausschlägen berechneten maximalen Torsionsmomente an den verschiedenen Meßstellen 
in Abhängigkeit von der Zeit 
a) r=1(cn); b)r=25(cm); o)r=ö5(em); d) r=17,5 (cm) 


Bild 7. 
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oszillographen gelegt, wodurch zwischen Beginn der Zeitablenkung und dem Einsatzpunkt des 
1. elektrischen Kontaktes die Zeit (£,) für eine bestimmte Durchbiegung abgelesen werden kann. 
Durch Verändern des Abstandes zwischen Rasierklinge und der Meßspindel läßt sich die Durch- 
biegung bis zum Maximum in Abhängigkeit von der Zeit gewinnen, wie sie im Bild 13 wiederge- 
geben ist. 

4. Auswertung 


Aus den an jedem Meßpunkt gewonnenen s,-Kurven und s;-Kurven wird nach Gleichung (19) 
zunächst für ein Maximum die optische Phasendifferenz berechnet. Da man nun weiß, welche 
optische Phasendifferenz A,., zum maximalen Ausschlag (s;)nı. gehört, kann man für die ver- 
schiedenen anderen Punkte der s,-Kurve das Verhältnis s5/(S3)mar bilden, das dem Ausdruck 


sin? (A/2 + 45°) — 1/2 
ins Mn + 45°) SE iR Er olalgeeee BZ Pod a oe 
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Bild 8. Die maximalen Torsionsmomente an weiteren Plattenorten 
a) r=10(cm); _b) r=125(cm); c)r=15(cm); ..d) r = 17,5 (cm) 
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i i ießli ie j öri i Phasendifferenz zu be- 
roportional ist, um schließlich die jedem s,-Wert zugehörige optische Zi 
hen Dazu zeichnet man sich die sin? A/2 Kurve auf und wertet die eben genannte er 
am besten graphisch aus. Auf diese Weise kann man die ganze s,-Kurve bezüglich ihrer optischen 
Phasendifferenzen festlegen und letztlich, da 


- h z 10-2. 40 

To = 6 K, 360° 4° = 8,26 . 10-2. 4° (kp), ee 
cm 
kp 
ist, das maximale Torsionsmoment in Abhängigkeit von der Zeit erhalten. 


(Kr = 0,005 ( ). h= 0,893 m) 


Bein INRr 
“ 4 Se] Nu 5 
‚LI 
REN 


Bild 9. Der örtliche Verlauf der maximalen Torsionsmomente zu verschiedenen Zeiten 


8 20 12 14 16 r(em)18 
Bild 10. Die Abhängigkeit des 1. Maximums der in Bild? und 8 fest; 
T pay" Kurven vom Plattenort 


gehaltenen 


Eine Zusammenstellung der an verschiedenen Plattenstellen erhalte 
findet man in den Bildern 7 und 8. Diese Auswertung ist zu empfehlen, da 
nach Gleichung (19) für die kleinen Ausschläge zu großen Fehlern Anlaß g 
Bildern 7 und 8 wiedergegebenen Ergebnisse wurden so umgezeichnet, da 
Zeiten die maximalen Torsionsmomente längs des Radius erhält. Auf 
Ausbreitungsverhältnisse des Biegespannungszustandes besser sichtbar 
Bild 9 zu entnehmen ist. Ferner ist zur weiteren Veranschaulichung in Bild 10 die Abnahme des 
1. Maximums von T,,., mit zunehmender Entfernung vom Plattenmittelpunkt festgehalten. 
deren Richtungen parallel bzw. 
schlagen worden. Bei der Kreis- 


nen Torsionsmomente 
die Berechnung von A 
eben kann. Die in den 
B man für verschiedene 
diese Weise werden die 


EI WETE 
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platte stehen für den rotationssymmetrischen Biegespannungszustand die einzelnen Haupt- 
biegungsmomente mit den Ableitungen von w in folgendem Zusammenhang: 


row nn 
M=M=—B(v.% Be a a on (22), 
OB av a 
M =M, =— Tan. a En 5 (23), 
Eh: 
B= PA — 4,152. 10*(kp- cm), (E = 6,56 10° (kp/cm?), v»— 0,25, h= 0,893 (cm)) (24). 
Bildet man die Differenz der Momente, so erhält man bekanntlich das maximale Torsionsmoment 
9:w 1 ow 
Tu ME er © (25), 
wobei man nach Umformung der Differentialausdrücke setzen kann 
o (1 0ow 
Ba—mrZ |) 2 Inu se 9 er WR . (26) 
oder 
o/1 ow 2, I max Beer ER 
A ee ae en 


Diese Gleichung bildet nun den Ausgang, um schrittweise durch zweimalige Integration längs des 
2 i j : 1b 
Radius die Plattendurchbiegung zu erhalten. Im weiteren Verlauf zeichnet man C- ve für 


t = const. in Abhängigkeit vom Radius auf Millimeterpapier auf, errechnet zwischen bestimmten 
Grenzen, nämlich zwischen einem beliebigen Plattenpunkt und dem Plattenrand den Flächen- 
inhalt, um dann 


a 

1 ow 1 ow ie e. BT 

e ne) = Zl — c. | = dr (C=0642.104&p-tem) . ... (28) 
a 

zu erhalten. Führt man die Integration zunächst nur für solche Zeiten aus, in welchen die Biege- 


Ä al ‚28 10 
wellen den Plattenrand noch nicht erreicht haben, so verschwindet in diesem Fall a = N 
Es ist somit 2 


1 dw Aue 
ae 1 2 5 Eee  SRSEEE 2 
u =-c-| r (29) 
oder 
ow f Tone 
dr = rc. | Iawar rs FE (30), 
woraus sich nach nochmaliger Integration schließlich die Plattendurchbiegung 
[Times 
wir, 9=C-| | ne ir + ua en 7 Were (31) 


errechnet. Hier ist, solange die Biegungsstörung den Plattenrand noch nicht erreicht hat, w(a,) =0. 
Die auf diese Weise gewonnenen Plattendurchbiegungen sind für verschiedene Zeitpunkte in dem 
Bild 11 festgehalten. Dies sind die Durchbiegungen, die der Wirkung der Biegungsmomente zu- 
zusprechen sind. Der Durchbiegungsanteil, welcher von der Querkraft hervorgerufen wird, kann 
also durch diese Auswertung nicht erfaßt werden. Unter Verwendung vorstehender Differential- 
ausdrücke von w lassen sich die einzelnen Hauptbiegungsmomente wie folgt 


a 


Et, | an dr + Tonelt, ) Pa ER. 2). 


r 


Mind =.B.C. 


m = mc ea... 0 
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MAD) =; x 6 +») | u dr + Tyaz(t, 0 wege: STRRERAN 


nn = la tn | Teldarrruue | 


aus den spannungsoptisch bestimmten Torsionsmomenten berechnen. Zur Veranschaulichung 
sind mit » — 0,25 die Hauptbiegungsmomente für die Zeit = 0,5 - 10”* (s) und != 0,75 - 10% (s) 
in Bild 12 festgehalten. 


W: 10* (cm) 
10 


f= 0,75:10°%(s) t= 05-1015) 


Bild 11. Die durch Integration erhaltenen Plattendurchbiegungen zu verschiedenen Zeiten 


Die Abschätzung der Biegungsmomente M, und M, im Stoßpunkt, in welchem im allge- 
meinen die Momente am größten und gleich groß sind, erfolgt nun in der Weise, daß man die 
Wer — 20 + Wo 
Ar? 
unterlagen (für r = 1 (cm)) entnimmt. w; ist hier die Durchbiegung unter der Stoßstelle und 
W; +, und w;_, sindim Abstand Ar benachbarte Durchbiegungen, die hier wegen der Symmetrie 
gleich groß sind. Für verschiedene Zeiten erhält man mit » = 0,25 


02 1 
Krümmung 52 durch ersetzt und den Wert = = aus den Integrations- 


2 

[= 0,5 10-469): En = 3,3. 10-+ (em); M, = 17,5 (kp), 
1: om 27. 
3, = 35. 104 em); M, = 18.(kp), 
ww 


t = 0,75 10 (s): Ba 5. 10-*(em-1); M,= 26,1 (kp), 


Pia 5,1 10% (cm); M,= 26,4 (kp). 

Ferner kann man die Querkräfte abschätzen, wenn man die für eine statische Belastung bestehen- 
den Beziehungen hier als angenähert gültig ansieht. Bei einer Kreisplatte steht die Querkraft 
(bezogen auf 1 (cm) Schnittlänge) mit den Biegungsmomenten in folgendem Zusammenhang: 


Kr 
r 


“ 
« 


a un h ai 2 a 


e 

# 
u: 
x 

= 

e 
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Für, den Kreisschnitt von r = 1 (cm) ergeben sich aus den Versuchsergebnissen 
= 0,5 *10=4 (5): 'Q,, = 12,4 (kp - em), 


0, = 277: Qr, = 78 (kp), 
2L=,0,75 107% (68) 0, = 15,6 &p-cm , 
0, = 277: Q, = 9% (kp). 


Setzt man diese Querkräfte den in den jeweiligen Augenblicken herrschenden Stoßkräften gleich, 
so erkennt man, welche beachtliche Wirkung diese Fallmasse beim Biegestoß hervorrufen kann. 
Zur größenordnungsgemäßen Gegenüberstellung mit obigen Werten wird zunächst die Stoßdauer 
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Bild12. Der Verlauf der Hauptbiegungsmomente M, = M, und M,.=M 
zur Zeit t= 0,5 .10-?(s) und i = 0,75 « 10% (s) 


r 


berechnet. Dazu gibt man einige £,-Werte vor und berechnet mit der aus dem Eigenwert p, be- 
stimmten Eigenkreisfrequenz a, = 1,641 - 10? (s-!) die zugehörigen Eigenwertparameter A,. Aus 
dem Bild 1 entnimmt man dann die R,- und R,-Werte und versucht mit den Gleichungen (15) 
und (16) die Gleichung (14) zu befriedigen. Auf graphischem Wege erhielt man mit a* =5,704- 10° 
(k = 0,7341 - 10 (p=2? cm), v, = 88,59 (cm/s), m = 200 (g), T, = 1(s)) eine bezogene Stoß- 
dauer von £, = 0,92 - 10%. Schließlich errechnet sich mit £, = 0,92 - 10-* die Größe S’ (£,/2) zu 
1,179 - 10-* und daraus A® = 0,8482 . 10%. Als Endergebnis bekommt man ein Maximum der 


Stoßkraft von Ft) = A®- Ben 58 (kp). Bei den Versuchen ist also die Stoßdauer 


länger als sie nach der Theorie berechnet wird und liegen die Querkräfte in der zu erwartenden 
Größenordnung. Diese können für einen anderen Zeitpunkt als die hier herausgegriffenen Mo- 
mente noch größer sein. Ein Unterschied ließe sich dadurch erklären, daß die Aluminiumfolie 
den Stoß gegenüber einer unmittelbaren Berührung zwischen Stahl und Glas ‚weicher‘ macht. 
Auch sind in der Theorie die Rotationsträgheit und die Querkraftverformung der Plattenelemente 
nicht. berücksichtigt worden. (Innerhalb der Stoßdauer ist aber gemäß Bild 13 die von den 
Querkräften hervorgerufene Plattendurchbiegung bald von gleicher Größe wie die von den Bie- 
gungsmomenten hervorgerufene Durchbiegung.) 
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Die maximale Durchbiegung in der Plattenmitte wurde hier zu Wnaz = 1,62 10-2 (cm) 
bestimmt. Dieser Wert stimmt gut mit dem nach 


Wa; = 2 )p 


6 7 8 9 10 -10*(s) 

Bild 13. Die mittels der Anordnung in Bild 6 gemessene mittige Plattendurchbiegung in Abhängigkeit von der Zeit. (Die Zeit bis zum 

Maximum entspricht etwa 44 Schwingungsdauer.) Die Werte x sind durch Integration erhalten und der Wirkung der Biegungs- 
momente zuzusprechen 


berechneten überein. Vorstehender Ausdruck wurde von Eschler [8] für kurzzeitige Stöße aus 
dem Impuls- und Energiesatz abgeleitet, worin die Größen folgende Bedeutung haben: 


G„: Fallgewicht, u : Massenverhältnis, 


g: Erdbeschleunigung, 
16x ,(l-+») 

6: ou B +») 

RR ı nach Szab6 [9] für eine 

"484+9) freigestützte Kreisplatte. 

a 
Ey) 


v,: Auftreffgeschwindigkeit, 


Zum Schluß sei noch mitgeteilt, daß zur Kontrolle auf elektrischem Wege die radiale Dehnung g, 
an der Stelle r = 4 (cm) gemessen wurde. Für { = 0,75 - 10—: (s) stehen sich gegenüber: 


Elektrisch: e&, = ca. — 0,4 : 10%, 


Optisch: &, = — 0,45 - 104, 


2. we. 
6 = hj2 en + De Ci | mie dr + C- Tas N) 


r 
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Berechnung des Spannungsverlaufs in Rohrverzweigungen 
Von H. Neuber 
Für Rohrverzweigungen wird ein Verfahren zur Berechnung der an den gefährdeten Stellen auftretenden 
festigkeitsmäßigen Beanspruchung angegeben. 


A method is presented which for pipebranchings in the form of Y-pipes allows the calculation of the strain 
which occurs at the endangered points. 


Pour des branchements de tuyaux est indique un procede pour la calculation des tensions a l’egard de la 
resistance du material aux points en danger. 


‚Ina pasBeTBiteHnä TPyÖ0oIIPOBONOB, HMEIOINMX BUN ÖYKBEI Y, HAeTCAH MeTOA AA BbIUUCHeHUH 
HaTPY3RM, KOTOPOH TMONBEePTaloTca HamOoJee CITAÖbIe MecTa. 
1. Einführung 


Die bei Rohrleitungen (z. B. für Heißdampf, Luft, Gas oder Flüssigkeiten) auftretenden 
Verzweigungsstellen werden vielfach als „Hosenrohr‘ ausgeführt. Meist erweist sich eine ge- 


Schnitt C-D | 


Ba 


Bild1 Bild 2 


schweißte Ausführung entsprechend Bild 1 als zweckmäßig. Das hosenrohrartige Verbindungs- 
stück besteht hierbei aus zwei schalenartigen Teilen (Preßteilen), welche an den drei Kanten 
verschweißt werden. Die drei anschließenden Rohrenden sind ebenfalls angeschweißt. Diese 
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Ausführung scheint den praktischen Erfordernissen zu entsprechen, wobei jedoch bezüglich der, 
festigkeitsmäßigen Dimensionierung der Wandung eine gewisse Unsicherheit besteht. Nach- 
stehend wird zunächst die exakte Theorie für derartige Schalen unter Innendruck aufgestellt. 
Die Integration der auftretenden Differentialgleichung führt bei genauer Einhaltung der Form 
des Hosenrohres zwar auf mathematisch unüberwindliche Schwierigkeiten; es wird aber gezeigt, 
wie sich bei sinngemäßer Idealisierung praktische Formeln für die Beanspruchungen des Hosen- 
ohres an den gefährlichsten Punkten aufstellen lassen. 


2. Mathematische Grundlagen 


Entsprechend der allgemeinen Schalentheorie des Verfassers!) lauten die Gleichgewichts- 
bedingungen für den Membranspannungszustand der Schale unter Benutzung der in der all- 
gemeinen Schalentheorie angegebenen Schreibweise folgendermaßen: (MP=0, mt—0, t—=(, 
n= 0): a 
' L#,=0, Pb trp=0.% nah 


Diese Beziehungen setzen die Anwendung der sogenannten „‚kovarianten Differentiation“ voraus. 
Bei Einschaltung der Größen c‘, welche die Transformation von den kartesischen Ausgangs- 
koordinaten zu den krummlienigen Koordinaten regeln, sowie der Größe /g, welche mit der 
Determinante der c/ übereinstimmt, läßt sich Gleichung (1) auch in folgender Form schreiben, 
welche nur noch die gewöhnliche partielle Differentiation erforderlich macht: 


0 P = 
Dan LPG) + PN“ g=0 Are RE 


Zur zweckmäßigen Behandlung des vorstehenden Problems dient die Einführung solcher Schalen- 
koordinaten, daß ihre Projektion in der x, y-Ebene mit den x, y-Koordinaten zusammenfällt. 
Es wird daher gesetzt: 


er =u=al gg? rer WERTEN Pi 
Für die Transformationsgrößen ergibt sich dann folgende Tabelle: 


+ (6% 


In der dritten Zeile dieser Tabelle sind die Komponenten des Normal-Einheitsvektors 
angegeben, welche sich aus den Bedingungen der Orthogonalität zu den Tangentenvektoren c* 
und aus der für den Betrag 1 dieses Vektors erforderlichen Bedingung 


SUN ER 2A Se 
errechnen. Die Größe g ergibt sich zu 
galt a tl en re 
Die in Gleichung (2) noch auftretenden Komponenten des Krümmungstensors ergeben sich zu 
1 
Di = vg CR. De = \ö RR de = ga uk a ae re 


Die Gleichungen (2) nehmen nunmehr folgende Form an: 
0 1) e 
11 @ 12 ı/\ __ u 
„em +, rm) —pzu=0, | 


a4 MW—=0, ae 
I Ygzu +2 yo +2 Ya +pg=0 
Die ersten beiden Gleichungen lassen sich mit Hilfe des Ansatzes 
Ryg=Futpz, I®Wg=-FuHtpz, Da yg = FW NR) 
!) H. Neuber, Zeitschr. f. angew. Math. u. Mech. 29 (1949) S. 97, 
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durch eine Spannungsfunktion integrieren, was bereits durch Pucher?) angewendet wurde. 
Die dritte der Gleichungen (8) liefert folgende Differentialgleichung für die Spannungsfunktion 


FoZu —- 2Fuzut huzo PER tz) tät t1=I. . (10). 
Auf Grund der Identitäten 


9 dran = @zulu = ZZuut Zu USW. At he re) 


erhält man hieraus 


a le, uu un —2F, u r 4 a 


2 22° 
2) TEN 
+ Bud l-pl-1- 440] ..... a2. 
Mit Einführung der Größe 
a 3 De En ee le 
folgt hieraus die Differentialgleichung 
1 W’, [ WW; W’, } 
2yw [Fon Wu >. >) en —2F, un Wiuo 2 9) w| ae Wir 5, zw} =D Be zAw 


(14). 
3. Integration der Differentialgleichung 


Die Integration dieser verwickelten Differentialgleichung zweiter Ordnung mit veränder- 
lichen Koeffizienten führt im Falle des Hosenrohres auf mathematische Schwierigkeiten, so daß 
eine Lösung in geschlossener Form mit übersichtlichen Endformeln bei exakter Einhaltung der 
Randbedingungen nicht zu erwarten ist. Es soll daher zunächst diejenige Lösung der Differential- 
gleichung diskutiert werden, welche sich in geschlossener Form darstellen läßt. 


Diese Lösung entspricht dem Ansatz 


Be WER ER ER n. RIS) 
Es ergibt sich hieraus 
3 3 i 
Fu=zZ C.WiR.W, .; Fu = 5 c wır|w er USW. De, (16). 
Unter Umgehung der Zwischenrechnung geht Gleichung (14) über in 
3 Ci w’ | Wi >| West 
a Me Fo Meng 2 Wet Me 
Wu wi 1 
+ [Wu + ae |w. Sl=pl1— 5 Aw . (17). 
Die Ausrechnung liefert folgende nichtlineare Differentialgleichung zweiter Ordnung: 
2 i 
ap 
Wen Wo Wu=zc 1m, am N 


Lösungen dieser Differentialgleichung werden offenbar durch die Kegelschnitte gegeben, d.h. 
durch Flächen von der Form 


u? v2 Ge: 


Bei positivem Vorzeichen handelt es sich um ein Ellipsoid (einschließlich des Spezialfalles der 
Kugel), bei ein oder zwei negativen Vorzeichen um ein Hyperboloid. 

Wir berücksichtigen zunächst nur die für das Ellipsoid geltenden oberen Vorzeichen. 
Dann folgt: 


c? (= 
Die Differentiation liefert: 
2 2 , 
gene = — Fr R W, DO en De 3 W, uvm 0 Bee He ACER (21). 


2)A. Pucher, Beton u. Eisen 33 (1934), S. 298. 
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MD RENT Een ee Ser SE RE 


Durch Einsetzen in Gleichung (18) folgt: 
_pebie 00 
rast U Anne Ach an Ver (22). 


Weiterhin werden die Komponenten des metrischen Fundamentaltensors benötigt. Wir be- 


rechnen hierzu zunächst: 
W;u eu ev 


Ze ee a Bee en m ee (23). 
“ESyW ayW ” ByW Trrs 
Dann wird r A ar 
e u 
-1+2,-4]1+(5-1)5-5 I (4), 
B u. 
| +2, - 1-5 1 ac ee (25), 
Auv 
Ap = Zur = a2? W nr eh ea RT at ee (26). 
Ferner ergibt sich 
b c2 es u? e \ 
hir 141) Re (27). 
Die kontravarianten Komponenten des metrischen Fundamentaltensors folgen hieraus in der 
Form: 
2%. 5 % 
Ay en 
a m u (@ ) 2 2 ) Pre (28), 
ea ara! 
2 \ 172 v2 
rt G Gh 
Ay a a b? 
ER m ne a a a un ne EN 
i Ba e 1 
at yr 
Aauv 
a a? b? 
aa — U _ _____ re EnZEe 2 (30) 
x @ AaUTz za ial eu Ä 
Br ET 


Die Komponenten des Längskrafttensors ergeben sich entsprechend den Gleichungen (9), (15) 
und (16) zu: 


3 2 
LM vg - 5 C WıR2 In. 2 +p wur, ] 
= 3 2 W°’, 
IyVg= 5 GC wıR2 Im. +3 +. PWERN 20, ws (31). 
mn Fewelw..4 ww 
Wird zur Abkürzung 
@ u? & u? 
14 (51a) 6») 
gesetzt, so folgt 
ce 
y ala- 1 ia Ice Rare (33) 
W 


und die Gleichungen (31) gehen über in 


W c? cd v2 
LA zn SA in ge 
el nremtn| 
f W ce ct u? 
22 < \ 
1 be-s+an) tel. ur rt, 
12 _ 3C Auv 
\g ab 
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a = Ne ergeben sich aus den Längskräften mit Hilfe der Beziehungen (h = Wand- 
stärke): 


EN RT a en 


Die Ausrechnung liefert: 


a 
On wu Dear 2 Er = wg nr me a (36), 
2ehli tl ) [7 


| b2 c& Di e le 
ler rer 
SETRBRLIN RN) EIKE EIN ANEBE SL SUR UNE 2 
C u v j 
4a 


0, = 


Zeh 


Für die Anwendung auf das vorliegende Problem interessieren insbesondere die Spannungen 
an der Stelle u=a und v=0 (Punkt A, Bild 2). Es ergibt sich: 


a € = b 
De er 
pa Dem abe 
OgA = In 1 no FE F2 ner (39). 


Entsprechend Bild 2 betrachten wir nun ein Hyperboloid mit elliptischem Querschnitt. Die 
Achsen a und c in Richtung der u- und z-Achse bleiben erhalten, während an die Stelle der bis- 
herigen dritten Halbachse 5b (in v-Richtung) des Ellipsoids die Größe b* tritt (imaginäre Halb- 
achse). Für diese Größe gilt die Beziehung (vgl. Vorzeichenwechsel in Gleichung (19)1): 


Dame ED LE a Bere = 2 (40). 


Wird gleichzeitig noch beachtet, daß der Krümmungshalbmesser an der engsten Stelle des 
Hyperboloids im Längsschnitt am Ende der in u-Richtung liegenden Halbachse der Beziehung 
(Punkt A, Bild 2): 

b*2 


Bu 


ea 


entspricht, so folgt 
MRIISOATITERIT Immer. ee a2), 


Mit Hilfe dieser Gleichung erhält man aus den Gleichungen (38) und (39) unmittelbar die Span- 
nungen für das Hyperboloid. Es ergibt sich 


pa (a 

Eee 14 lo +1). tr), 
pal 0, (a? 

Sant Alle wat. ER, Arabien.” (Ad): 


Wird ferner noch der Krümmungsradius des engsten Querschnittes des Hyperboloides (u, z-Ebene) 
entsprechend der Beziehung 


EEE Rh Oak le ya Aär; (45) 


= a 


eingeführt, so nehmen die Beziehungen (43) und (44) folgende Form an: 


218 ZAMM 39 (1959) Heft 5/6, Seite 218—236 


Schließlich interessieren noch die Spannungswerte im Endpunkt der kleinen Halbachse c der 
Schnittellipse (an der Stellez=c, u=v=0). Aus den Gleichungen (36) und (37) folgt mit 
Der; 


FL ne ak: a 


ee N a er 


Mit Einführung der Krümmungsradien entsprechend den Gleichungen (42) und (45) ergibt sich 
mithin für die Spannungen im Punkte B 


IR a 


RR) 8 


Man erkennt, daß diese Spannungen mit den Spannungen im Punkte A im einfachen Zusammen- 
hang stehen: Es gilt: 


nn (oBl 


Hierbei wurde zur Abkürzung 


ee a a en 


gesetzt. 


Manuskripteingang: 12.7.1958 
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Einige strenge Lösungen für den Spannungszustand 
in ebenen Verbundkörpern 
Von Hans Bufler 


Es wird der Spanmungszustand in einem Verbundkörper (bestehend aus einer starren endlichen Scheibe 
und einer elastischen Halbscheibe oder aus zwei verschiedenen elastischen Halbscheiben) untersucht. Das 
Problem läßt sich auf zwei gekoppelte singuläre Integralgleichungen für die Druckspannung p und die 
Schubspannung q längs der Berührungsfläche (diese beiden Größen interessieren hauptsächlich) zurück- 
führen. Es werden die strengen Lösungen für p und q sowie für die tangentiale Randspannung oz ermittelt 
und auf fünf verschiedene Belastungsfälle (reiner zentrischer Druck bzw. Zug, Schub ohneVerdrehungs- 
möglichkeit, Schub mit Verdrehungsmöglichkeit, reine Verdrehung, Wärme- und Vorspannungen) ange- 
wandt. Die dafür erhaltenen Ergebnisse sind in dimensionsloser Form in Diagrammen niedergelegt. 


The state of stress in a compound body (existing of a rigid finite sheet and an elastic semi-infinite sheet 
or of two different elastic semi-infinite sheets) is investigated. T’he problem can be reduced to two coupled 
singular integral equations for the pressure p and the shearing q stress (these two quantities are of the main 
interest) along the contact-surface. T'he rigorous solutions for p and q and for the tangential-stress oz along 
the surface are found and applied to five different load-cases (pure central pressure resp. tension, shear 
without twisting, shear with possibility of twisting, pure twist, thermal- and initial stresses). The results 
obtained for these cases are laid down dimensionless in diagrams. 


L’&tat de tension dans un corps compound (consistant d’un disque rigide et fini et d’un demidisque 
elastique ou bien de deux differents demi-disques Elastiques) est examine. On peut reduire le probleme & deux 
equations integrales, singulaires liees ’une & l’autre, pour la tension de compression p, et la tension de cisaille- 
ment q le long de la surface de contact (ce sont ces deux quantit&s qui sont surtout interessantes ). Les solutions 
exactes pour p et q aussi bien que pour la tension marginale tangentiale 0, sont determinees et appliquees & 
cinqg differents cas de charge (pure pression centrique, respectivement traction, cisaillement sans possibilite de 
torsion, cisaillement avec possibilite de torsion, torsion pure, tensions thermiques et tensions initiales). Les 
resultats y obtenus sont donnes comme diagrammes sous une forme sans dimensions. 


n 
» 
> 

. 
. 
> 
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Hecaenyerca COCTOAHUE HANPSUKEHUA B COCTABHOM Tele (COCTOAIEM U3 FKECTKOU KOHEYHOH 
IIIACTHHKU MU YIACTUYHOH HOJIYINIACTHHKU HIN ;ke U3 MBYX Pa3uyHBbIX HITACTUYHBIX IIOJIY- 
INIACTUHOR). 3anaya CBONUTCH K IBYM CBA3AHHBIM M&KIY CO00U CHHTYJIAPHBIM UHTETPAIBHEIM 
YPABHeHUAM NIIA HalpsBKeHHA NaBlJIeHuUA pP MU HAaIPsBKeHHA CcABHUTA q BNONB HOBEPXHOCTU 
CONPMHKOCHOBEHUA (9TH ABe BEJIHYUHBI IPeNCTABIAIOT HAMOONBMUH HHTEePec). OnpeneIstotch 
CTporTue pelleHuA MIA p Mg m INIA HaupsBRKeHUA B HaupaBlleHum Kpası 0,, KOTOpbIe 3aTeM 
UPHMEeHAWTCH K IIATM CIIYyYaAM HATPy3KH (UMCTOE NMEHTPAJIBHOE JABJIEHHE COOTBETCTBEHHO 
TATA, CABUT 6e3 BO3MOFKHOCTU KPy4eHHuA, CHBUT C BO3MOFKHOCTBIO KPYYeHUA, YUHCTOE KPy4eHNHe, 
TeIlIOBble U IIPelBApHTelIbHbIe HanpsskeHun). IloıyueHHbIe pe3yJIbTaTu HaHeceHbI B ]IHA- 
TpaMmMe B Öe3pa3MepHOM Bupe. 


1. Einleitung 


— Ein Verbundkörper besteht aus zwei oder mehreren Teilkörpern aus verschiedenen Stoffen 
(d.h. mit verschiedenen elastischen und thermischen Konstanten), die fest miteinander in 
Verbindung stehen. Der Kraftfluß von einem Teilkörper in den anderen erfolgt längs der je- 
weiligen Berührungsfläche durch Normal- und Schubspannungen (p; und q,): Für die allgemeine 
Lösung des Problems ist zunächst die Kenntnis der Verformung der Teilkörper als Funktion 
der p; und g; erforderlich. Aus den durch die Haftung gegebenen Bedingungen für die Ver- 
formung an den Berührungsflächen erhält man die Integralgleichungen zur Bestimmung der 
p; und q.. 

Im allgemeinen Fall scheitert die Lösung der Aufgabe bereits an der Schwierigkeit der 
Berechnung der Verformungen. Daher wurden erst für gewisse Spezialfälle einige Lösungen ent- 
wickelt, siehe z. B. [1], [2] und die unter [2] angegebene weitereLiteratur. Die exakten Untersuchun- 
gen beziehen sich hierbei auf eine unendliche Voll- oder Halb-Scheibe, welche eine unendliche 
oder halbunendliche elastische geradlinige Versteifung besitzt und verschiedenen Beanspruchun- 
gen unterworfen wird. Die Ergebnisse werden meist in Form unendlicher Reihen angegeben. 


In der vorliegenden Arbeit werden als Teil- 
körper eine elastische Halbscheibe und eine starre 
endliche Scheibe oder 2 elastische Halbscheiben (siehe 
Bild 1und2), dieim Intervall— asx = afest mit- 
einander verbunden sind, zu Grunde gelegt. Die Auf- 
gabe führt dann auf zweigekoppeltesinguläre Integral- 
gleichungen, deren Lösungen sich streng und in ge- 
schlossener Form angeben lassen. Integralgleichungen 
desselben Typs spielen übrigens auch in der Theorie 
der rollenden Reibung [3] und für die Bodenpressung 


Bild1. Aus einer endlichen starren Scheibe und einer 


elastischen Halbscheibe bestehender Verbundkörper von Staumauern [4] eine wichtige Rolle. 
- + — o= — u 9 
X— 
05 a a) 05 
Y 
Bild2. Aus zwei elastischen Halbscheiben bestehender Bild3. Spannungen an der Halbscheibe 


Verbundkörper 


2. Ausgangsgleichungen 


Wir betrachten eine im Bereich —a=<x=.a durch den Normaldruck p und die Schub- 
spannung q beanspruchte elastische Halbscheibe, die den Elastizitätsmodul E, die Querdehn- 
zahl v und die Wärmeausdehnungszahl « besitzt. Die Halbscheibe sei ferner mit der Spannung o, 
in x-Richtung vorgespannt (siehe Bild3) und habe im Vergleich zum Ausgangszustand die 
„Übertemperatur“ t, so daß auch die Wärmedehnungen mit erfaßt werden. 

Für die Dehnung eg, längs des Randes gilt dann bei Voraussetzung eines ebenen Spannungs- 


zustandes: 
— IR t 1) 
I aa) (EBD alla: 22 et N aa ERICH: 


Die Rand-Normalspannung o, (in (1) bis (4) bezieht sich o, auf den Rand der Halb- 
scheibe) setzt sich dabei aus drei Bestandteilen zusammen, die von der Vorspannung 0, dem 
Normaldruck p und der Schubspannung g herrühren. Es ist 

ee N ci 
I6- 


F; ü u; u: U 
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d 
2 De EN (3) 


(hydrostatischer Spannungszustand, s. [5], S. 37). Da für eine mit der positiven x-Achse gleich- 
gerichtete, bei x = 0 angreifende Schubkraft Q 


2 1 
Be 95 


gilt (siehe [5], S. 29), so folgt für die kontinuierliche Schubspannungsverteilung q — indem man 
Q durch den Laststreifen q(u) du und x durch x — u ersetzt sowie über den gesamten Lastbereich 
integriert — 


A | 
LER? Mu u (4). 


7 
Das Zeichen $ bedeutet, daß das Integral wegen der Singularität des Integranden bei 
u=x für x®< a: als Cauchyscher Hauptwert zu verstehen ist. Aus (1) erhält man somit 
mit (2), (8), (4): 


+a 
0E 1 2 q(u) 
rn an ale, are (5). 
; 
Im weiteren benötigen wir die senkrechte Randverschiebung 7 bzw. die Ableitung = 
Zu diesem Zweck berechnen wir erst die Dehnung 
nr] 

a Ei 9 Verne ee (6) 


(n zählt positiv ins Innere der Halbscheibe).!) 


Der Beitrag, der vom Normaldruck p herrührt, wird wie folgt ermittelt (s. [5], S. 61): Man 
geht zweckmäßig von einer im Koordinatenursprung angreifenden Normalkraft P aus, welche 
die Normalspannungen 


2 a 
Ne N Be sl 
q; = P + ii (7) 
und 
2 n® 
er it Een eh, A ae 5 re 
O, TC 15 (23 + y2)° ( ) 


in der Halbscheibe zur Folge hat. Setzt man diese Ausdrücke in (6) ein und integriert zwischen 
den Grenzen y = 0 (obere Grenze) und y = h (untere Grenze), so erhält man für die Durch- 
senkung des Randes relativ zur Linie y = h: 


pP 2 h? 
Msn. mi & . (9). 
n 3: hazatl may (9) 
Läßt man h im Vergleich zu x sehr groß werden, so folgt: 
P zn 
Na — ; it 
n TE In (4 +1+v (10) 


Um auf eine kontinuierliche Lastverteilung überzugehen, ersetzen wir P durch den Laststreifen 
p(u) du, x durch @— u und integrieren über den gesamten Belastungsbereich: 


+a 
nd NER da 5 sa A 

Daraus folgt für den Differentialquotienten ii. s 
07) EN E Plu). em (12) 


EEE _ an 


!) Die Wärmedehnung sowie die Vorspannung brauchen hier nicht berücksichtigt zu werden, da sie 
; r on 
keinen Einfluß auf m haben. 
X 
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| Zur Ermittlung des Beitrages zu n, der von der Schubspannung q kommt, sei von einer 
im Bereich— c = x = ckonstanten, mit der positiven x-Achse gleichgerichteten Schubspannungs- 
verteilung q, ausgegangen. Diese führt zu den Normalspannungen (s. [5], S. 39) 


2 T. ıl 
= — ln 3, 9 (08 20, — cos 20) N? (13) 
und 
1 
y=—5,% (0082.07 7.003.200)... u ee (14), 
wobei 
een 2? 
n=yy +(c+2)%; cos2u =1 ArC+a‘ 
a EEE (15). 
RIM. RA 2 y: 
nr = Yı? — 3 Pe ee es 
2 Yy NER) COS 2.0, y + (c— a)? 


Nach Einsetzen von (13) und (14) in (6) unter Benützung von (15) und anschließender Integration 
über y bekommt man 


A +cte! 
„=, |( I), 


en ee ae ee SL ee a ur (16). 


IR= (—a)arctg I —(c+x)arct 2. 
e—ıx Sctz 


Erstreckt man das Integral zwischen den Grenzen y = 0 (obere Grenze) und y = oo (untere 
Grenze), so folgt unter Umgehung der Zwischenrechnung 


ea) (1—»)x fürn <’c, 
A E 


eg nr ne (17) 


und nach Differentiation: 


On‘ go F= 
EB SE ER Fe (18). 
—=0 fürie| > ce 
Bei einer veränderlichen Lastverteilung qg kann man eine Aufteilung der Belastung q = q(«) in 
einzelne Laststreifen q(x) dx vornehmen. Für jeden einzelnen Streifen gilt dann die Beziehung 


18) (nämlich, daß z roportional der Belastungsordinate an der betrachteten Stelle ist und 
a 8 


unabhängig von den Belastungen außerhalb davon), so daß man allgemein schreiben kann: 
On'® (1 — v) 
= —. Es N a WR ee 019) 
97 FE (0 (19) 


Gleichung (19) wurde von G. Heinrich auf andere Art abgeleitet (siehe [6], S. 367); allerdings 
enthält seine Formel einen Vorzeichenfehler. 


Für den in Bild 3 betrachteten Belastungsfall gilt somit wegen (12) und (19) nach Über- 
lagerung der von p und g herrührenden Anteile: 


+a - 
nnd p(u) iu ELTERN 0 0 3 


t—u 
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3. Aufstellung der Integralgleichungen 
a) Elastische Halbscheibe — starre endliche Scheibe 
Das hier wirkende Kräftespiel ist aus Bild 4 zu ersehen. Für die Halbscheibe (Index 2) 


gilt gemäß (5) und (20) | 


+a 
0 1 2 u 
ae >E mm Pd —_ . gen a Be: ET: Sa 
+a 
Mm _ Mm_Ilıa 24 rw 99 
m - le v,) 4(%) Eu a 


Bild 4. Spannungen an einem aus einer endlichen starren Scheibe und einer elastischen Halbscheibe bestehenden Verbundkörper 


und für die starre endliche Scheibe (Index 1), die dieselbe Temperatur besitze, 


088 | 

(&,)y=0 = dr, = — = u Uri uch: Beh ae See 
a ip 

a“ ze DEM. a wen se a 


Infolge der Haftung gelten längs der Berührungsfläche die Bedingungen (e,,),-o = (&.)y=0 


(% x \ : i 

und = = 9 (wobei p den Verdrehungswinkel der starren Scheibe relativ zur Halbscheibe in dem 
2 

in Bild 4 angegebenen Sinn bedeutet).2) Dafür kann man auch schreiben: 


0 

„tt BE ie ar (AN 
und 

0) 2 6} 

an Na == © x = EEE 5 ‘ . . (26). 


Unter Verwendung von (21) bis (24) gehen die Bedingungsgleichungen (25) und (26) über in 


-+a 


. q(u) 7 7 Fee 

q RT du + 5 1—3): pe) = 5 9 a 
—ı 

+4 ( ) 

p(u 7 Te ? 
q RER Ar 9 (1 — 9) ga) = — 9E2 p a E. 
wobei 
= ya la) En niart + JE ao arme Da 


b) Elastische Halbscheibe — Elastische Halbscheibe 


Für das Kräftespiel ist Bild 5 maßgebend. Es wird vorausgesetzt, daß nur im Bereich 
Ss54 Berührung stattfindet. Damit das auch im Fall einer Verdrehung oder Zusammen- 
drückung der beiden Scheiben möglich ist, muß die Begrenzungskurve mindestens einer Halb- 

?) Der ebene Spannungszustand erfährt in unmittelbarer Umgebung der Haftzone wegen der Behin- 


derung der Querdehnung (senkrecht zur Scheibe) eine Störung, ‚die jedoch bei ü j Io 
f vi ’ e (0) ” 
hier voraussetzen wollen, unbedeutend ist. x J i sehr dünnen Scheiben, die wir 


e 


a zn 
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scheibe etwa so verlaufen, wie sie in Bild 5 gestrichelt angedeutet ist. Es gelten für die beiden 
Halbscheiben die Gleichungen (5) und (20) (der Index 1 bezieht sich wieder auf die obere, der 
Index 2 auf die untere Halbscheibe). Dabei ist lediglich zu beachten, daß für die untere Scheibe 


%=T, Wu, =: P=Pp 
und für die obere Scheibe 
4=%, eu g4=eh—ph=Pp 
zu setzen ist. Dann folgt 
+a 
0E 1 2 u 
(&,)y=0 = = = E, In —Aen) DR) = £ ut. m TMml..... (30), 
—a 
PR ER ae 
Na __ Ma _ ”E a p(u : 
I ß v,) (X) „$ ERaRT m a TE (31), 
dE | ae (u) 
RITTER L a 
(&,)y=0 z Ox, == dr 73 I (di v,) P(&) Si zz zu m +at ie (32), 
ö et 2 Kia) 
A NETTE N 
Fr Erg ß )q&) + E $ Be au ERRTZEEEENEE (33). 
Die durch die Haftung bedingten Forderungen lauten (e,,),-0= (&,)y-o und + 2 = n ; 
oder in anderer Form geschrieben: 7 7; 
0 
BE (&, +$&)=0 DEI be ER ETE (34), 
0 2 2 
a mtmW=P a EBENE (35). 


Bild5. Spannungen an einem aus zwei elastischen Halbscheiben bestehenden Verbundkörper 


Benützt man (30) bis (33), dann erhält man aus (34) und (35) folgende Integralgleichungen für 
p und g: 


+4 
au) auUmwW)HhTlgsmWE, ll EE | Zune Cr 
+7 Eos A a 
Bea. (00) 
und 
+a 
Beeren, ren. base han (1 < ugi 
en 1 EHE, =) ART Eu 2. 


In (386) und (37) sind übrigens auch die Lösungen (27) und (28) als Spezialfall enthalten. Man 
erhält sie, indem man E, = ® setzt. 
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4. Allgemeine Lösung der Integralgleichungen 
Wir schreiben die oben abgeleiteten Integralgleichungen in der Form: 


dr 2 u. K)-G| "7 Aa ne (38) 
er oe (39), 
wobei im Fall zweier elastischer Halbscheiben 
a ar ER 27T 
G=-5 Eo Pe (40) 
en an 


gilt. Liegt der Fall „elastische Halbscheibe — starre endliche Scheibe‘ vor, so braucht man 
in (40) lediglich E, = w zu setzen. 

Das Integralgleichungssystem (38), (39) kann in eine Integralgleichung 2. Art umgeformt 
werden, welche sich mit Hilfe der Neumannschen Reihe näherungsweise lösen läßt (s. [7], 
S. 427). Zur strengen Lösung führen wir die beiden Integralgleichungen (38) und (39) in eine 
komplexe Integralgleichung für die gesuchte Funktion 


f{@) = p(x) + iq(@) Peir ar = Be bl > (41) 


über. Zu diesem Zweck addieren wir Gleichung (38) und die mit i multiplizierte Gleichung (39) 
und multiplizieren die so gewonnene Gleichung mit i. Sodann erhält man: 


+a 
em wre es 
Nach Einführung der dimensionslosen Koordinaten 
i=-T; ü=. 4 
er BE nn Be el, als An 
wird 
f@) = p@) +iga@) = FÜ) =PÜM)+iII@)......2.... (4) 
und obige Gleichung geht über in 
+1 
« i F(u) C ‚c 
F(&) — TR: _ ee rar Ma. ao 
(X) or Ta—ı du ro cc Ri 1 EEE. ° * 
Diese Integralgleichung ist vom Typ 
+ 
RE a 
F@) +74. dü=6@) <1.... (6) 


— 
der von H. Söhngen [8] ausführlich behandelt wurde. 
Durch Koeffizientenvergleich bekommt man: 


= —iR (wobei 7 > 2wegen (0)... 2.2... . (6) 
En 0 “ol 
ee 
a en wine Re 


Nach Satz 7 der oben zitierten Arbeit gilt: Ist G@) in (—1, 1) eine Funktion der Klasse 
L? (p > 1) und werden solche Lösungen von (45) gesucht, die in (—1, 1) einer Klasse La > 
angehören, so gilt, wenn für} # +i 
1— iA N re 
rm mit — oRl u Tee Be 
ist (Re bedeutet Realteil): 
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Ist A=i4A, mit —1<A, <s1, so hat die homogene Gleichung nur die triviale Lösung. Liegt 
, Bar in diesem Bereich, so wird die allgemeine Lösung der homogenen Gleichung gegeben 
urc 


N a + z\% 
Keen; —) NETT. (49), 
wobei 
+1 
Ges) Polc)de 
| 
ist und der Ausdruck \ 
GR A 1 ea a Gu) 
F LIE : 
Dr 1+% 1+% nr 1—ı\l+u 1—r ı— ea 0 20) 
liefert die allgemeine Lösung der inhomogenen Gleichung. Die Konstante 
a 
GET, Player a a er (51) 


kann dabei beliebig vorgegeben werden. 
> 2 ist, die Lösung (50) zuständig. 


In unserem Fall ist also, da | “= 


Für & erhalten wir wegen (48) und (46): 


ur EIER Me ee (52) 
mit 
a ee 2 (53). 
| 
Co 


Die Konstante C errechnet sich nach (51) und (43) zu 
+1 h +a 1 
Gi | r@«@ - r | ine) + ia@)ı a ee) a (54), 
= un) 


wobei P die Normal- und Q die Schubkraft pro Scheibendicke ist, die längs der Berührungsfläche 
übertragen wird (P wird positiv als Druckkraft gezählt). 
In (560) wird — da G(u) = const ist — folgendes Integral a 
+1 
er, «- (17 en du 
irn z—u 1) dü+(—3) Fun: 


A pa Ne N 
Das erste Integral geht mit der Substitution u = = über in 


e2 +1 
2 
a u\® -2[ ee) 
" 1+u (e@ + ja % 
Zur Lösung benützen wir die in [8] Gleichung (19) und (21) angegebenen Integralformeln 
+® 
1 s2 —1 2 
er al) 0<Res<n.... (5) 
PLA 2 (1 7 raw =. ) — Vsinzs n ganzzahlig, positiv ( 
und 
pn 1—u\ du 1—ır i 
1 —u\ ey 
EB ae, | Zr an R ze 
= a na Ban el - (66) 
—ı 
Damit erhält man: 
+1 
1—ue1—u,- al — (+20) 57 
er an | ar ze a 
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Die allgemeine Lösung (50) lautet nun mit unseren Bezeichnungen unter Berücksichtigung von 
(57), (46), (47), (82), (53), N eh 


zirR 
ak + il, ) 
= —ı —— || — 26 — GE +il= +2 — G8]|. . (58). 

Trennt man (58) in Real- und Imaginärteil, so bekommt man wegen (43) nach einigen 
Umformungen die allgemeine Lösung für die Druckspannung p und die Schubspannung q als 
Funktion von t: 


wobei 
P 3 Q 
A= 72GP—-G8; BT 


VICBue, Br = Bun (ı +2) (61). 


Der Einfachheit halber ist hier — wie auch im folgenden — der Querstrich bei p und q wegge- 
lassen. 


Es soll nun noch die Relation zwischen dem Moment M (je Scheibendicke) 
+a +a 
M=f2-Dea)yds a u STB Fa Fra (62) 
—a —a 


und dem Verdrehungswinkel » bzw. der Konstanten C, abgeleitet werden. Wir gehen dazu von 


(58) aus und bilden 
+1 


h +1 E „= 2)" 
S 2 1 72 .[Q 1—V _ 
z. Piz) de = = —2C +2 +2 )- —— te 
J Ba a2 rilz r2aß Te 
+1 6 Zeile 
1 1-2 _ 
u Se FT G-+tiG me er a ee EN 63). 
ze (63) 
Mit der Substitution 
en el 
"Oer+l 
erhält man 
nr ar ee Pe 
Fa R ei (& +0) ie (2 +) 
—rd= ek dz — Be, 
a at Y er „ern 
und 
+L/1 + 3% +8 A +® E +% j 
=). (ee ee), [ehr 
—T!d= — dz—2 —_ dz + : Ba: 
J yı—22 KOITERT), KA )\, 5 e@+D% 
Zur Lösung vorstehender Integrale verwenden wir die Formel (55). Dann ergibt sich 
Pl -F ar 
Ir. 3, Br 
— tt Malen einen 64 
| ia WET, N 
und 
di + | iß 
1—:/ _ n (1—4) 
a re iS 4 
Hi =» x ee (65) 


= 
Somit folgt mit Berücksichtigung von (53): 
+1 


EF@) di = — 


—1 


er ch 
arm +28] -ilSa+am—ae2]|.... (o 
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der Vergleich der Realteile von (66) liefert wegen (62) die gesuchte Beziehung zwischen 
und C;: - 


5. Ermittlung der tangentialen Randspannung o, 


Die längs des Randes verlaufende Spannung o, beträgt nach (2), (3), (4) (für die untere 
Halbscheibe) 


+1 
na TiDeka Tr ra . (68) 
Zt 
bzw. 
+1 Ä 
= 0: | en Tre I. . (69), 


denn bei 2? > 1 ist wegen der freien Oberfläche p = 0. Zur Berechnung der vorstehenden Inte- 
grale greifen wir auf die komplexe Lösung (58) zurück und bilden zunächst 


4 4 ( er ar dü + (i 2% ya du 
PIUN- = 1 1—u 1—u 
= om du en D; a = Fe 3 —D, L —_ = (70). 
Ben Ym® — Q} 4 Y1—-m@w—i) 4 Y1—-% (ü—2) 
Hierbei wurde zur Abkürzung 
pP ß ; 
= 26) +2 +20), De tserti@r rest (71) 
gesetzt. Durch Umformen von (70) bekommt man 
+1 
Fe... 1 = 
9 Fra N ar (72) 
mit ; ER 
ee Aa 
1—u 1—u 
Ir q — J= ——— 
| yı — m (wmw— 2) | yı—ız 
Nimmt man die Substitution u — — vor, so läßt sich auf das Integral 
ER (5 +8) 
Js a e+1 f 
die Formel (55) anwenden und es folgt 
ER Zul. . (73). 


jap 
Bei der Berechnung des Integrales J, hat man zwischen den Fällen 2? < 1 und 7? > 1 zu unter- 
scheiden (nur für ?? <1 ist (72) ein uneigentliches Integral). Im Fall x? < 1 setzt man 


= nn ; Fe 4 Ereliinr th) vieler: (74) 
und erhält 
a Pe i 
l+e [ee (a ) d 
LETTER ee 


—o 


Zur Lösung verwenden wir die in [8], Gleichung (20) angegebene Integralformel 


1 7 ine 
8. de= e7’tet0ns WEN hes<t. 
n e— e 


>= 


Sie liefert das Ergebnis 


J; =in- Tan - er &o) 2 Sr EN. 


td |m 
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Im Fall > 1 schreiben wir 


1 ZEN > 
u fe = 1 . . . . “ z . 76 
rg And Fun (£ reell für 82 > 1) (76) 


H | 
we et") 
Be ir 


Dieses Integral kann mit Hilfe der Residuen-Methode gelöst werden. Wir betrachten dazu 
vorübergehend z als komplexe Veränderliche. Der Integrand f(z) besitzt dann in der oberen 
Halbebene bei 


U= 


und bekommen 


z=£E+in(l+2n); = 0,152: 
Pole 1. Ordnung und es folgt nach dem Residuensatz 
ur: En, 3 
(-3 +iB) &+i2) 5 e: +0) eu 
n=0 
wobei das Integral sich auf die gesamte obere Halbebene erstreckt. Durch Aufteilung des Inte- 
grationsweges erhält man 


Öle) de = ] K)de + S HREn) Rierdp. 


(reelle Achse) Ro» 
(Halbkreisbogen) 


fe) de=2ni:e 


Wie man leicht nachweisen kann, verschwindet der Halbkreisbogenanteil und es wird daher 
nach einigen Umformungen 


+® 
n 
ei denn rer a 
gi I() dz 076) z TIET; e e 
Das Ergebnis für J, lautet somit 
= gr nz et! I BT ee, nr So 


Setzt man nun (73) und (75) bzw. (77) sowie (71) in (72) ein und beachtet, daß wegen (53) und 
(74) bzw. (76) 


Be a er Torß=-—; 
E 1 ” 
$ a : 
En] 9 Nm für @<1 
2 
© = — 
US m; für ?>1 


+1 
Pin) er G A cos ß* — B sin ß* 
I ____@  Acosß*— Bsinß' N: 
Z =_54 MG HB +6 dr Ei. 
bzw. 


+ 
Fine x Asin ß** + B cos ß** 
I 3 ur 
| De ei +G für a>ir 2. 2 
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Wir bekommen nun aus (68) und (78) bzw. (69) und (79) unter Berücksichtigung von (59) für 
die untere Halbscheibe 


= G+r@-(2 0-1) fürs Eee a 


und 


für B>1.... (8). Br 
= |le ye—1 k 


2 A sin B** + B cos P** 
0; = Ogg — | 3 ai = p +6 
In entsprechender Weise gilt für die obere Halbscheibe (hier ist end I, au 2 du 
mit umgekehrtem Vorzeichen zu verwenden; vergleiche (30) und (32)) 


2 2 
Ga ten G-rOl6+1) 


für ®<1....(8) 


und 


T 1 ** *r* 
€ ‚A sin ß**+Bcosß +Gl 
x| y®—1 


für B>1... (84). 


2 
+ 


Es sei hier noch bemerkt, daß man zu den For- 
meln (81), (83), die für 22 < 1 gelten, nn sofort 


dadurch kommt, daß man aus nd ae 


entnimmt und in (68) einsetzt. 


6b. Spezialfälle 
Wir wenden die allgemeinen Lösungen (59) 
und (60), sowie (81) und (82) (bei der Span- 
nung o, beschränken wir uns auf die untere 
Halbscheibe) auf fünf verschiedene Belastungs- 


r, ” 
4 
E 

* 


Mr fälle an: 

£ Fall I: Es wirkt nur die Normalkraft P 

fi (s. Bild 6a). 

b Fall 11: Die Beanspruchung erfolgt durch die 

4 Querkraft Q, wobei keine Verdrehung 

- auftreten soll (p = 0). Hier ist also 

R nach (67) das Moment M= — 2ßQa yJV 

R erforderlich (s. Bild 6b). 

i Fall III: Es ist nur die Querkraft Q vorhanden 0%, 2 
3 (s. Bild 6c). e) 
3 Fall IV: Es wirkt nur das Moment M 6. % 05, 
= Bild 6d). 

r Fall V: Die Beanspruchung erfolgt durch Wär- 

FF me- und Vorspannungen (8. Bild 6e). Bild 6. Die betrachteten Spezialfälle 

{ I. Zentrischer Druck bzw. Zug 


rn . n II. Schub ohne Verdrehungsmöglichkeit 
Tabelle 1 gibt eine Zusammenstellung der Tr. Schub mit Yerärehungsmöglichkeit 
. . „ “ IV. ehun 
in den einzelnen Fällen maßgebenden Größen. len und Vorsnäfinungen 
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& Tabelle 1 
I u Im | IV V 
JRR VEREIN. VERBR 2 —  — — _—_ -  — _ Ren 
| a .0|_ 2 ,M| 
at 0. |TTraR a | Ip ‘ 
ER f) | N) | CO, nach (40) 
—— en en a un Von ln 
IR 48-0 
we 0. rap IRRTB I | aan 
ER. ei. — IM 6 RE DER 
| | Q 1-48 Q 4 M E 
3 0 a 1+4ß a 1+4P a | 0% 


gentiale Randspannung o, in der unteren Halbscheibe folgende Formeln: 
Fall I: 


N P 1 [® li 


q am—ca yı —z |sin p* 
9 x N. 
er = SZ ERS RL. I WE ee (85). 
Be = 
2P 5% Br 
TEaaı sin ß e>1 
Fall II: 
= RB RN a 
d) aym—G yı— cos ß* 
2 
ER 2 
+6 Co ) P@) = ie 00 ee (86). 
z 
2 el —.cosß** M>1 
n.a yv®—1 
Fall Ill: 
1 2 2 
A OM 1 +4 Ko, 
a) aym—c yı—a® (1 —4P®) cosß*+4Pßxsin B* 
Au = 2 
fr E G— ı) ee) | (87). 
Are er 
2 Or ee Kl n Er a S 
TEN fe [1 — 42) cosß** +4ßzsinß**) ®>1 
Fall IV: 
1 : 
Pie MT Br +2ßsin ß* 
as ay—C y1ı —= (—2ßcosß* +% sin ß* 
Br. nr 
£ E G— ) ° P(x) a<1 (88). 


= 2 
4 M 1 | rt: 
ae Tran 2>1 


22. se 


u ze A Be u 
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Fall V n- 

’\- CH 2 ß cos ß* + Xsin ß* 

q Yr® — C} Ge ji» 1 —2 — [2 [x cos $* + 2 ß sin ß*] 

0a + p@) a -. si B 
5 = x h 
+: re Besetaeen a | 
wobei E 
Pr; ent) 
und daher 
1 
net 


Die Konstanten C, C, und ß sind dabei nach (40) und Es zu berechnen. 


Da wegen (40) U <c„<, ist der für 8 mögliche Beeich 
5 a eg 


Bestehen beide Scheiben aus demselben Material, so gilt 
G=0 und daher P=0. 


Ist dagegen eine der beiden Scheiben starr und besitzt die andere die Querdehnzahl » = 0, 
dann hat man 


G=+7 nn, 


27 


Für 8 = 0 vereinfachen sich die Formeln (85) bis (89) wesentlich. Man erhält für 
Fall I: 


SE N en (90) 
a 1 TE 
0% 
0 er 
Fall II und Fall Ill: 
Plauen d2run Ale 0 N 
q an ı—® |! 
0 \ 
= | a re RID: 
05, — = 
229. Kl | 
na y®a—1 ) 
Fall IV: 2 
= an T ie 
er ıı _= 10 
: y1 Bulseneltak ati RR) An she aut ern (92). 
— pa) ı@<I1i 
0% | 0 2>1 
Fall V: = 
P_G. &E [0 
q an yı—2 nn 
(or Se ee ee (93). 
Or e er 
z 2 Ei a 
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Der Grund für die Einfachheit dieser Formeln liegt darin, daß die Integralgleichungen (38) und 
(39) für C, = 0 (d.h. $ = 0) nicht mehr gekoppelt sind und jede Gleichung für sich gelöst wer- 
den kann. Man erhält nämlich auch direkt obige Lösungen aus den beiden linearen Integral- 
gleichungen 1. Art 


a 


+4 ” 
£ an du = u und d Au) du = [9A ’ 


—Uu LU 


(welche auch in der Tragflügeltheorie eine wichtige Rolle spielen) mittels der Umkehrformel 
(siehe [7], S. 55) 


+a 
C* e "Ya 
Pu) m 
Y@a—ı mya—ır J u 
bzw. 
+a 
gl C, Ya — ı? 


q(u) Pi rl: — ® nr 
1 
ya? —  m2 Va: — u2 m. 


(Die Ausrechnung des Integrals ergibt 


Y yai a x? 
d —— dk=nzufür ?<a;?°). C* C* und C, 
u— x 
sind durch die Bedingungen z 


+4 +a 

Perl ner (daraus“C* =’P); Q= f qdr (daraus C**= 0) 
+a 9M 

M= | x p dx (daraus = -. 


festgelegt. 
Gleichung (90) ist übrigens auch als Spezialfall in der Formel für die Spannungsverteilung 


über den engsten Querschnitt in einer unendlichen Scheibe mit beiderseitiger hyperbolischer 
Außenkerbe bei Zugbelastung enthalten (siehe [10], S. 39). Es folgt nämlich aus 
a a 


0, = 3 -Coju-cosn.[2 + m 


h2 
mit 
p sin dy 


e > . h = y©Sin?® u + cos?» 
2a vy+t sind, cosp,. 


=tg?v,; A= 


a |e 


(u, v bedeuten hier Hyperbelkoordinaten) für o = 0 (o Krümmungsradius der Kerbe) und u = 0 
(engster. Querschnitt) wegen cos » — Y1 — x? sofort 
;“ ei l 
ie un. i—a' 

(In analoger Weise lassen sich auch die Formeln (91) [Schub] und (92) [Biegung] herleiten.) 

Dieselbe Druckverteilung bekommt man für einen starren Stempel, der auf eine elastische 
Halbscheibe gedrückt wird, wenn man die Schubspannungen längs der Berührungsfläche ver- 
nachlässigt (siehe [5], S. 63, oder [11]). Wegen g=0 und C,—=0 (als Folge von M—=0—0) 

a 
ist nämlich hier die aus (38) folgende Integralgleichung a du = 0 maßgebend, deren 
J z—u 

Lösung sofort zur genannten Druckverteilung führt. a 

Wie aus den allgemeinen Lösungen hervorgeht, besitzen derDruck p und dieSchubspannunggq 
sowie die tangentiale Randspannung o, an den Eckpunkten ? — + 1 wegen des Ausdrucks 
yl—a? bzw. Ye—1 im Nenner Singularitäten. (Der Zähler nimmt für = +1,d.h. 
P*, B** = + 0, infolge der Unbestimmtheit der trigonometrischen Ausdrücke sin © und cos © 
einen unbestimmten Wert an). Bekanntlich bringt aber eine — wenn auch noch so geringe, in 
Wirklichkeit jedoch stets vorhandene — Abrundung der Ecken die Singularitäten zum Ver- 
schwinden (siehe [10] und [11]) und für sehr kleine Krümmungsradien (o— 0) spielen die Ge- 
fügeeigenschaften des Materials eine Rolle (siehe [10], S. 164). Ferner wird im Bereich höherer 


3) Zur Ableitung siehe [9], Gl. 5a. 
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Bild 7. Verlauf von »,q und 0, für Falll 


233 


Bild 8. Verlauf von 9,q und 0, für Fall II 


DR 


Bild 9. Verlauf von p,q und o, für Fall IIl 


u 
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0,.497 


nach6L (a5) 5, 
— — —— nach OKUBO 
a gemäß (12) 


Bild 12. Vergleich der Näherungslösung von Ökubo mit der strengen Lösung (Druck eines starren 
Stempels auf eine elastische Halbscheibe der Querdehnzahl » = 0,257 bei vollkommener Haftung) 
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Beanspruchung die lineare Spannungs-Dehnungs-Beziehung ungültig. Aus diesen Gründen 
können unendlich große Spannungswerte in der Praxis nicht auftreten; allerdings lassen die auf 
Grund der linearen Elastizitätstheorie gewonnenen strengen Lösungen auf erhebliche Spannungs- 
spitzen in den Ecken schließen. P 


. _ Bemerkenswerterweise wechseln im allgemeinen sowohl p als auch q und o, Kurz vor Er- 
reichen der Eckpunkte unendlich oft das Vorzeichen. Im Fall I z.B. liegt die erste Nullstelle 
des Zählers von p bei 


a er 
I = E13 == | n für Co 3 ° 
+1 (s1-2.0 8#=1—- 25.104 5 
und die erste Nullstelle des Zählers von g bei 
En 
= Dan . he a > 
ee +1 a ae en na > 


Sie liegen also praktisch bei |x| = 1. Der oszillierende Verlauf der Spannungen in unmittelbarer 
Nähe der Ecken läßt sich physikalisch folgendermaßen erklären: Wird beispielsweise eine starre 
Scheibe auf eine elastische Halbscheibe gedrückt und wird Gleiten verhindert, so ist dazu eine 
Schubspannung gq erforderlich, die für £ > 0 positiv, für © < 0 negativ verläuft. Nach (19) ver- 
wölbt sich dann der Rand der Halbscheibe infolge von g nach außen, und zwar um so mehr, je 
größer q wird. Da jedoch eine Verwölbung durch die Starrheit des Stempels verhindert wird, 
ist durchaus eine Stelle denkbar, an der auf diese Weise p = 0 wird und dann das Vorzeichen 
wechselt, was wiederum die Schubspannung g veranlaßt kleiner zu werden und ebenfalls das 
Vorzeichen zu ändern usw. 


In den Bildern 7 bis 11 ist der nach (85) bis (89) berechnete Verlauf des Drucks p, 
der Schubspannung q und der tangentialen Randspannung o, in der unteren Halbscheibe 
als Funktion von X für die Beanspruchungsarten I bis V aufgetragen, und zwar für die 


beiden Extremfälle C,=0 und GC, = 5: Die Ordinaten wurden dadurch dimensionslos ge- 


j P j 
macht, daß die Spannungen p, q und o, bei FallI mit dem mittleren Druck p„ = rk bei 


Fall II und III mit der mittleren Schubspannung q,„ = = bei FallIV mit der Spannung 


Di = (dieser Wert gilt bei linearer Spannungsverteilung an der Stelle = 1), und bei 
Fall V mit der Vergleichsspannung o, 2 C, dividiert wurden. Man erkennt das schnelle An- 
IT 


wachsen der Spannungen gegen die Ecken. 


Bild 12 zeigt den nach (85) ermittelten Verlauf von p und g, der für einen starren, auf eine 
elastische Halbscheibe der Querdehnzahl » = 0,257 drückenden Stempel bei vollkommener Haf- 


tung gilt (hier ist $ = 0,124). Zum Vergleich sind die für dieses Beispiel von Okubo [12] 
durch Überlagerung verschiedener Spannungszustände in Form unendlicher Reihen gewonnenen 
Ergebnisse eingetragen, nebst den von ihm später verbesserten Werten [13]. 

Die Behandlung des Problems eines mit einer elastischen Halbscheibe verbundenen 
starren Stempels bei den Beanspruchungsarten I und IV findet man auch in [14]. Die 
dort nach ganz anderen Methoden hergeleiteten Lösungen für p und q stimmen mit den 
Formeln (85) und (88) überein. 
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Fehlerabschätzung für das Ritz-Galerkinsche Verfahren 


bei Eigenwertproblemen 
Von G. Bertram 


Auswertung und Verbesserung der im I. Teil der Arbeit!) gewonnenen unteren Eigenwertschranken 
beim Ritz-Galerkinschen Verfahren für spezielle Klassen gewöhnlicher Eigenwertprobleme. Numerische 
Beispiele. 

Application and improvement of the lower bounds for eigenvalue approximations, got by the Ritz- 
Galerkin method in the I. part of the paper !) for special classes of ordinary eigenvalue problems. Numerical 
examyples. 


Application et amelioration des limites inferieures pour les approximations des valeurs propres obtenues 
par la methode Ritz-Galerkin dans la I. partie de cet expose!) pour quelques classes speciales de problömes 
"ordinaires de valeurs propres. Exemples nume£riques. 


‚Haiotca 06pa60TKa un yılyyeHue NOJIyYeHHbIX B UHePBOli yacrı paöoTkI!) HE5KHHX TPaHHL 
COÖCTBEHHbIX 3 HayeHHMü B MeTone Puua-l'arepkuHa IA CHeIHaIbHEIX KJIACCOB OÖBIKHOBEHHEIX 
3anay 0 COÖCTBEHHBIX 3HaYeHUAX. IIlpuBonAtca YHucHeHHEIe IIPHMep&I. 


II. Teil: Spezielle Problemklassen und Beispiele 


8. Eine Klasse von Problemen 2. Ordnung 
Das bei vielen technischen Fragestellungen auftretende Problem für z(f) 
—kdz—-$2=I,dz; d=z(b)=0; ah >0, H>0inasit<sph 
(s. etwa [2], S. 410ff.) läßt sich zurückführen auf (s. etwa [6], S. 116 und [2], S. 105£f.) 
Ya) = Ag) y@; y(0) = yl) = 0. 
Die orthonormierten Ansatzfunktionen gewinnen wir aus 
—y'’=ıy, y(0) = y(l),=0''zu 
9 V2sinonx bei 0, = o!n® Belag 


Mit dem Ansatz (der Normierungsfaktor /2 kann bei der Ermittlung der Näherungen, nicht aber 
bei der Fehlerabschätzung unterdrückt werden! Die a, sind noch von v» und r abhängig!). 


r 
Dur: 2% sinonz (daran, een ra use (a 
g= 


gehen wir in die Galerkinschen Gleichungen 
1 
S met +Agd,nde=0 | 
ein. Aus der Säkulargleichung 
1 1 
Det | 9am de tR [gap pda] = 0 (a, nal 20) 


denken wir uns die den für v=|],...,r und aus den Galerkinschen Gleichungen die a, für 
e=1,...,r jeweils ermittelt. Ein Faktor bleibt frei für die Normierung. i 


!) Siehe Z. angew. Math. Mech. 37 (1957), S. 191—201. Dort auch Literatur. 


Pe U 
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Die Fehlerformel (5.7) lautet dann zunächst 


N ee, (8.2) 
1 1 A, u) 922 2 
’ 2 e 
E# u 8 m 
mit 
1 ax . - 
Hi =/f 2) y2* sin on Bela) dr nr (8.3). 
In folgenden kennzeichnen wir Maximalbeträge durch „Dächer“, z.B. |q - q’Imar = 4:4 
1. Fall: g(&) > 0 und stetig m 0 sr <sI1. 
Analog zu (5.5) und (5.8) erhalten wir dann 
d 1 
nesaesworedr Haute (8.4) 
0 0 
und 
Are Qi; 
As, ee A, = 2 rT TR I A eo Fer lie (8.5) 
mit 
bei 
So 1 1 101 
A ABB 
werz er 
x Tabelle 1 
£ r 105 - Ar | 106. Br | 106-0, | 106 -D, | 102: 9, | 10:2, 
1 65 345,48 28311 | 118 22 1,026 598 1,013 211 8 
2 40 015,18 74 29 5 1,090 761 1,266 514 8 
3 28 757,27 297 11 2 1,103 435 1,379 093 9 
4 22 424,70 14 5 ıl 1,107 445 1,442 419 6 
5 18 371,85 Se 3 0 1,109 087 1,482 948 1 
6°| 15473,88 5 2 0 1,109 879 | 1,511 927 8 
a TE : Sr ae RER j 
oo 0 0 | 0 0 1,1 = 1,6 =5 
- (Die B, usw. werden später erklärt.) 


2. Fall: g(x) > 0 und stetig differenzierbar nOsSr =sI1. 


I, wird zunächst partiell integriert; dabei verschwindet der ausintegrierte Bestandteil 
wegen y,(0,&) = y,(1,&) = 0: 


I. = A, Ka era dr Ar: „1:4 48.0) 
0 


Die Schwarzsche Ungleichung liefert 
1 1 1 1 
2 PR ee 
ee m | [g(&) y,&, ©)]'? dx | tor on rdr = Eye j [g(&) y,(z, 9)? dx. 
Ö Ö ö 


Wir setzen y,(x, &) vermöge der Integralgleichung (3.1) mit 
Liz, &, 6) = K,(a, &) + &{K(&, 9 — K,@, 8}, 


4 x(1— 8) für DEE -B. SinENT-sinonE 
Kad= una) für ı<ı=22 ee? 


Bas. S merr: nem & 


2 r2 
e=1 % 1 
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und bekommen 
12 


1 
R< A i j [9% 1,6) | La, 8,8) y,& 8) dE| dr. 


Mit (5.9) erhalten wir bei Benutzung der Schwarzschen Ungleichung, der Orthonormierungs- 
relation und des verallgemeinerten Mittelwertsatzes der Integralrechnung zunächst die Ab- 


schätzung: 
2 


1 1 
R< “ | j (L- KO yu& de) dr 


0 


za I (L» Q)'? g(£) dE | ge) u(& Pr d£ | dx 


1 


<” 2 1! [Le &, &) ga)? d£ dr. 


Differenzieren wir aus, quadrieren wir und setzen wir L(x, £, &) ein, dann erhalten wir nach kurzer 


Umrechnung mit Rücksicht auf &d —&) S — bei Anwendung der Dreiecksungleichung eine 


Abschätzung der Form 4 
yH a 10 
Rn BT Ta erde oe Sr (8.7) 
o=1 
Die J; sind darin die folgenden mit (vgl. Tab. 1) 


at 


En: : sinon&-sinonE 
4 2 = == 
ı) sinonx Inn ananeosome) dee a 


51 1 
je9e#=n: HEART a 
1. || A gi 7A’ || Ra ddr), a 
= aut 


leicht abschätzbaren a 


ı 4 ht 
v,9 e 1 
; > Pdide<q:Z 6 a K': Kr: g’ qd£dx SL ayz: P, 
00 
103 213 
nk K' a en ol [een 1,4 
u ae, IK Pd sz eV — 7-18 K':K-d’gd£d| s2g ag /—_.— 
RR S 6 dt 6 9 


5) Don 
4 Be 
3 Air q?* ddl T 8 [re rraa sg':.P, 
0.0 0.6 
ibagil heat 
er: 2 S A 
4 E23 Kr'g’’qd£dx s27 4: Vr-V, ) | fe «K .g’2 .dede on Yrr 1 
aa | 2 90 
Ö 
i 11 1 
vere BR FA | ER Sl 
ö ae 0 [mer sr 
00 


00 


(8.9). 
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Damit erhalten wir 
= ln (z AB ; Yz 4 " m 
lan hear) 
+ az + lv. . at 


u | 1 ; 
ren 
und bei Einführung von 


« 1 1 1 
2=(2.+4 2 -4+4):W. 
5 1 1 1 1 I 
-[21Z +3 Va +3 Ya v+2]E 5) W 


’ 1 1 1 
D-(V, _ 3) 90 + a) W, 
nimmt die Abschätzung (8.2) die Form an 


A, r 20 . 
= r II 
ea Shut gu de nn (8.10). 


Darin bedeutet (siehe Tab. I.) 
3. Fall: g(&) > 0 und zweimal stetig differenzierbar nO<xr= 1. Wir integrieren ], 
in (8.6) erneut partiell: 


Rh e (a) 3a, 91” sin emx- de. 


04 


Die Schwarzsche Ungleichung, S 9) und (3.2) ergeben jetzt 


1 
Bee = [a(a) (2, 3]? de j sin? gmx de 
0 


'2 


Sen | | [L@@, 8%) ga) Ey, 5) der dr 
0 0 


Nun führen wir wieder L(x, &,&) ein und differenzieren die eckige Klammer aus. Mit Rücksicht 
auf die Sprungrelation der Greenschen Funktion 


a I IE 
ist dabei allgemein 


| IR, 918 zZ = FRWHNOdE baw. i K’@,91& «| 22 


und damit 
N RK’, &) 1) q(&) y.(&, &) «| = FR, 4/) 10 WED) dE— P) u 3) 


Berücksichtigen wir diese Beziehungen bei der zweiten Differentiation, dann erhalten wir für 
den Integranden in (8.11) 


el) htAartrItlr Ir) 
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mit 
= f Kr» 4(&) q(&) Y,(&, 5) dE, A= 2/ K': (©) 4(E) v(E, &) dE, 
E2SR WED, IT K-4’) ao WE) de 
RK SONEHE, = — a) uw. 
Wir setzen in (8.11) ein, multiplizieren aus und schätzen die Integrale einzeln ab. Mit 
f / Kr'%z,&) dede=r 


ist zum Beispiel (abgekürzt geschrieben!) 


SJl= If: 4() q(E) Y,(&, &) «| dx 


SS fRE Ho) a8 ae a 6 9 ar) de 


0 


11 r, 
<P:-[SKrdd=@:r 


ia 0.0 


und ähnlich mit (8.8) 


IAJl=2 


i IK 1a) g(&) y,(&, &) de f Ki (0) 48) 9,6 &) ae) de. 


9 


1 2 1 2 1/2 
<2 [j | I} K? g(x) (8) y,(&, &) 2 a li Re 0) 4 v6 &) «] ae 


1/2 
<2 | FIRE: 06a) 08) de dr [IRA d%@) a) de dr 


2.7 .Yr-Z,. 


IN 
DD 
BIER 

ES 
7) 
en 
un 
17 
I 
N 
Ze 
pa‘ 
Io} 
I 
16) 
"> 


ol: 


Analog ermitteln wir die anderen J,, 2 |f I, J,|: 


| Pr | 220 | Ag yr.v, ad yr-z Be yr-i Pr 


| 


AA an Aay/ ER 7 > 
3 > a | gg F, 297 9°] Ve q 7° Vr- n Va@g':-V, 
| | | @ Mnık ee TE So 
4 "3 49:g°7% 2qq q ra) 2yarg:-ı 
i a 1 N N, 
g:q ’ 0 Vi ®.— 


6 u | j 5 | £ = R" 


Wir schätzen J& mit der Dreiecksungleichung ab, benutzen 


US. BUS azı (8.12) 


a OT 0 5 


und führen ein 


a I an! ns 
=.8 = Mo-Sa Rn 
a 08.76 76 W ) Al 5(6) 945 


ol 


| NM | Pe a I Aa —A ‚ Er a N a u N PN 
— | 4 zZ 200: - Vz V az 1 Yrara ae alg-7°:2, 


Pe 8 ns 


ug 
4 
4 
5 
> 
7 
\ 

ü 

\ 


und damit 


yg Be 
E=sletgfetn), | nes Vger), 


eh fuVz Hr g 1/1 1 
F,— S,\4Z, —Z+-=)|; = - En 17 Be 
| +2]/42 + 5 Ser, s.(vZ v+la z.+)/4 v+4l7 


\ 11 #1 1 , 24. a 
ee nt); 2 s{2+4/4). 


j} ei Ar: 0 gun 
sr vz’y,r 1 «U 2er 3 gi — Pe a a Te 0 
mit 
TnL=_7)79 2 #72 ER?) a; 
en Etrg Ft tg Hr+gg”l, 
IR: q° A | / gig” 2 
++ | | an ir 
; u q 
bei 
Tabelle 2 
r 108 - E, 108- Fr 108 - G, 108. A, 108 - I; 108- K, | 108-2, 108 - M, 
ass 
1 4510 2402 48 1516 230 503 3518 563 
2 663 262 5 232 33 52 355 56 
3 175 55 1 59 se 10 72 11 
4 66 17 0 21 3 3 22 3 
5 35 8 0) 11 1 1 10 | 1 
6 16 3 0 1) 0 0 4 0 
: > : : s Zn N : 
0) 0 0 0) 0 or% De Der 0 
Beispiel: 


Die beiden ersten Eigenwerte /,, /, des Problems 


—y" =Aß2+cosza)y y()=y()=0 


G. Bertram, Fehlerabschätzung für das Ritz-Galerkinsche Verfahren bei Eigenwertproblemen 241 


lassen sich praktisch genau nach der Kettenbruchmethode ermitteln (s. [2], S. 397££.). Wir 
vergleichen sie mit den aus 2- bzw. 5-gliedrigen Ritzschen Ansätzen (8.1) gewonnenen Er- 
gebnissen A,» bzw. A,,, und schreiben die zugehörigen Fehlerspannen A}, 2 bzw. A; für o=], 


II, III nach (8.5), (8.10) und (8.13) dazu an: 


| v=]1 | v=2 

Aa | 4836688 | 21,482 18 
des 1,776 14,416 
Are 0,964 20,920 

III 
a | 0,254 17,60 
Mr 4,836 5 20,325 
Aus 1,016 11,0 
Am: 0,118 | 132 

IIL 
An 0,010 2 2,8 

}y 4,836 46 20,324 5 
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Die besten Schranken A}! liefern also: 
4,8263 <A, = 4,8365, 
17,525 SA, 20,325. 


Ferner stellen wir für r = 2; 4; 6 die Fehlerprozente zusammen in der Form 


lm | 
D ge 
un a et $ 100%: 
nDI 
| vis Ta 9 A,, r 
r | v1 | = 3 
2 2,7% 69,4% 
4 0,2902% 18,24%, 
6 0,0430), 3,5% 


Die erste Tabelle zeigt den Einfluß der Abschätzungsverbesserungen. Beim Übergang 
von der ersten zur zweiten Abschätzung werden die Verbesserungen für },, erkauft durch Ver- 
schlechterungen für A,,; ähnliche Aussagen gelten zum Teil für den Übergang zur dritten Ab- 
schätzung. Die zweite Tabelle zeigt den starken Einfluß der Ansatzgliederzahl r auf die Güte 
der Schranken. Daß bei r Ansatzgliedern der r-te Eigenwert schlecht angenähert wird, ist eine 


generelle Erfahrung. 


Nach den Erfahrungen an anderen Beispielen scheint das vorgeführte Problem die Wirk- 
samkeit der Abschätzungen weder besonders günstig noch besonders ungünstig wiederzugeben. 


9. Eine Klasse von Problemen 4. Ordnung 
Wir behandeln das Problem 
Dre AAO: oe) > 0m 03:31; 
y(0) = y’(0) = y) =’) =. 
Dem vereinfachten Vergleichsproblem 


yY=Ay; y(0) = y’(0) = y(l) = y’(l) =0 
entnehmen wir 


9,=Y2sinenz bei = 0o*.n* =1l,2,.3 
Der Ansatz 


„,= N at-sinonz (at a,-y2) 
e=1 
liefert aus den Galerkinschen Gleichungen für die q, 


1 14 - 
[vet — gun de — 0 e@=1...,n 


die Säkulargleichung 
1 1 l 
Det RR 19 TE a 
0 0 


zur Ermittlung der 4, , fürv=|,...,r. 


Die Fehlerformel lautet zunächst 


mit 


2 / 4) -y2-sinonz- y,z, &) de (vgl (8.3))). 
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1. Fall: q(&) > O0 und stetig nO sr <s1. 
Mit (8.4) erhalten wir 


RL) 
— 4 EEE er RE 9.1 
A Ve »ETHA,, 3 vr ( ) 

bei 
2 1 1 1 
A = — — 
o Pr “RN Di m“ 
Tabelle 3 

r 106 - A, | 10:0, 106. P, |1010. B, 1010. 0,1010 Dy|101: 1100-77 1010 . G,[1010. 7, 1010 - 7,11010. K 


845,128 9 105,39 | 1040,16 | 470 50 115 45 1 146 24 15 


203,505 0 105,80 | 1056,41 50 30 5 5 0,8 | 0,05 6 il - 0,6 
76,764 5 105,82 | 1057,84 1 6 0,1 0,3 | 0,1 | 0,003 0,4 0,07 0,04 


1 300 

2 

3 

E a ee N. (Alle weiteren Koeffi- (Alle weiteren Koeffizienten < 0,1 und 
6 


20,237 4 105,82 | 1058,16 : 
H 3 4 zienten <2 und schnell abnehmend!) 
z.2 3 105,82 Da. schnell abnehmend!) 


© 0 105,82 | 1058,20 | 0 0 0 IR Pr 0 0 0 | ) 
|= 1/9450 | = 1/945 | 
(Die O, usw. werden später erklärt! Buchstaben A,, B, ..... usw. in neuer Bedeutung!) 


2. Fall: g(x) > 0 und stetig differenzierbar nO sr sI1. 
Die Abschätzung läuft wie im 2. Fall des Abschnittes 8, allerdings mit (s. [2], S. 426) 


Be D@+e—29 für z<E$ 


Se 


K(z, &) — 4,4 
o!n 
GR dDE+M— 22) Eur 22 or 
und 
N nee! sinond 
Ka, &) = 2 > Er 
Mit den leicht ermittelbaren Integralen (s. Tab. 3) 
il 13 1 
n : 
2 — ’2 en te Tea «A)y 
Irre da aeg: | |[ees«@« 945 (9.2) 
00 00 
. ie 11 - 
r 1 - r 
0, = Keen —— 5; P,= Ra, &).dede = 
ji Ed I 9 Da 
00 00 
erhalten wir genau wie früher (8.7) mit den J* aus (8.9), allerdings jetzt neu abgeschätzt, bei 
Einführung von 
el 1 
— = — P=S 
Be ee 
und 
B,=[P, u Ve 95T +) Dr 
@=|20P..0, EP. Pe A ! 5 
ln FEAR 95072) 98 945 9450) 
0) 0, : r .S 
alt); 
die Abschätzung 


nen 2 ea 


| Be 
3 At e® Var ij a2 1. . Qu — A,, 
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Darin bedeuten (s. Tab. 3) 
n AN » 
on =glP B.+7g-Co+g®?-Dy. 


3. Fall: g(x) > 0 und zweimal stetig differenzierbar nO<r=sI1. 

Hier kommt im Gegensatz zum 3. Fall des Abschnittes 8 die Sprungrelation der Greenschen 
Funktion K’(E +0,99 —K”’(E—0,&)=1 nicht zum Tragen. Nach einiger Zwischen- 
rechnung erhält man nach zweimaliger Integration die Abschätzung 


1 1 2 
na | || ED no ED art ar. 
a 


Mit den jetzt neu definierten Integralen | | 
nESKE)O EHE, ASK" aa) a WE, 
KO ME DEE,  Is=2[K'd) ae) Ede, 
=] R,d’o) 06) 1.65) dE, ech FRE’) ae) ne, 8) ae. 
erhält der Integrand die Gestalt 


M-)Atht+IS)ta+JI+JIIP- 
Die Beziehungen (9.2) und (s. Tab. 1) 


171 
1 1 
9 she. RB = Dt en 
I de de= on; Zr en Sam: T,= 9450 0, 
00 


erlauben ähnlich wie oben Abschätzungen, die hier jedoch nicht angeschrieben werden sollen, 
und welche die Einführung der folgenden Größen nahelegen: 


E,=(V, + sv ar n)' % 

Fr=[4 Pa 2\/>, en 4, u) -T, 

(0 ea 3-0 Ri | 2 

1:8 ar YV,-P,+ 22 . . + Ve-% + eo: u) 
n=[2 90.4 Vr- rn A Vo. + 2 a) Z 


en i 
K,-[4yP,:6,+ |/p.: Ri 
| / + )/ 9450 7 jo Zeh 945 a) 


Die endgültige Abschätzung lautet dann 


A: En: III 


ge Br y s A, es nr“ a ee — a Di PR Re ee u (9.4) 


mit (s. Tab. 3) 
ae RA.E an «, 8 
gie E + Fr +4": +99 HM, +40, +gQgt-K). 
Beispiel: 
Für 
y!v —/ (2 Be x) T, y(0) — y’’(0) — y(1) = y’(1) Ef] 
ergibt der Ansatz 


3 
I . ee: | - 
Im=2a sinonx (,- v2) 


» 
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aus den Galerkinschen Gleichungen 


1 
[snonsi AR — Di) de=0 (e==1,2, 3) 


die Säkulargleichung (# = .) : 


A 
3 16 
Er = 
2 Ir x 
16 3 48 
Fr ne , er Er 
Ir Hl 35% 5) 
48 a 
en y 4, 
0 5 5 + 81mtı 
Ihre Lösungen A,,; und die zugehörigen Fehlerschranken nach (9.1), (9.3) und (9.4) sind 
vV=]1 v=2 v=3 
1,5 64,324 1291,8 | 4354 
A: 0,61 211 1754 
A; 0,055 1130 4360 
Ans | 0,007 613 | 4261 


r 1 1 2 
Setzen wir etwa /), = (r = Ai.) Fr > 4A,y,,, dann ist z.B. 


64,019 + 0,305 (+ 0,48% Fehler) für o=1 
%ı = 164,296 + 0,028 (+ 0,043% Fehler) für = II 
64,320 + 0,004 (+ 0,006% Fehler) für o= II. 


Die Verbesserungen lassen sich ähnlich diskutieren wie im 8. Abschnitt. Man überschlägt 
leicht, daß bei Hinzunahme nur eines weiteren Ansatzgliedes die Fehlerprozente für A, im Falle 
o=1l etwa mit ;sc= II mit a undo= III mit & multipliziert werden. 


10. Ausblick 
In ähnlicher Weise lassen sich Fehlerformeln für weitere Problemklassen bereitstellen, 
auch für partielle Differentialgleichungen. 
Die Ansatzfunktionen fallen leider oft unhandlich aus. Schon für das Problem 


yY=igda)y gm >O0 in Osısi 
y(0) = y’(0) = y’() =’) =0 
sind die Ansatzfunktionen dem Problem 
„W=iy; 30 =-y0)=y’W)=y"l)=0 
zu entnehmen. Der Säkulargleichung 
cos YA - Coj YA =—] 
entnimmt man 
12419; 14862; weil, = 104l0E. .. 
Setzen wir 
s,=sin Yo; 6, c0s Yo; & = Sin Yo; GC, Coj Yo, 


und 


SEO, Yeltu 


; 1 1 an, 1 1 4, — 
a E 1 I &+ )— = (2 ct s) +% © &,— ) +6 @ SG — s) -r je! 
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dann wird 


ZAMM 39 (1959) Heft 5/6 


98) = d, y. [Sin Y @, 2 — sin Vo, | —o, [Co Yo, x — cos Vo, «| i 


Ferner ist (s. [2], S. 426) 


1m@E—a) fürase 
Ka, &) = Ä 
78:9) für st 


o 
ll 
fen 

o 


Zuweilen läßt sich die explizite Ermittlung von K(x, &) dadurch umgehen, daß man die o, nach 
unten, also die Summe nur nach oben abschätzt. Eine Skizze zur Säkulargleichung bietet dafür 


oft einen Anhalt. 


Die Fehlerformel lautet im einfachsten Falle 


ER 


mit 
AN 


I) 


0:= 
bei 


u! 1 
ed. 
0e= 


ine ae U 
1 + A, z Q! 2 
Tabelle 4 
f | 106 Ar 
| 

1 | 2813 

22 | 756 

a 494 

4 398 

(0) 0 


Man überschlägt leicht, daß A, bei g=1 sich im Falle r=4 etwa mit einem Fehler 


von + 0,25% wird angeben lassen. 


Vielfach wird es zweckmäßig sein, vor der Durchführung der Rechnung einer derartigen 
Überschlagsrechnung die Zahl der Ansatzglieder zu entnehmen. 


Manuskripteingang: 29.7.1958 


Anschrift: Dozent Dr. G. Bertram, Hamburg 39, Krohnskamp 72, I 


KLEINE MITTEILUNGEN 


Stagnation Point Pressures on Symmetrical 
Blunt-Nosed Bodies Entering an Infinite In- 
compressible Ocean 


In reference 1, a method is presented for solving 
the wedge entry problem. One application described 
therein is an approximate method for determining 
the pressures on bluntnosed bodies entering an in- 
finite imcompressible ocean. These may be two or 
three dimensional, symmetrical or unsymmetrical 
bodies, at either constant or variable velocity. A 
brief outline of the method follows: 

Consider a straight sided, infinitely long, two 
dimensional wedge moving at constant vertical 
velocity V and entering an infinite ocean. (A theory 
identical to that which follows holds also for three 
dimensional wedges.) Equivalent in every way, but 
more convenient for our purposes, is the case of the 
wedge stationary and the ocean moving vertically 
upward with constant velocity V, and engulfing the 
wedge. 

In the physical x, y, t, plane, this unsteady pheno- 
menon maps as shown in Figure la for times i = (, 
t=t, t=t. If we introduce new coordinates, 


then in the (£,n) plane the phenomenon 


is shown in the single map, as in Fig. 1b. From this 
point on all statements refer to the (£,n) plane, 
which, in analogy to a similar transformation with 
coneshapes, we call the „conical‘“ plane. 

Because of continuity, the following two conditions 
hold, 

a) The fluid displaced by the wedge must appear 
above the undisturbed free streamline position, or 


Volume0ab =Volumbcoo ....(l) 


b) The lenght of the free surface line remains 
constant with time, or 


Lengthaboo = Lengthecw ....(2). 


A fundamental quantity in this formulation is 
the „reduced velocity, Q“ defined as 


Dir. en) 
in which 


q isthe actual physical velocity at a point in the 
(&,n) plane 
r is the distance, in the (&,n) plane, of this point 
from the origin. 

Because the flow is irrotational, a potential func- 
tion, 9, exists, such that 


-7=- wu =40=(4,2,)0 + (4). 
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In addition, it is shown that (because of continuity) 


A 02 0° 
ro=tzn)e=0 0. Bet) 
everywhere on the field, and 
FT HR 
= + DE Dane ae (6) 


everywhere, with p and o the pressure and density 
respectively. ; 


Y. 


(2) (b) 
Fig. 1 


Finally it is shown that on the F.S., at infinity 


DEN N TEBREE (T): 
There are many other relations and properties 
which hold for this formulation of the problem — 
. see Ref. 1. However, for our purposes, the above are 
sufficient, and we summarize them as follows: 
The solution for the unsteady flow problem of 
Fig. 1b requires finding a function, g(£,n), which 
satisfies 


(a) V?p9=0, everywhere in the field 


(b) op 0 along the wedge face, where n is 
on the direction of the normal to the 
wedge 
op Pe 
(e) Free V at infinity 


4 _@?® r? along the free streamline (since 
(d) P=77059p= 0 here). Hence 
2 


(e) ee on the F.S. at o. 


g 


Uatoo 


| 
| en 
Fig. 2 


It is possible to obtain p functions which satisfy 
all of the above except (d) by utilizing the potential 
flow solutions for classical, irrotational steady state 
motion. These may be taken as approximate solu- 

tions to the problem (2). Having 9, using equations (4) 


c p=0 
ree Streamline 


- point and this is used for AR < about 


and (6) we can obtain the pressure acting on the 
body. The method may also be used to give approxi- 
mate solutions for blunt-nosed shapes — simply apply 
it to the corresponding potential function. 

One complete family of blunt-nosed bodies may be 
obtained from the source-sink pair in a steady flow 
stream (reference 3), see Fig.2. When a =(, the 
shape is a circle, and when a — b, the shape becomes 


Circle 
U atoo 


Stagnation Point 


8 
L 
s 
Q 


Fig. 3 


an elongated cigar-type. If we define the aspect 
ratio ARin terms of the undisturbed freesurface level 
see Fig. 3, then this single family of shapes will 
include the complete range, 


Verlorene) 


For AR<]1 (corresponding to very small immer- 
sion), the flat plate potential solution is very nearly 
equal to the blunt-nosed solution near the stagnation 


= . 
For AR> A the source-sink potential function 


solution is used (4). 


Computations using the above noted methods then 
indicate that the approximate stagnation point 
pressure ?9s, for twodimensional, symmetrical, con- 
stant-velocity entry is given by the single equation, 


u 2 
: een N 
eV 
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Das Hertwigsche Iterationsverfahren zur Auf- 
lösung linearer Gleichungssysteme als Ge- 
samt- und Einzelschrittverfahren 


Bekanntlich sind die meisten Iterationsverfahren 
zur Auflösung eines linearen Gleichungssystems 


ae Da 


von der Form 
Ber) HE dev, B+E=N, det B-+0 (2). 


Hierbei stellt 2 — (=?) die v-te Näherung für eine 
exakte Lösung r = (x;) dar. 

Zur Konvergenz des Verfahrens ist im Falle 
det X = 0 notwendig und hinreichend, daß die Folge 
der Potenzen der Matrix 


& = (,)=—-HE=-E- DEU . 
gegen die Nullmatrix konvergiert, wofür z. B. 
G=M << 

ae 


. (8) 


gu, ER IR 


hinreichend ist (Zeilensummenkriterium). Ist es er- 
füllt, so erhält man bekanntlich in 


Mer (aa © 


eine rekursive Fehlerabschätzung. 


Für die Konvergenz des Verfahrens (2) wird es 
nach (3) also nötig sein, daß (mit einer Formulierung 
von Wittmeyer) die Matrix 8 der Matrix \ „hin- 
reichend benachbart“ ist (d.h. 8914 > &). Im Falle 
B—= Da = (dr air) mit dem Kroneckersymbol öik, 
d.h. beim üblichen Gesamtschrittverfahren, wird das 
der Fall sein, wenn die Hauptdiagonalelemente von X 
groß gegen die restlichen Elemente sind. 

Das in der Literatur seltener erwähnte Verfahren 
von Hertwig besteht nun darin, aus der Matrix WU 
durch Weglassen kleiner Elemente eine neue Matrix ® 
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Matrix X noch besser „‚benachbarte‘‘ Matrix sein als 
z.B. Da. In der Tat: ist z. B. die Anzahl » der Glei- 
chungen und Unbekannten gerade, und wählt man 
die W; von der Ordnung 2, so folgt, wenn man mit 


Anm = %r m — um %k 


die Unterdeterminanten zweiter Ordnung der Matrix Di 
bezeichnet, 


12 dus, EEE 

Da 1225 X Zr 

j Ku Q a: 

a a Bu xx0 1 x x 
A,ani - . u 

Be 10 

ee | 


G=i+(- Yin) 


9-14 stimmt also nicht nur (wie beim üblichen 
Gesamtschrittverfahren) in der Diagonalen mit € 
überein, sondern sogar in den Diagonalkästchen. 
Überwiegen die Diagonalelemente von Y, so sind auch 
hier die restlichen Elemente von 812% klein. Das 
Verfahren wird daher im allgemeinen besser konver- 
gieren als das Gesamtschrittverfahren. (4) liefert 
auch für diese Methode wieder ein hinreichendes 
Konvergenzkriterium und (5) eine Fehlerabschätzung. 
Die gegenüber dem üblichen Gesamtschrittverfahren 
größere Vorarbeit vor dem Beginn der eigentlichen 
Iterationsrechnung (sie erfordert etwa n? Additionen 
und Multiplikationen mehr) wird bereits ausgeglichen, 
wenn man zur Erzielung einer gewissen Genauigkeit 
beim Hertwigschen Verfahren mit nur einem 
Iterationsschritt weniger auskommt. 

Ausführlich geschrieben ergibt sich x” aus 7» 
beim Hertwigschen Verfahren (Ordnung der 9; 
gleich zwei) in folgender Weise: 


7 1 1) e-1) er 
a) = ee | . an * +A,s02%  tAase a tt A naar, 
—1 1) (—1) : 1 

y=7, a = = Tr « Ay +4, 411% Heise )+w, 
o_ 1 1), „ol ‘ e ES 1 

er Pewi UL PIEWFIGE Gute up. FEW REE Ge u r E teten) tn 
6) 1 —ı —1) A Ar DE 

4.227 As + A933 % 4 E T Aunss )+w 
) _ ee Ba . „e—1) „el 

Fa An, n, n—1,n—1 in un, eg 1 ni Ana Alf ne } +4 n,3,n—1,n—ı "73 


derart zu bilden, daß (nach evtl. Umnumerierung 
der Gleichungen des Systems (1)) ®B aus diagonal 
aneinandergereihten quadratischen Untermatrizen von 
U besteht und sonst nur Nullen enthält: 


A 0 
U, 
B= 
0 Uni 
Dann folgt 
Dr 0 
1 n Ur 
0 UzE, 


Diese Matrix läßt sich leicht berechnen, wenn die 
Ordnungen der W; nicht zu groß, etwa nicht größer 
als 2 sind. ® wird nun im allgemeinen eine der 


w—1 
ur An ' u > *) + Um 


mit (w) = Bd. 
Wesentlicher als der Vergleich des Hertwigschen 
Verfahrens mit dem üblichen Gesamtschrittverfahren 
ist der mit dem Seidelschen Einzelschrittverfahren. 
Die Erfahrung lehrt, daß die Seidelsche Methode 
der Hertwigschen überlegen ist. f 

Es liegt nun nahe, das Hertwigsche Verfahren, 
das ja auch in Gesamtschritten fortschreitet, nach 
dem Muster des Seidelschen Verfahrens in Einzel- 
schritten durchzurechnen. 

Bei dieser Methode lautet die Iterationsvorschrift 
in Matrizenschreibweise: 


BBIMEN+ EBEN — Ede, 


wobei ® die vorher beschriebene Matrix des Hert- 
wigschen Verfahrens bedeutet, und unter ©, (bzw. 
©;) eine Matrix verstanden werden soll, die oberhalb 
(bzw. unterhalb) und in der Hauptdiagonale die 
Elemente von & und sonst Nullen enthält. Auch 


. k E 
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dieses Verfahren ist von der Form (2), so daß (4) 
wiederum ein Konvergenzkriterium und (5) eine 
Fehlerabschätzung liefert. 

Verfasser, dem aus der Literatur kein Hinweis auf 
dieses abgewandelte Hertwigsche Verfahren und 
etwaige Erfahrungen mit ihm bekannt ist, hat damit 
oft schnellere und bisher bei noch keinem Beispiel 
nennenswert langsamere Konvergenz festgestellt als 
bei der Seidelschen Methode. 
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H. Wittmeyer, Über die Lösung von linearen Gleichungssystemen 
durch Iteration. Z. angew. Math. Mech. 16 (1936), S. 301—310. 


H. Sassenfeld, Über lineare Gleichungssysteme und ihre Anwen- 
dung bei Summenverfahren für Rand- und Eigenwertaufgaben. 
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Ein Sicherheitscode für Fernschreibgeräte, 
die zur Ein- und Ausgabe an elektronischen 
Rechenmaschinen verwendet werden 


Neben dem internationalen Fernschreibeode können 
bei Rechenmaschinen auch andere Codes für Loch- 
streifen angewandt werden, falls man darauf ver- 
zichtet, die Maschinen an das Fernschreibnetz an- 
zuschließen!). Insbesondere wurden Codes mit 
tetradenverschlüsselten Dezimalzahlen verwendet. 

Sinnvollerweise wird man an einen speziell für 
Rechenmaschinenzwecke zu entwickelnden Code die 
Forderungen stellen: 

1. Er soll ein Sicherheitscode sein, d.h., die den 
zehn Ziffern zugeordneten Lochkombinationen sollen 
sich jeweils in mindestens 2 Dualstellen unterscheiden, 
wodurch Fehler im Ablesegerät, im Ein- und Aus- 
gangspufferspeicher und bei der Ausgabe weitgehend 
kontrolliert werden. 

2. Der (formale) Dualwert einer Lochkombination 
soll sich aus der dargestellten Ziffer durch eine ein- 
fache Formel ergeben, wodurch sich je nach Rechen- 


anlage entweder für die Entschlüsselungsgeräte tech- 


1) Fernschreibmaschinen, die nach einem anderen als dem In- 
ternationalen Code arbeiten, werden von der Post nicht zum An- 
schluß an das öffentliche Fernschreib-Netz zugelassen. 
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nische Vereinfachungen ergeben, oder die Konver- 
tierungsgeschwindigkeiten gesteigert werden können. 

Man könnte sich vorstellen, daß die Verwendung 
eines solchen Codes auch auf anderen Gebieten als 
der Rechenmaschinentechnik Vorteile aufweisen 
würde. 

Wir haben noch die weitere Forderung gestellt, daß 

3. Leerband (Dualwert 0) und Buchstabenumschal- 
tung (Dualwert 31) nicht als Bilder von Ziffern auf- 
treten sollen. 

Die Diskussion zweier möglicher Wege zur Kon- 
struktion eines den obigen Forderungen genügenden 
Fernschreib-Codes — Quersummenverfahren und 
Kongruenzenverfahren — ergab, daß der folgende den 
Vorzug verdient: Die Ziffern werden als Loch- 
kombinationen mit dem Dualwert 


=zkı+b (mod 32) 


dargestellt. Bei ungeradem %k wird die Rücktrans- 
formation geleistet durch 


x=kt(y—b) (mod 32). 


Die Forderungen 1. und 3. ergeben eine Ein- 
schränkung der zulässigen Werte für k und 5 auf die 


. folgenden Möglichkeiten 


K=3 b>=1,2,.3,18 
k=—3 b = 13, 28, 29, 30. 
Unter diesen sind 
yz=z370+2 und = 30 —3x 
c=11ly+10 x=10—11y 


durch das Auftreten der 10 in der Umkehrungsformel 
besonders ausgezeichnet, weil hierbei die Rücktrans- 
formation auch gleich ohne weitere Rechenoperation 
zu einer Unterscheidung zwischen Ziffer und Nicht- 
ziffer führt. 

Herrn Th. Fromme von der Fa. Zuse K.G. bin 
ich für die Anregung des Problems zu Dank ver- 
pflichtet. -Die Fa. Siemens und Halske AG hat sich 
freundlicherweise bereit erklärt, für die Fernschreib- 
maschine T 100 besondere Wählschienen anzufertigen, 
mit denen der Code y = 30 — 3 x realisiert wird. 


Verfasser: Dr. E.Nuding, i. Hs. Siemens-Schuk- 
kertwerke AG., Reaktor-Entwicklung, 
Erlangen, Schillerstr. 62 


BUCHBESPRECHUNGEN 


Dr.-Ing. R. Zurmühl, Matrizen. Eine Darstellung 
für Ingenieure. Zweite, völlig neubearb. Aufl. XV + 
467 S. m. 76 Abb. Berlin/Göttingen/Heidelberg 1958. 
Springer-Verlag. Preis geb. 33.— DM. 

Der Verfasser legt mit der zweiten Auflage seiner 
„Matrizen‘‘ ein fast völlig neues Buch vor. Schon 
ein flüchtiger Vergleich mit der ersten Auflage zeigt, 
daß die Darstellung der Matrizentheorie und des 
Matrizenkalküls durch die gründliche Neubearbeitung 
noch gewonnen hat, sich noch mehr als bei der ersten 
Auflage durch große Klarheit und Eindringlichkeit 
auszeichnet und von einem hervorragenden didak- 
dischen Einfühlungsvermögen des Verfassers zeugt. 
Es ist nicht nur für den Ingenieur und Physiker 
wesentlich, sondern auch für den Mathematiker 
nützlich, wenn die mathematischen Begriffsbildun- 
gen aus Fragestellungen der Anwendungsgebiete 
heraus entwickelt werden, wie dies an vielen Stellen 
des Buches geschieht. Im ersten Kapitel werden 
die Grundbegriffe und der Matrizenkalkül bis zur 
Algebra der linearen Abbildungen und Transfor- 
mationen entwickelt, wobei auch die komplexen 
Matrizen berücksichtigt werden. Das zweite Kapitel 
bringt — in interessanter Weise vom Gaußschen 
Algorithmus ausgehend — die Theorie der Gleichungs- 
systeme und dringt bis zu den Orthogonal- und 


Biorthogonalsystemen vor. Durch den Paragraphen 
über Polynommatrizen und die Transformation auf 
die Smithsche Normalform wird die Elementarteiler- 
theorie vorbereitet, die im fünften Kapitel über die 
Struktur einer Matrix eine besonders liebevolle und 
didaktische wohldurchdachte Darstellung gefunden 
hat. Die dazwischenliegenden Kapitel III und IV 
gehen auf quadratische Formen und ihre Anwen- 
dungen in der Ausgleichsrechnung, bei der Berechnung 
von Fachwerken und bei der Behandlung von Schwin- 
gungsaufgaben und sehr ausführlich auf das Eigen- 
wertproblem ein. Der Einstellung des Verfassers 
entsprechend werden die den Matrizenkalkül mit 
großem Vorteil benutzenden numerischen Verfahren 
zur iterativen und zur direkten Lösung von Eigen- 
wertaufgaben und zur iterativen Behandlung linearer 
Gleichungssysteme mit großer Gründlichkeit darge- 
legt. Das letzte Kapitel, das fast ein Viertel des ganzen 
Buches umfaßt, ist den Anwendungen der Matrizen- 
rechnung in der Elektrotechnik, in der Statik, in der 
Elastomechanik und in der Schwingungstechnik ge- 
widmet. 

Es führt bis zu dem modernen Problem der Be- 
handlung von Systemen linearer Differentialglei- 
chungen. Gerade dieses Kapitel zeigt, mit welcher 
Sorgfalt und Liebe sich der Verfasser immer wieder 


18 


Aa a 

, ee) ’R H a = 
; a 250 r 
Be, 
ex _ mit der überaus wichtigen Frage beschäftigt, wie 
mathematische Theorien in zweckvoller Weise — da- 


: Ertl 
NE 


Ye 


rauf darf wohl der Ton gelegt werden — den Anwen- 
dungen dienstbar gemacht werden können. i 
Viele — meist bis zum zahlenmäßigen Endergebnis 


durchgeführte — Beispiele und zahlreiche übersicht- 


liche Schemata erleichtern die Lektüre dieses wert- 
vollen Buches, das sich auch ohne besondere Empfeh- 
lung einen großen Freundeskreis schaffen wird. 


Dresden H. Heinrich 


J. W. Dungey (Lecturer in Mathematics at King’s 
College, Newcastle upon Tyne), Cosmic Electro- 
dynamics. IX + 183 S. m. Abb. Cambridge 1958. 
Cambridge University Press. Preis geb. 32s. 6d. 

Wie alle Monographien dieser Veröffentlichungs- 
reihe stellt das Buch einen kurzen Abriß des betreffen- 
den Fachgebietes dar. Bei der Darstellung geht es 
dem Verfasser vornehmlich um die Klärung und 
Darstellung der theoretischen Grundlagen und weiter 
um eine Darlegung der physikalischen Gesichts- 
punkte und mathematischen Methoden. Zu detail- 
lierten Anwendungen läßt der enge Rahmen nur 
wenig Raum. Behandelt werden hauptsächlich die 
elektromagnetischen Erscheinungen auf der Sonne, 
magnetische Stürme und das Nordlicht, ferner werden 
stellare und interstellare Felder behandelt, die Be- 
schleunigung von Höhenstrahlteilchen und elektro- 
magnetische Probleme der Ionosphäre. Das ganze 
Gebiet, obwohl rein klassischen Ursprungs, ist noch 
verhältnismäßig neu und befindet sich in rascher 
Entwicklung. Es hat insbesondere im Zusammen- 
hang mit den Bemühungen um kontrollierbare ther- 
monukleare Reaktionen eine Bedeutung, auch für 
viele Gebiete der modernen Laboratoriumsphysik und 
dient dem Zweck, solche Physiker anderer Fachgebiete 
verhältnismäßig rasch zu informieren. 


Dresden 


W. Macke 


L.D. Landau and E.M.Lifshitz, Quantum Me- 
chanics, Non-Relativistic Theory. Volume 3 
of a Course of Theoretical Physics. Authorized Trans- 
lation from the Russian by J. B. Sykes and J. S. Bell. 
XII + 5158. m. 51 Abb. London-Paris 1958. Per- 
gamon Press. Preis geb. 80,— s. 

Für eine Darstellung der Quantenmechanik bieten 
sich wegen der bekannten Mittelstellung des Quants 
zwischen dem Korpuskel einerseits und der Welle 
andererseits im allgemeinen drei Wege an, nämlich 
beim Teilchenbild oder Wellenbild zu beginnen und 
die Quantentheorie korrespondenzmäßig zu erfassen, 
oder schließlich drittens, die Quantenmechanik aus 
sehr allgemeinen Erfahrungen und Postulaten heraus 
zu entwickeln, ohne sich dabei auf das Korrespondenz 
prinzip zu stützen. Dieser dritte Weg ist hier be- 
schritten. Die Darstellung beginnt mit den Prin- 
zipien der Unbestimmtheit und der Superponierbar- 
keit quantenmechanischer Zustände und entwickelt 
von hier aus deduktiv die quantenmechanischen 
Erscheinungen bis in ihre letzten Details. Diese 
streng logische Deduktion der Quantenmechanik in 
hohem Niveau liefert eine Darstellung von großer 


ec Bi 73 Seue 
Geschlossenheit. und Zwa keit, nicht 
allerdings auch vom Leser ein ge\ iveau 
auszusetzen. Dem Theoretiker wie auch dem 
schlagenden Experimentator wird das Buch iı 
Fragen eine wertvolle Stütze sein. _ % 
Der Inhalt konzentriert sich auf die nichtrel ’ 
scheQuantentheorie. Das Buch stellt den dritten Band. hr 


einer Reihe von insgesamt neun Bänden der 
Verfasser dar. Die einz« Bände dieser 1 
behandeln 1. Mechanik, 2. Klassische Feldt 
3. Quantenmechanik, 4. Relativistische Quante 
rie, 5. Statistische Physik, 6. Mechanik der F 
keiten, 7. Theorie der Elastizität, 8. Elektrodyna: 
kontinuierlicher Medien, 9. Physikalische Kinetik. 


Dresden W.Mäcke 


Neor} 


Grenzschichtforschung — DBoundary Layer 
Research. Symposium Freiburg/Br. 26. bis 29. Au- 
gust 1957. Herausgegeben von H. Görtler. XI + 
411 S. m. 206 Abb. Berlin/Göttingen/Heidelberg 1958. Br 
Springer-Verlag. Preis geb. 67,50 DM. R 

Der vorliegende Berichtsband enthält sämtliche 
Vorträge und die wichtigsten Teile der zugehörigen x 
Diskussionen, die auf einem von der Internationalen 
Union für theoretische und angewandte Mechanik 
(IUTAM) veranstaltetem Symposium über Probleme 
der Grenzschichtforschung in Freiburg/Breisgau ge- 
halten wurden. Thematisch bedecken die insgesamt 
32 Vorträge den gesamten Bereich von der klassischen 
Grenzschichttheorie bis zu den modernsten Pro- 
blemen der Strömung in verdünnten Gasen und in 
Medien, die aus verschiedenen Komponenten be- 
stehen. Einen großen Teil des Bandes nehmen Unter- 
suchungen zur Stabilitätstheorie der laminaren Grenz- 
schichten ein. Hier sind sowohl von der theoretischen 
als auch von der experimentellen Seite her beacht- 
liche Fortschritte erzielt worden, über die auf diesem 
Symposium wohl zum ersten Male in dieser Ausführ- 
lichkeit berichtet wurde. Behandelt werden ferner 
Temperatur-Grenzschichten, Scherströmungen, drei- 
dimensionale Grenzschichten, ähnliche Strömungen 
und nichtlineare Effekte — ganz abgesehen von den 
zahlreichen Problemen, die vor allem in den erfreu- 
lich gestrafft wiedergegebenen Diskussionen ange- = 
stoßen werden. 

Etwa 80 Teilnehmer von 17 Nationen aus Ost und 
West waren zum Symposium zusammengekommen. 

Als Ergebnis liegt nun ein Berichtband vor, der den 
Stand der Grenzschichtforschung so widerspiegelt, 
wie dies überhaupt in freier, Länder und Grenzen 
überschreitender Kommunikation erreichbar scheint. 
Der eingeweihte Leser wird unschwer erkennen kön- 
nen, daß diese großartige Übersicht vor allem auch 
der planenden, ordnenden und klärenden Arbeit des 
Herausgebers, der gleichzeitig der wissenschaftliche 
Leiter des Symposiums war, zu danken ist. Für den 
auf dem Gebiet der Grenzschichtforschung tätigen 
Wissenschaftler — sei er nun theoretisch oder experi- 


mentell orientiert — ist dieser Berichtsband eine 
Fundgrube. 
Stuttgart K. Magnus x 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


NEUE ZEITSCHRIFTEN 


1. A Nehezipari Müszaki Egyetem Közlemönyei. 
Mitteilungen der Technischen Universität für das 
Schwerindustriewesen, Miskole/Ungarn. 

Der Schriftleitung liegen Band 1 (1957) und Band 2 
(1958) vor. An Aufsätzen, die den Leserkreis der 
ZAMM interessieren, sind enthalten 
in Band l: 

Dr. G. Petrich: Graphische Lösung der kürzesten 
geradlinigen Verbindung dreier windschiefer Achsen. 
Dr. G. Gäspär: Die axiomatische Begründung der 

Determinantentheorie. 


I. Raisz: Interpolationspolynome. 

L. Huszthy: Bestimmung der minimalen Beschleuni- 
gungszeit rotierender Massen. 

M. Hosszu: Gemeinsame Verallgemeinerung einiger 
Funktionen von mehreren Veränderlichen. 

Dr. A. Nikod&musz: Über die graphische Lösung par- 
tieller Differentialgleichungen. 

E. Vineze: Über eine Verallgemeinerung des Begriffes 
der monotonen Folge. 

Dr. I. Sälyi: Die Erregungsschwingungen gedämpfter 
Systeme. 
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I. Huszär: Die optische Spannungsprüfung im Dienste 
der wissenschaftlichen Forschung und der Industrie, 


in Band 2: 


Dr. I. Sälyi: Funktionen mit komplexen Veränder- 
lichen in der theoretischen Elastizitätslehre. 

R. Vankö: Die Verwendung der allgemeinen Lagrange- 
schen Gleichungen zur Berechnung von Hebezeugen. 

Dr. S. Geleji: Eine neue Theorie des Walzverfahrens. 

Dr. M. Hosszü: Die Auslegung und Charakterisierung 
der Eulerschen homogenen Funktion in einem all- 
gemeinen algebraischen Saystem. 

Dr. G. Petrich: N. I. Lobatschewskij’s Geometrie. 

Dr. Z. Terplän: Die Entwicklung eines viergelenkigen, 
ebenen Getriebes für vorgeschriebene Bewegungs- 
gesetze mit komplexen Veränderlichen nach der 
Blochschen Methode. i 


Sämtliche Aufsätze sind in ungarischer Sprache ge- 
schrieben. 


2. PMM — Journal of Applied Mathematics and Me- 
chanies. Vol. 22, No. 1, 1958. Pergamon Press Ltd., 
London. Preis £ 12.19.0. per volume, bi-monthly. 

Bei dieser Zeitschrift, deren erstes Heft vorliegt, 
handelt es sich um eine englische Übersetzung der 
sowjetischen Zeitschrift ‚„Prikladnaja Matematika i 
Mechanika‘‘, die vom Institut für Mechanik, Abtei- 
lung Technische Wissenschaften, der Sowjetischen 
Akademie der Wissenschaften zwei-monatlich her- 
ausgegeben wird und seit zwei Jahrzehnten eine der 
führenden sowjetischen Zeitschriften auf diesem Ge- 
biete ist. 


3. Numerische Mathematik. Herausgegeben von 
R. Sauer, E. Stiefel, J. Todd, A. Walther. Berlin/Göt- 
tingen/Heidelberg 1959. Springer-Verlag. Preis pro 
Band maximal 100,— DM. 

Auf dem Programm dieser neuen Zeitschrift steht 
die Veröffentlichung von Arbeiten, die sich mit all- 
gemeinen Problemen des digitalen Rechnens, mit der 
Diskussion bestehender und der Entwicklung neuer 
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numerischer Verfahren beschäftigen. Die numerischen 
und programmierungstechnischen Gesichtspunkte des 
Einsatzes von Rechenautomaten werden dabei im 
Vordergrund stehen. Auch Aufsätze aus dem Gebiet 
der Informationstheorie sollen aufgenommen werden. 
Die Persönlichkeiten der Herausgeber und die Zu- 
sammensetzung des Kreises, der bei der Herausgabe 
mitwirken wird, gewährleisten, daß sich die Arbeit. 
der Zeitschrift auf breiter internationaler Grundlage 
entfalten und sich auf einem hohen wissenschaftlichen 
Stand befinden wird. Rz 

Es liegt bereits das erste Heft des ersten Bandes Bii 
vor mit Aufsätzen von 


Stiefel, E., Über diskrete und lineare Tschebyscheff- 
Approximationen; 0 

Householder, Alston S., and Friedrich L. Bauer, On KR 
certain methods for expanding the characteristic 
polynomial; 

Weisfeld, Morris, Orthogonal polynomials in several Mi. - 
variables; Du 

Perlis, A. J., and K. Samelson, Report on the Al- Rx 
gorithmie Language ALGOL. s 
Die Hefte werden in zwangloser Folge erscheinen 

und je 5Hefte zu einem Band zusammengefaßt 

werden. 


4. Mathematik- Technik -Wirtschaitt (MTW). re 

Zeitschrift für moderne Rechentechnik und Auto- 
mation. Herausgeber: R. Inzinger, Wien. Stiasny- 
Verlag Graz. Erscheinungsweise: 4mal jährlich. 

Die bisher vom Mathematischen Labor der Tech- 
nischen Hochschule Wien herausgegebenen M-T-W- 
Mitteilungen werden vom Beginn des Jahres 1959 an 
als Zeitschrift erscheinen, die in Zusammenarbeit mit 
wissenschaftlichen Instituten und Gesellschaften so- 
wie Körperschaften und Verbänden aus Technik und 
Wirtschaft gestaltet wird. Dem Herausgeber und der 
Schriftleitung stehen als Mitwirkende namhafte Per- 
sönlichkeiten aus Österreich, aus Deutschland und 
der Schweiz zur Seite. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


W. Rehwald, Elementare Einführung in die 
Bessel-, Neumann- und Hankel-Funktionen. 
(Mathematische Funktionen in Physik und Technik. 
Herausgeg. v. Prof. Dr.-Ing. O. Zinke. Erster Band.) 
III + 478. m. 19 Abb. Stuttgart 1959. S. Hirzel 
Verlag. Preis kart. 6,30 DM. 


Anwendung von Rechenmaschinen bei der 
Berechnung von Regelvorgängen. Vorträge ge- 
halten bei einer Tagung des Fachausschusses Re- 
gelungsmathematik der Gesellschaft für Angew. 
Mathematik u. Mechanik (GAMM) in Düsseldorf am 
8.11.1957, zusammengestellt von W. Oppelt. 
128S. m. 121 Abb. München 1958. R. Oldenbourg 
Verlag. Preis geb. 16,80 DM. 


F.E. Hohn (Ass. Prof. of Mathematics, University 
of Illinois), Elementary Matrix Algebra. XI+ 
3058. New York 1958. The Macmillan Company. 
Preis geb. $ 7,50.. 


W.W. Solodownikow, Grundlagen der selbst- 
tätigen Regelung. Band I: Allgemeine Grund- 
lagen der Theorie linearisierter selbsttätiger Rege- 
lungssysteme. Deutsche Bearbeitung unter Prof. Dr. 
H. Kindler. XVI + 727 S. m. 460 Abb. u. 23 Tafeln. 
München 1959. R. Oldenbourg Verlag. Preis geb. 
65, — DM. 


L.Robin, Fonctions Spheriques de Legendre 
et Fonctions Sphö6roidales. Tome II. VIII + 
384S. m. Abb. Paris 1958. Gauthier-Villars. Preis 
brosch. 5000 F (alt). 


D. Hilbert f und W. Ackermann, Grundzüge der 
theoretischen Logik. (Grundlehren der Mathe- 
matischen Wissenschaften, Band XXVII) 4. Aufl. 
VIII +183S. Berlin/Göttingen/Heidelberg 1959. 
Springer-Verlag. Preis geb. 36,60 DM. 


J. H. Weihnacht, Prinzipien zur Lösung 
mathematischer Probleme. VIII + 1168. m. 
45 Abb. Braunschweig 1958. Verlag Friedr. Vie- 
weg & Sohn. Preis kart. 8,80 DM. 


Dr. E. Burger (apl. Prof. d. Mathematik a. d. Uni- 
versität Frankfurt/Main), Einführung in die 
Theorie der Spiele. 169 S. Berlin 1959. Walter 
de Gruyter & Co. Preis geb. 28,— DM. 


Elektronische Datenverarbeitung. (Fachbe- 
richte über programmgesteuerte Maschinen und ihre 
Anwendung, Folge 1). Redaktion Dr. H.K. Schuff. 
IIT+55S. m. Abb. Braunschweig 1959. Friedr. 
Vieweg & Sohn. Preis brosch. 6,80 DM. 


Dr. phil. H. Schlichting und Dr.-Ing. E. Trucken- 
brodt, Aerodynamik des Flugzeuges. Erster 
Band: Grundlagen aus der Strömungsmechanik. 
Aerodynamik des Tragflügels (Teill). XV + 4558. 
m. 260 Abb. Berlin/Göttingen/Heidelberg 1959. 
Springer-Verlag. Preis geb. 52,50 DM. 


Dr. A. Vogel, Vierstellige Funktionentafeln. 
VIII+1578. m. 28. Differenzentafeln u. 248. 
Mathematische Formelsammlung. Stuttgart 1958. 
Verlag Konrad Wittwer. Preis geb. 6,80 DM. 
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M. Deuring, Die Klassenkörper der kom- 
plexen Multiplikation. (Enzyklopädie der mathe- 
matischen Wissenschaften, Band I,, Heft 10, Teil II. 
Herausgegeben im Auftrage der Akademien der Wis- 
senschaften zu Berlin, Göttingen, Heidelberg, Leipzig, 
München und Wien sowie unter Mitwirkung zahl- 
reicher Fachgenossen). 608. Stuttgart 1958. B.G. 
Teubner Verlagsgesellschaft. Preis brosch. 15,— DM. 


W.. Specht, Algebraische Gleichungen mit 
reellen oder komplexen Koeffizienten. (Enzy- 
klopädieder mathematischen Wissenschaften, Band I,, 
Heft 3, TeilII. Herausgegeben im Auftrage der 
Akademien der Wissenschaften zu Berlin, Göttingen, 
Heidelberg, Leipzig, München und Wien sowie unter 


Nachrichten 
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Mitwirkung zahlreicher Fachgenossen). 808. Stutt- 
gart 1958. B.G. Teubner Verlagsgesellschaft. Preis 


brosch. 19,— DM. : 


R.E. D. Bishop und D. 6. Johnson, Schwingungs- 
technische Tabellen. Tafeln zur Berechnung von 
Schwingungen. Deutsche Übersetzung von Prof. Dr. 
H.J. Menges. VIII + 608. Baden-Baden 1958. 
Verlag für angewandte Wissenschaften. Preis geb. 
14,80 DM. 


oa 

S.W. MeCuskey, Introduction to Advanced 
Dynamics. VIIL + 263 8. m. Abb. Reading, Mass. 
1959. Addison-Wesley Publishing Company, Inc. 
Preis geb. 8,50 $. 


NACHRICHTEN 


Am 30. November 1958 verstarb der ungarische 
Mathematiker Prof. Dr. Eugen Egervary. Er wirkte 
bis zu seinem Tode als Professor an der Technischen 
Hochschule Budapest und war Mitglied der Akade- 
mie Müszaki Egyetem Budapest. Die Fachwelt ver- 
liert in ihm einen hervorragenden Vertreter der ange- 
wandten Mathematik und Mechanik. 


Der Zehnte Internationale Kongreß für An- 
gewandte Mechanik findet vom 31. August bis 
7. September 1960 in Stresa (Italien) statt. Neben 
einer Anzahl von Hauptvorträgen auf Grund beson- 
derer Einladungen sind Arbeitssitzungen in zwei Sek- 
tionen Dynamik der Flüssigkeiten (Hydro- und 
Aeromechanik) und Festkörpermechanik (Dy- 
namik der starren Körper, Schwingungen; Elastizität, 
Plastizität und Strukturtheorie) geplant. Die Thermo- 
dynamik soll nicht, Rechenverfahren nur so weit be- 
rücksichtigt werden, als sie sich auf einschlägige Pro- 
bleme beziehen. Um eine ausreichende Diskussion zu 
ermöglichen, wird für jede einzelne Mitteilung eine 
Gesamtzeit von 45 Minuten (30 für den Vortrag und 
15 für die Diskussion) vorgesehen. Dies erfordert, 
daß Auszüge aus den beabsichtigten Vorträgen bis 
zum 1. Januar 1960 an 


The Secretary oftheInternationalCommittee 
Prof.-Mekel-Weg 2, Delft, Netherlands 


in vierfacher Ausfertigung eingereicht werden (zwei 
bis höchstens vier, zweizeilig geschriebene Schreib- 
maschinenseiten, möglichst in zwei der offiziellen 
Kongreßsprachen Deutsch, Englisch, Französisch, 
Italienisch). Die Entscheidung über die Auswahl der 
Vorträge wird von einem Programm-Ausschuß ge- 
troffen und ist endgültig. 

Anschrift des Italienischen Organisations-Komitees 
(Präsident: Prof. G. Colonnetti, Sekretär: Dr. F. Rol- 
la): 

Consiglio Nazionale delle Ricerche, 

Ufficio relazioni internazionali, 

Piazza della Scienza 7, Roma. 


Das Exekutiv-Komitee des Kongresses setzt sich zu- 
sammen aus den Herren C. B. Biezeno, Präsident; 
Richard V. Southwell, W. T. Koiter (Sekretär). 


Fourth Congress on Theoretical and Applied Mechanies. 

The Fourth Congress on Theoretical and Applied 
Mechanics was held at the Bengal Engineering College, 
Howrah from December 28 to December 31, 1958, 
under the Presidentship of Dr. $S. R. Sen Gupta, 
Director, Indian Institute of Technology, Kharagpur. 
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About two hundred and fifty scientists and engineers 
registered themselves for the Congress. These inclu- 
ded members from Australia, Burma, Czechoslova- 
kia, Egypt, Hungary, Italy, Japan, Poland, U.S.A. 
and the U.S.S.R. Sri A.C. Roy, Principal, Bengal 
Engineering College, Howrah welcomed the delegates. 
Messages were read from prominent workers all over 
the world including those from von Karman, H.L. 
Dryden, S. Goldstein, L. Rosenhead, K.S. Krishnan, 
J.C. Ghosh and M. S. Tihacker. 

The Congress received 75 original papers of which 
40 were read. The subjects dealt with included Finite 
deformation, Visco-elastieity, Stress waves, Stresses in 
strips, columns and discs, Plasticity, Elasto-porous 
problems, Vibration and Stability, Fluid flow, Ballis- 
tics and Statistics. 


In his Presidential Address Dr. S. R. Sen Gupta 
stressed the importance of experimental methods in 
the solution of engineering problems. He considered 
in detail the causes of cracking of some capitals of 
stone columns in the Arts Faculty building of the 
Osmania University. By rejecting a few of the plau- 
sible hypotheses he concluded that the cracks might 
be due to the localized concentration of stress due to 
unevenness of bearing between the capital and the 
column. 

Half hour addresses were delivered by S. K. Chakra- 
varty, 8. I. Paiand M. Lunc. Dr. Chakravarty spoke 
on propagation of elastic waves across continents and 
oceans. He dwelt on the existence of small distur- 
bances, called microseisms superposed on Rayleigh 
waves through multi-layered media. Prof. Pai dis- 
cussed the propagation of a cylindrical shock wave 
produced by the instantaneous release of energy from 
infinite wire into the surrounding medium. Prof. Lunc 
outlined the molecular aspects of gas flows. 


Popular lectures on “The present concept of the 
universe” and “Mechanics and the engineer” were’ 
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respectively delivered by Prof. N.R. Sen and Dr. » 


J. C. Morrison of Glasgow University. 


Dr. A.N. Khosla was elected President for the 
next two years. Dr. S.K. Chakravarty was elected 
Vice-president in the vacancy created by the retire- 
ment of Sri V. Cadambe. Prof. B.R. Seth was re- 
elected Secretary-Treasurer for the next three years. 
Dr. Jai Kishan, Prof. B. Sen Gepta, and Prof. D. 
Banerjee were elected new members of the Executive 
Committee. 

The Executive Committee accepted the invitation 
of University of Roorkee to hold the Fifth Congress in 
Roorkee in December 1959. 
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